Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 755,400 G→A 100% W75* (TGG→TAG)  sdhC → succinate dehydrogenase, membrane subunit, binds cytochrome b556

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913755,4000GA98.7% 252.1 / ‑4.3 74W75* (TGG→TAG) sdhCsuccinate dehydrogenase, membrane subunit, binds cytochrome b556
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (41/32);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCTCCATCGCGTTTCCGGTGTGATCACCTTTGTTGCAGTGGGCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTGGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGCTTTCACTT  >  NC_000913/755260‑755542
                                                                                                                                            |                                                                                                                                              
tctcCATCGCGTTTCCGGTGTGATCACCTTTGTTGCAGTGGGCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATc                                                                                                                                         >  8:148245/1‑148 (MQ=255)
 ctcCATCGCGTTTCCGGTGTGATCACCTTTGTTGCAGTGGGCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATc                                                                                                                                         <  7:148245/147‑1 (MQ=255)
              ccGGTGTGATCACCTTTGTTGCAGTGGGCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTgg                                                                                                                            <  4:224142/147‑1 (MQ=255)
                     gATCACCTTTGTTGCAGTGGGCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCAc                                                                                                                    >  7:435735/1‑148 (MQ=255)
                           cTTTGTTGCAGTGGGCATCCTGCTGTGGCTACTGGGTACCAGCCTCTCTTCCCCGGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGACAAATTTATCATGTAGGGCCTCCTTACCGCTCTGGCGTATCAcgtcgtc                                                                                                              >  4:52312/1‑148 (MQ=255)
                               gttgCAGTGGGCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCAcgtcgtcgca                                                                                                           <  1:7647/147‑3 (MQ=255)
                                  gCAGTGGGCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGt                                                                                                        <  8:435735/147‑1 (MQ=255)
                                      tGGGCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTGGGGCATCCTTACCGCTCTTGCGTATCACGTCGTCGTAGGTATTCg                                                                                                   >  1:262600/1‑148 (MQ=255)
                                        ggCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGc                                                                                                  >  6:309715/1‑147 (MQ=255)
                                         gCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCGTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCca                                                                                                <  5:309715/148‑1 (MQ=255)
                                             ccTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACat                                                                                             >  5:93325/1‑147 (MQ=255)
                                                 ctgTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACatgatg                                                                                         <  1:294990/147‑1 (MQ=255)
                                                   gtgGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCAcg                                                                                                                   >  8:109865/1‑119 (MQ=255)
                                                   gtgGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCAcg                                                                                                                   <  7:109865/119‑1 (MQ=255)
                                                    tgGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACatgatgatg                                                                                      >  7:337231/1‑147 (MQ=255)
                                                         tCTGGGTAACAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATATATCATGTAGGTCATCCTTAAAGCTCTGGCGTATCACGTCGTCGTAGGTATACGCCACATGAAGATCGAtttt                                                                                >  1:178705/1‑148 (MQ=255)
                                                           tGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTgg                                                                              >  5:112118/1‑148 (MQ=255)
                                                            gggTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCcaca                                                                                              <  5:301601/131‑1 (MQ=255)
                                                            gggTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCcaaa                                                                                              >  6:301601/1‑129 (MQ=255)
                                                             ggTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCt                                                                            <  1:59034/148‑1 (MQ=255)
                                                              gTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTa                                                                                                                                     >  6:387657/1‑90 (MQ=255)
                                                              gTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTa                                                                                                                                     <  5:387657/90‑1 (MQ=255)
                                                              gTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCg                                                                                                                          >  4:58382/1‑101 (MQ=255)
                                                              gTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCg                                                                                                                          <  3:58382/101‑1 (MQ=255)
                                                              gTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTa                                                                           >  2:165427/1‑148 (MQ=255)
                                                                  cAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATc                                                                                                                      >  2:231742/1‑101 (MQ=255)
                                                                  cAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATc                                                                                                                      <  1:231742/101‑1 (MQ=255)
                                                                    gCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACatg                                                                                            <  1:67633/125‑1 (MQ=255)
                                                                    gCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACatg                                                                                            >  2:67633/1‑125 (MQ=255)
                                                                     ccTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCAACCTTACCGCTCTGGCGTATCAAGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGATATCTGgat                                                                    >  4:367155/1‑147 (MQ=255)
                                                                       tctctTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGAtt                                                                                  >  3:176558/1‑132 (MQ=255)
                                                                       tctctTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGAtt                                                                                  <  4:176558/132‑1 (MQ=255)
                                                                        ctctTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGgaaga                                                                  >  4:344863/1‑147 (MQ=255)
                                                                         tctTCCCCTGAAGGTTTCGAGCAAGATTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATACTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGgaagaa                                                                 >  1:266051/1‑147 (MQ=255)
                                                                          ctTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTgg                                                                              <  3:195301/133‑1 (MQ=255)
                                                                          ctTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTgg                                                                              >  4:195301/1‑133 (MQ=255)
                                                                            tCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGAGGATGGATTTTGGCTATCTGGAAGAAACAt                                                             >  6:379496/1‑148 (MQ=255)
                                                                             ccccTGAAGGTTTCGAGCAATCTTCCGCTATTATTGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACAt                                                             <  2:178705/147‑1 (MQ=255)
                                                                             ccccTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACAt                                                             <  5:405428/147‑1 (MQ=255)
                                                                             ccccTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACAt                                                             <  6:112118/147‑1 (MQ=255)
                                                                              cccTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTgg                                                                                                                            >  7:481/1‑83 (MQ=255)
                                                                              cccTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTgg                                                                                                                            <  8:481/83‑1 (MQ=255)
                                                                              cacTGAAGGTTTCGAGCAAGCTTCGGCGATTATGGGCATCTTCTTCGTCAAATTTTTCATGTAGGGCATCCTTACCGCTCTGGCGTATCAGGTCGTCGTAGGTATTCGCCACATGATGATGGAttt                                                                                 <  6:78589/124‑1 (MQ=255)
                                                                               ccTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGAtttt                                                                                >  5:78589/1‑126 (MQ=255)
                                                                                 tGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACatgatga                                                                                        <  2:180469/116‑1 (MQ=255)
                                                                                 tGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACatgatga                                                                                        >  1:180469/1‑116 (MQ=255)
                                                                                   aaGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGctct                                                                                                                              >  3:120528/1‑76 (MQ=255)
                                                                                   aaGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGctct                                                                                                                              <  4:120528/76‑1 (MQ=255)
                                                                                        ttCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGt                                                  >  4:364789/1‑147 (MQ=255)
                                                                                         tCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCATTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTa                                                 >  3:248483/1‑147 (MQ=255)
                                                                                          cGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTg                                                                       >  1:91084/1‑124 (MQ=255)
                                                                                          cGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTg                                                                       <  2:91084/124‑1 (MQ=255)
                                                                                               aaGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTc                                          >  8:247746/1‑148 (MQ=255)
                                                                                                aGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACatgatgatg                                                                                      >  1:79611/1‑103 (MQ=255)
                                                                                                aGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACatgatgatg                                                                                      <  2:79611/103‑1 (MQ=255)
                                                                                                    tCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGcc                                      >  1:6501/1‑147 (MQ=255)
                                                                                                    tCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGcc                                      >  3:33139/1‑147 (MQ=255)
                                                                                                       gcgATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCaaaa                                  >  2:4372/1‑148 (MQ=255)
                                                                                                         gATTATGGGCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAAtc                                <  3:364789/148‑1 (MQ=255)
                                                                                                               gggCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCtt                           <  4:33139/147‑1 (MQ=255)
                                                                                                                ggCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTg                         >  6:62203/1‑148 (MQ=255)
                                                                                                                 gCAGCTTCTTCGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGt                        <  7:247746/148‑1 (MQ=255)
                                                                                                                        cttcGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTaa                                                <  5:269434/117‑1 (MQ=255)
                                                                                                                        cttcGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTaa                                                >  6:269434/1‑117 (MQ=255)
                                                                                                                        cttcGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCaaa                                   >  3:399123/1‑130 (MQ=255)
                                                                                                                        cttcGTCAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCaaa                                   <  4:399123/130‑1 (MQ=255)
                                                                                                                              cAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTg           >  5:264393/1‑148 (MQ=255)
                                                                                                                              cAAATTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTg           >  2:109784/1‑148 (MQ=255)
                                                                                                                               aaaTTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGc          <  3:344863/148‑1 (MQ=255)
                                                                                                                                 aTTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGCtt        <  1:165427/148‑1 (MQ=255)
                                                                                                                                 aTTTATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGCtt        >  6:339776/1‑148 (MQ=255)
                                                                                                                                    tATCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGCTTTCa     >  4:360063/1‑148 (MQ=255)
                                                                                                                                     aTCATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGCTTTCa     <  8:337231/147‑1 (MQ=255)
                                                                                                                                       cATGTAGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGCTTTCACtt  >  5:222860/1‑148 (MQ=255)
                                                                                                                                         tgtaGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTg                         >  6:209908/1‑123 (MQ=255)
                                                                                                                                         tgtaGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTg                         <  5:209908/123‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
TCTCCATCGCGTTTCCGGTGTGATCACCTTTGTTGCAGTGGGCATCCTGCTGTGGCTTCTGGGTACCAGCCTCTCTTCCCCTGAAGGTTTCGAGCAAGCTTCCGCGATTATGGGCAGCTTCTTCGTCAAATTTATCATGTGGGGCATCCTTACCGCTCTGGCGTATCACGTCGTCGTAGGTATTCGCCACATGATGATGGATTTTGGCTATCTGGAAGAAACATTCGAAGCGGGTAAACGCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGCTTTCACTT  >  NC_000913/755260‑755542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: