Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2750TC50.7% 53.2 / 82.4 67pseudogene (335/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base C (10/24);  total (26/41)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-07

TCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTT  >  NC_000913/2228128‑2228415
                                                                                                                                                         |                                                                                                                                            
tctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTTTCGCGTTTCGCCCCCCCGGCCTTTAACCGCCCCCCGGGGGCCGTAACCTcc                                                                                                                                            >  2:480434/1‑150 (MQ=255)
 ctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTTCCCCCCCCCCGGTTTAAACGTACCCCGGCGGGGAGGAACcccc                                                                                                                                            >  1:1085478/1‑149 (MQ=255)
   ttttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGGTCCGCCGCCCCCGCTTTTACCCCACCCCGGCGGGGAGGAACcccccccc                                                                                                                                        >  2:406211/1‑151 (MQ=255)
     ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCg                                                                                                                                       <  1:607729/150‑1 (MQ=255)
     ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCg                                                                                                                                       <  2:1278638/150‑1 (MQ=255)
     ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCCCCCGGGCTTTTAACCGACCCCGGCGGGCAGGAACCTcccccc                                                                                                                                        >  1:685465/1‑149 (MQ=255)
      tCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTTCGCCCCCCGCGCTTTTAACCCAACCCAGGCGGGGGGAAACCCCCCCCCgt                                                                                                                                      >  2:1242831/1‑150 (MQ=255)
         aTTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGAtt                                                                                                                                    <  2:630909/149‑1 (MQ=255)
         aTTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGAt                                                                                                                                     <  2:404564/148‑1 (MQ=255)
              ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTTGCCGCGGTTCCCCGCCCCCGCCTTTAACCTACCCCAGCGGGGAGGAAACCCCCCCCCGCTTCTTccc                                                                                                                              >  2:140088/1‑151 (MQ=255)
                gcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTcccc                                                                                                                            <  1:969332/150‑1 (MQ=255)
                  aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTAAGCTGCCAGCCAGCAATTGGCGGGCTTTTGCTCGTTTCGCCGCCCCCGCCTTTACCCGACCCCCGGGGGCAGTTACCCCAGCCCGATTCtttttcc                                                                                                                            >  2:628036/1‑144 (MQ=255)
                       cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCg                                                                                                                    <  2:226844/151‑1 (MQ=255)
                       cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTTTCCCGTTTCGCCCCCCCTGCTTTTACCCTACCCCGGGGGGGCGTAACCCCCCCCCCCTTTTTTTCCCACatat                                                                                                                      >  2:916588/1‑149 (MQ=255)
                        gtgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCga                                                                                                                   <  2:374222/151‑1 (MQ=255)
                              aaCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGc                                                                                                              <  1:1210028/150‑1 (MQ=255)
                                cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAg                                                                                                           <  1:691275/151‑1 (MQ=255)
                                cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCCCGCGTTCCCCCGCCCGCGGCTTTACCGCACCCCCGCGCGGCGTAAACCCCGCCCCCTTTTTTCCCCCGATATCGGGATGcc                                                                                                             >  2:135789/1‑149 (MQ=255)
                                  ccgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTTCGCCGCCCCTGCTTTTACCGGACCCCGGCGGGGAGGAAACCCCCCCCGCTTTTTTTCCCCGATTTCGAGATGgccgcc                                                                                                         >  2:405183/1‑151 (MQ=255)
                                   cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                        <  1:971808/151‑1 (MQ=255)
                                   cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGGGTTCCGCCGCCCGTGCTTTATCCGTACCCAGGCGGGGAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                        <  1:1102931/151‑1 (MQ=255)
                                      cTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCGCCGCCCGGGCTTTTACCGTACCCAGGCGGGGAGTAACCCCCCCCCGATTCTTTTCCCCGATATCGAGATGCCAGcccccc                                                                                                     >  1:952141/1‑151 (MQ=255)
                                             tGCCACTCCCAGTACAGGCGGGGTACGCTGCCAGCCAGCAATTTGCCGGTTTGCTCCCGCTCCCCCGCCCCGGCGTTTACACTGCCCCGGGGGGCAGGAACCCGCCCCCCCTTCTTTCTCCCCATATCGCGGTGCCACCccgccccgaag                                                                                               >  2:457709/1‑144 (MQ=255)
                                                           cAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCGCCCGTGCCTTAACCCTACCCCGGCGGGCCGGAACCCCCCCCCCGCTCTTTTCCCAGATATTGCGATGGCACGCCCCCCTTAAGGCCAAAggaac                                                                                    >  2:426519/1‑143 (MQ=255)
                                                           cAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAAACGTACCCAGGCGGGCAGTAACCTCCCCCCGATTCTTTCCCCAGATATCGAGATGGCAGCCCGCCGTTAAGCCCAAAGAACCGtt                                                                                 >  2:653400/1‑150 (MQ=255)
                                                              gCGGGCTCCGCTGCCGGCCAGCAATTGGCGGGTTTGCTCGCGTCCCGCCCCCCGTTCTTTACCCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGgtg                                                                             <  2:11938/151‑1 (MQ=255)
                                                                  gCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCCCCCGCTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGCTGGTGTAcc                                                                         >  2:200727/1‑151 (MQ=255)
                                                                    tACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGAACCGTTGGTGTAc                                                                          >  1:948955/1‑148 (MQ=255)
                                                                      ccctccccgccaccAATTGGCGGGTTTGCTCGCTTCCCCCCGCCCGTGCTTTAACCGTCCCCAGGGGGGCAGTCACCCCCGCCCGATTTTTCCCCCATATCTCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTAcc                                                                         <  2:170177/139‑1 (MQ=255)
                                                                      cGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCGCCGCCCGTGCTTTAAACGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGGACCGTTGGTGTACCACgg                                                                     >  1:58106/1‑151 (MQ=255)
                                                                      cGCTGCCACCCAGCAATGGGGGGGTTTGCTCGCTTCCCCCCGCCCGTTTTTTACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGTTTCTTCCCCCAGATTTCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCa                                                                        <  2:808163/148‑1 (MQ=255)
                                                                       gCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGCGCTTTAACCGTACCCAGGCGGGCAGTAACCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCCTTGGTGTACCACggg                                                                    >  2:934896/1‑151 (MQ=255)
                                                                        cTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCCCCGCCCGTGCTTTAACCCTAACCAGGCGGGCAGTAAACTCCGCCCGAATCTTTTCCCCGATATCCAGATGCCAGCCCGCCCTTAAGCCCAAAGTACCGTTGGTGTACCACcgg                                                                    >  2:782351/1‑150 (MQ=255)
                                                                           ccagccagCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTACCCGTCCCCGGGCGGGCAGTACCCCCCGCCCTATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                  <  1:344561/149‑1 (MQ=255)
                                                                            cggccgacAATTGGCGGTTTTCTTCGCGTTCCCCGGCCATGTCTTTAACCTTCCCCTGGAGTGCGATTACTCCGGCTCGATTATCGCGCAATTTCTCGAGACGCAGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCg                                                                 <  2:911859/142‑1 (MQ=255)
                                                                             agcccccAATTGGCGGGTTTGCTCGCGTCCCGCCGCCGGTGGTTTAACGGCACCCCGCCGGTAGGTACCACCCGCCGAATTCTTGCCCCAAATACGGAGATGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGc                                                                   <  1:311383/140‑1 (MQ=255)
                                                                               caagcAATTGGCGGTTTTGCGCGCTTTCCCCCGCCCGTGCTTAACCCGTACCCGGGCGGGTATTACCCCCCGCCCGTTTTTTTCCCCAGATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGc                                                                   <  2:937582/142‑1 (MQ=255)
                                                                                      tgggCGGTTTTGCTCGCGTTCCGCCGCCCGTGTTTTACCCGTCCCCAGGCGGGCAGTAACCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGc                                                      <  1:128426/148‑1 (MQ=255)
                                                                                       tGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCgg                                                    >  1:995675/1‑151 (MQ=255)
                                                                                        ggCGGGTTGCCTCGGGTTCCCCGGCCGGTCCTTTACCCTTACCCAGGGGGGCGGTACCCCCCGCCGGTTTTTTCCCCCGGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCg                                                     <  1:492526/149‑1 (MQ=255)
                                                                                         gggggttttGCTCGCGTCCCGCCGCCCGTGCTTTACCCGTCCCCGGGCGGGCAGTACCCCCCGCCCGTTTCTTCCCCCAGATATCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCgg                                                    <  2:165186/144‑1 (MQ=255)
                                                                                         gCGGGTTCGCCCGCGTCCCCCCGCCCGTGTTTTACCCGTCCCAGGGGGGGCAGTACCCCCCGCCCGTTTCTTTCCCCAGATAACGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGcc                                                  <  2:989269/151‑1 (MQ=255)
                                                                                          cGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCg                                                 >  2:711665/1‑151 (MQ=255)
                                                                                            gggTTTGCCCGGTTCCCGCCGCCCGTGCTTTACCCGTCCCCGGGCGGGCAGAACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCgg                                                <  1:197255/151‑1 (MQ=255)
                                                                                             gttttGCTCGCGTCCCCCCGCCCGTGTTTTACCCGTCCCAGGGCGGACAGAACCCCCCGCCCTTTTTTTCCCCCAGATATCGGGGTCCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCgg                                                <  2:191536/148‑1 (MQ=255)
                                                                                               tttGCTCGCTTCCCGCCCCCCGTGTTTAACCCGTACCCGGGCGGGAAGTAACCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGAt                                              <  1:753510/150‑1 (MQ=255)
                                                                                               tttGCTCGCGTTCCGCCGCCCGTGCTTTATCCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGCGATGCCAGCCCGCCGTTATCCCAAAAGTACCGATGACGTCCCTCGCGCCGGTCGTACCTGCGGCCGGa                                               >  2:975839/1‑149 (MQ=255)
                                                                                                   cccgcgTCCCCCCGCCTTTCTTTTAACCGTCCCAGGGCGGGCAGCACCCTCCCCCCGATTCTTGCCCAAGATATGGATATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATc                                             <  2:790102/145‑1 (MQ=255)
                                                                                                      gctttCCCCCCCCTGTGTTTTACCCGTCCCCGGGCGGGTAGTACCCCCCGCCCGATTTTTCCCCCAGATATGGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGt                                           <  1:369751/143‑1 (MQ=255)
                                                                                                      gcgTTCCGCCGCCGGGGTGTAACCCTTCCCGCGGCGGGAGGTACCCCCCGCCGGATTTTTTCCCCAGATAACGAGACGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGtt                                          <  1:1248511/147‑1 (MQ=255)
                                                                                                        gTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                    >  2:1267427/1‑151 (MQ=255)
                                                                                                           ccgccgCCCGTGTTTTAACCGTCCCCAGGCGGGCATTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                    <  2:1251360/148‑1 (MQ=255)
                                                                                                           ccgccgCCCGTGCTTTAACCGTCCCCGGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCa                                   <  2:203447/149‑1 (MQ=255)
                                                                                                           ccgccgCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCaa                                  >  1:344259/1‑150 (MQ=255)
                                                                                                               ccgccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACgg                              >  2:960347/1‑151 (MQ=255)
                                                                                                                  cccGTGTTTTACCCGTCCCCGGGCGGGCAGTACCCCCCCCCCGTTTTTTCCCCCGGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAAc                                <  2:683883/146‑1 (MQ=255)
                                                                                                                  cccGTCCTTTAGCCGTCCCGGGGCGGAAAGTGACCCCCGCCCGTTTTTTCCCCCAAATGTCGAAAGCCCAGGCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATc                                             <  1:1165899/133‑1 (MQ=255)
                                                                                                                      tGCTTAAACCGTACCGGGGCGGGGAGTAACCTCCGCCCGATTTTTCCCCCAAATACAAAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCTCCAt                        <  1:879550/150‑1 (MQ=255)
                                                                                                                            ttAACCGTCCCGAGGAGGGCATTACCCCCCGCCAGATTTTTTCCCAAGATATCGAGGTGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAg                    <  1:1226809/150‑1 (MQ=255)
                                                                                                                              accctaaCCCGGGGGGGGATTACCCCCCGCCCGATTCTTTCCCAAAATAGCGAGATGCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTa                      <  1:295523/140‑1 (MQ=255)
                                                                                                                              acccGTCCCCAGGGGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGcc                  <  1:1257026/148‑1 (MQ=255)
                                                                                                                                cggtccccAGGGGGGAAGTACCCCCCGCCCGTTTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCtt                <  1:533625/145‑1 (MQ=255)
                                                                                                                                cctcaCCCGGGCGGGCGATACCCCCCCCCCGATTCTTCCCCCAGATAGCGGGAGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCtt                <  2:304077/146‑1 (MQ=255)
                                                                                                                                ccGTCCCCAGGCGGGCAGAACCCTCCCCCCGTTTCTTTCCCCAGATACCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGcc                  <  1:480454/148‑1 (MQ=255)
                                                                                                                                 cGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTc               >  1:579096/1‑150 (MQ=255)
                                                                                                                                  gTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg             >  1:829407/1‑151 (MQ=255)
                                                                                                                                  gTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg             >  2:977701/1‑151 (MQ=255)
                                                                                                                                             gggCATTACCCCCCCCCCGTTTTTTTCCCCGGATATGGAGATGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt    <  1:1012241/150‑1 (MQ=255)
                                                                                                                                              ggtaGTACCCCCCGCCCGTTTCTTTCCCGACATATGGGAAGCCCATCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttttt  <  1:628060/148‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                            
TCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTT  >  NC_000913/2228128‑2228415

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: