Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,491,0980AC21.8% 259.8 / 20.3 124Y66S (TAC→TCC) cybBcytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (66/31);  new base C (0/27);  total (66/58)
Fisher's exact test for biased strand distribution p-value = 1.83e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 133 total aligned reads displayed.

TATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCG  >  NC_000913/1490961‑1491234
                                                                                                                                         |                                                                                                                                        
tATCGCAGCGTATTGCGCAATGGAGTTTCGTGTTTTCTCCCCCCGTAGCGATCGGCACCTCATCAACATGATTCATTTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTTTTTTGTTAAGGCTAAAATccc                                                                                                                                        <  2:592629/140‑1 (MQ=255)
 aTCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCAACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATAccc                                                                                                                                       >  1:592629/1‑140 (MQ=255)
       gCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTTTTAAGGCTGAAATACCCACcccc                                                                                                                                 <  2:1594039/140‑1 (MQ=255)
        cGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCccg                                                                                                                                <  1:103322/140‑1 (MQ=255)
        cGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCccg                                                                                                                                >  1:1326864/1‑140 (MQ=255)
         gTATTGCGCAATGGAGTTTCGTGTTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCccgc                                                                                                                               <  1:655753/140‑1 (MQ=255)
         gTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCccgc                                                                                                                               <  2:1422344/140‑1 (MQ=255)
         gTATTGCGCAATGGAGTTTCGTGGTTTCTCCCCACGTGGCGATCGGCCCCTCATAAACATGATTCATGTTTCCTGGGGCATCTCAATTCTCGTGCTGATGGTCGTTCTTCTTTTGTTAAGGCTGAAATCCCCACCCccgc                                                                                                                               <  2:1326864/140‑1 (MQ=255)
             tGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGAt                                                                                                                           >  2:1549732/1‑140 (MQ=255)
              gcgcAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATAGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGAtt                                                                                                                          <  1:43702/140‑1 (MQ=255)
               cgcAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCccg                                                                                                                                >  2:293535/1‑133 (MQ=255)
               cgcAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCccg                                                                                                                                <  1:293535/133‑1 (MQ=255)
                    tGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCt                                                                                                                    >  1:253978/1‑140 (MQ=255)
                              tGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAAc                                                                                                          >  1:324961/1‑140 (MQ=255)
                                  ttcttcCCGCGAAGCGATGGGCCACTCATCAACTTGTTTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTTTAAGGGCGAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCgat                                                                                                      <  2:241502/140‑1 (MQ=255)
                                  ttcttcCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCCCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCgat                                                                                                      <  2:1731628/140‑1 (MQ=255)
                                       cccACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTTTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGAtgatga                                                                                                 <  1:1729238/140‑1 (MQ=255)
                                          aCGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACgg                                                                                              <  1:194879/140‑1 (MQ=255)
                                           cGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCgat                                                                                                      <  2:1291353/131‑1 (MQ=255)
                                           cGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCgat                                                                                                      >  1:1291353/1‑131 (MQ=255)
                                           cGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACggg                                                                                             >  2:811603/1‑140 (MQ=255)
                                               gCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTg                                                                                         >  1:17634/1‑140 (MQ=255)
                                               gCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTg                                                                                         >  1:1870046/1‑140 (MQ=255)
                                                     ggCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCAt                                                                                   <  1:1753491/140‑1 (MQ=255)
                                                      gCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCAtt                                                                                  >  2:756065/1‑140 (MQ=255)
                                                      gCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCAtt                                                                                  >  2:1498175/1‑140 (MQ=255)
                                                        cACTCATCAACATGATTCATGTTTCCTGTGGGATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTg                                                                                >  1:830495/1‑140 (MQ=255)
                                                        cACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTg                                                                                <  1:1526390/140‑1 (MQ=255)
                                                        cACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTg                                                                                >  2:410363/1‑140 (MQ=255)
                                                          cTCATCAAAATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTggg                                                                              >  2:1513594/1‑140 (MQ=255)
                                                           tcatcaACATGATTCTTTTTTCCTGTGGCATCTAAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGa                                                                             <  2:533033/140‑1 (MQ=255)
                                                             atcaACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGgc                                                                                       >  2:1547436/1‑128 (MQ=255)
                                                             atcaACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGgc                                                                                       <  1:1547436/128‑1 (MQ=255)
                                                             atcaACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACa                                                                           >  1:1865241/1‑140 (MQ=255)
                                                             atcaACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACa                                                                           >  2:16916/1‑140 (MQ=255)
                                                             atcaACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACa                                                                           >  1:511982/1‑140 (MQ=255)
                                                               caACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGActgt                                                                         >  1:1445007/1‑137 (MQ=255)
                                                               caACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACAtt                                                                         >  2:1692468/1‑140 (MQ=255)
                                                                aaCATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACAttt                                                                        <  1:757663/140‑1 (MQ=255)
                                                                 aCATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTg                                                                       >  2:765647/1‑140 (MQ=255)
                                                                  cATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTgg                                                                      <  2:56715/140‑1 (MQ=255)
                                                                   aTGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCTTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGt                                                                     <  1:335150/140‑1 (MQ=255)
                                                                     gATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGa                                                                   >  1:1777790/1‑140 (MQ=255)
                                                                       ttCATGTTTCCTGTGGCATCTCAATTCTCTTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGAtt                                                                 >  2:461538/1‑140 (MQ=255)
                                                                       ttCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGt                                                                     <  1:1058545/136‑1 (MQ=255)
                                                                       ttCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGt                                                                     >  2:1058545/1‑136 (MQ=255)
                                                                        tCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGAttt                                                                <  1:135437/140‑1 (MQ=255)
                                                                          aTGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAt                                                              <  1:1312060/140‑1 (MQ=255)
                                                                          aTGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAt                                                              <  1:1487578/140‑1 (MQ=255)
                                                                           tGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGTGACTGGCGCATTTGGGACATTTGGt                                                                     >  1:492127/1‑132 (MQ=255)
                                                                            gTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAtct                                                            >  2:425672/1‑140 (MQ=255)
                                                                                 cTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGt                                                       >  1:87169/1‑140 (MQ=255)
                                                                                  tgtgGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGtt                                                      >  1:1368961/1‑140 (MQ=255)
                                                                                    tgGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTa                                                    >  1:787803/1‑140 (MQ=255)
                                                                                       cATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTg                                                 >  2:602852/1‑140 (MQ=255)
                                                                                         tctcAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAt                                                              <  2:432829/125‑1 (MQ=255)
                                                                                         tctcAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAt                                                              >  1:432829/1‑125 (MQ=255)
                                                                                         tctcAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCt                                                         <  1:2868/130‑1 (MQ=255)
                                                                                         tctcAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCt                                                         >  2:2868/1‑130 (MQ=255)
                                                                                          ctcAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTgcgc                                              >  1:590706/1‑140 (MQ=255)
                                                                                            cAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTg                                            <  2:1119656/140‑1 (MQ=255)
                                                                                             aaTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGTACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                           <  2:703710/140‑1 (MQ=255)
                                                                                             aaTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                           <  1:883462/140‑1 (MQ=255)
                                                                                              aTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGcc                                          <  1:1501999/140‑1 (MQ=255)
                                                                                               ttCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCa                                         >  2:1245348/1‑140 (MQ=255)
                                                                                                   cGTGCTGATGGTCTTTCGTTTTTTGTAAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGa                                     <  2:455405/140‑1 (MQ=255)
                                                                                                   cGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGa                                     >  1:1644393/1‑140 (MQ=255)
                                                                                                    gTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                                    >  1:455405/1‑140 (MQ=255)
                                                                                                     tGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                           <  2:1194647/132‑1 (MQ=255)
                                                                                                     tGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                           >  1:1194647/1‑132 (MQ=255)
                                                                                                     tGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAtt                                   >  2:1214020/1‑140 (MQ=255)
                                                                                                       cTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                                    >  1:1470346/1‑137 (MQ=255)
                                                                                                       cTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                                    <  2:1470346/137‑1 (MQ=255)
                                                                                                        tGATGGTCGTTCGTTTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggt                                <  2:1865241/140‑1 (MQ=255)
                                                                                                         gATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGTCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggtt                               >  2:1296194/1‑140 (MQ=255)
                                                                                                          aTGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggttt                              >  2:1764642/1‑140 (MQ=255)
                                                                                                            ggTCGTTCGTCTTTTGTAAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggtttgg                            <  2:511982/140‑1 (MQ=255)
                                                                                                              tcgttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtg                          >  1:22857/1‑140 (MQ=255)
                                                                                                               cgttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                         <  2:299933/140‑1 (MQ=255)
                                                                                                               cgttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                         <  1:16916/140‑1 (MQ=255)
                                                                                                               cgttcgtCTTTTGTTAAGGCTAAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                         <  1:425672/140‑1 (MQ=255)
                                                                                                                tttcgtTTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgat                        <  2:22857/139‑1 (MQ=255)
                                                                                                                gttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgat                        <  1:811603/140‑1 (MQ=255)
                                                                                                                 ttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCa                                         <  2:1108006/122‑1 (MQ=255)
                                                                                                                 ttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCa                                         >  1:1108006/1‑122 (MQ=255)
                                                                                                                 ttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgatg                       <  1:1613039/140‑1 (MQ=255)
                                                                                                                 ttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgatg                       >  1:1292314/1‑140 (MQ=255)
                                                                                                                   cgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTgatgat                     >  1:1899443/1‑140 (MQ=255)
                                                                                                                    gtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGatgatg                    <  1:1549732/140‑1 (MQ=255)
                                                                                                                     tCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGt                   <  1:300821/140‑1 (MQ=255)
                                                                                                                          tGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAAc              <  1:1686542/140‑1 (MQ=255)
                                                                                                                           gTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAAcc             <  2:1292314/140‑1 (MQ=255)
                                                                                                                              aaGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCggg          <  2:1385620/140‑1 (MQ=255)
                                                                                                                               aGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCgggg         >  1:1160752/1‑140 (MQ=255)
                                                                                                                               aGGCTGAAAACCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCgggg         <  1:1513594/140‑1 (MQ=255)
                                                                                                                                ggCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGc        <  2:1899443/140‑1 (MQ=255)
                                                                                                                                ggCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGc        <  2:17634/140‑1 (MQ=255)
                                                                                                                                 gCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCa       <  1:1220303/140‑1 (MQ=255)
                                                                                                                                      aaTACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCg  >  2:491880/1‑140 (MQ=255)
                                                                                                                                      aaTACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCg  >  1:1090919/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                        
TATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCG  >  NC_000913/1490961‑1491234

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: