Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6570AG41.2% 111.5 / 21.6 80F386S (TTC→TCC) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (18/29);  new base G (33/0);  total (51/29)
Fisher's exact test for biased strand distribution p-value = 1.36e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCG  >  NC_000913/3508514‑3508804
                                                                                                                                               |                                                                                                                                                    
aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGgaaaaat                                                                                                                                                >  2:884403/1‑148 (MQ=255)
 cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTc                                                                                                                                               <  1:298566/150‑1 (MQ=255)
 cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGCGGGAAAATTc                                                                                                                                               >  2:1058084/1‑150 (MQ=255)
  gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATca                                                                                                                                              >  2:800587/1‑150 (MQ=255)
    tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGGGGGAAAAATcaaa                                                                                                                                            >  2:916272/1‑148 (MQ=255)
    tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcaca                                                                                                                                            <  2:876995/150‑1 (MQ=255)
     gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCTGCGCGCGGGAAAATTcaaac                                                                                                                                           >  1:886712/1‑150 (MQ=255)
       ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATtacaca                                                                                                                                          >  1:348385/1‑149 (MQ=255)
       ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTacacac                                                                                                                                         >  1:18479/1‑150 (MQ=255)
           aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTacacaccg                                                                                                                                       >  1:1172928/1‑146 (MQ=255)
           aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTAATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGTCTGGCGCaaaaaatatcaccccgca                                                                                                                                     >  1:870848/1‑137 (MQ=255)
            aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGc                                                                                                                                   <  2:337192/151‑1 (MQ=255)
            aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTa                                                                                                                                     >  2:1245314/1‑149 (MQ=255)
                 ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCca                                                                                                                              <  1:1238767/151‑1 (MQ=255)
                  gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGCCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCcac                                                                                                                             <  1:155835/151‑1 (MQ=255)
                    gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACCGGAACCCCccca                                                                                                                            >  1:953273/1‑150 (MQ=255)
                      gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                         <  1:704003/151‑1 (MQ=255)
                      gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                         <  1:340566/151‑1 (MQ=255)
                         cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGaa                                                                                                                      <  1:264365/151‑1 (MQ=255)
                          aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGaaa                                                                                                                     <  2:1082335/151‑1 (MQ=255)
                          aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACCGGAGCCCCCACCACAGaa                                                                                                                      >  1:1138629/1‑150 (MQ=255)
                                      cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTACACACGGTAGCCCCCACATCAGAAAGCGTGGGgcg                                                                                                           >  2:896795/1‑149 (MQ=255)
                                         aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGTGGGGCGcaca                                                                                                       >  2:1051441/1‑150 (MQ=255)
                                           gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAAAGGTACGCCCCACATCAGaaaagggggggcgccc                                                                                                        >  1:744784/1‑134 (MQ=255)
                                            ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                                   <  1:297338/151‑1 (MQ=255)
                                                aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg                                                                                               <  2:1205863/151‑1 (MQ=255)
                                                 gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                                                              <  1:341501/151‑1 (MQ=255)
                                                   gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTTAAACACGGTAGCCCCCCCAACAGAAAGCGGGGTGGGCAAACCGAGCTGc                                                                                              >  2:1057907/1‑149 (MQ=255)
                                                   gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCACATCAGAAAGCGGGGGGCGCCCACCgcgc                                                                                                 >  1:367929/1‑146 (MQ=255)
                                                    tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGGCCCCACATCAGAAAAAGGGGGGCGCACACCGAGCTGCg                                                                                             >  2:825831/1‑149 (MQ=255)
                                                        aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCCCACCAGAAAGCGGGGGGCGCACACCGAGCTGCGTAg                                                                                          >  1:1136656/1‑148 (MQ=255)
                                                        aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTTACACACGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCGAGCTGCGTAGCg                                                                                        >  2:207504/1‑150 (MQ=255)
                                                            gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGgtc                                                                                     >  2:1051053/1‑149 (MQ=255)
                                                              cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACCGGAGCCCCCCCATCAGAAAACGGGGGGCGCACCCCCCGCTGGGGAGCGAACATcc                                                                                >  2:226188/1‑151 (MQ=255)
                                                                gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCCAGCTGCGTAGCGATCATccc                                                                               >  1:1042581/1‑150 (MQ=255)
                                                                gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGGGGGAAAAATCACCACCGGGAGCCCCCCCACCAAAAAACGTGGGGGGCACACCGAGCCCCGTGGCggacaaccc                                                                                >  1:679538/1‑142 (MQ=255)
                                                                 aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACCGTAGCCCCCACATCAGAAAGCGGGGGGCGCACACCCAGCTGCGGAGCGATCATCCCtt                                                                             >  2:402817/1‑151 (MQ=255)
                                                                 aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGCAGCCCCCACATCCAAAAGCGTGGGGCGCACACCGAGCTGGGTAGCGGTCATccc                                                                               >  2:53797/1‑149 (MQ=255)
                                                                 aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAATTCACACACCGGATGCCCCACAACAGAAAACGGGGTGGGGACACCGAACTGCGTAGGGGTCACTCCCt                                                                              >  1:459666/1‑151 (MQ=255)
                                                                  tCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGGGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                                                               >  2:54992/1‑148 (MQ=255)
                                                                   cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCtt                                                                             >  2:1100703/1‑149 (MQ=255)
                                                                   cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                                                            <  2:657556/150‑1 (MQ=255)
                                                                     aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTTACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                                                         <  1:605805/151‑1 (MQ=255)
                                                                      aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACGGTAGCCCCCCCATCAGAAAGCGGGGGGGGCACACCGCGCTGCGGAGCGATCATcccct                                                                             >  1:1166955/1‑144 (MQ=255)
                                                                       aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAt                                                                       <  2:242033/151‑1 (MQ=255)
                                                                         tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTCGCGCGCGGGAAAAATCACACAAGGTAGCCCCCACCTCAGAAAAGGTGGGGCGCACCCCGGGCTGGGTATCGGTCAAccc                                                                               >  2:433895/1‑141 (MQ=255)
                                                                          cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaa                                                                     <  2:348440/150‑1 (MQ=255)
                                                                            tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCTGGCGGAAAAAATGACACAAAGTGGGCCCCCCATCAGCAAGCGAGGGGCGCCCACCGAGCTACGTAGCGATCATCCCCTGCTGATaaa                                                                    >  1:231314/1‑149 (MQ=255)
                                                                             aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                 <  2:1145039/151‑1 (MQ=255)
                                                                             aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                 <  2:222943/151‑1 (MQ=255)
                                                                                      ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGGAACCCCCCCATCAGAAAACGGGGGGCGCCCACCGAGCTGCGTAGCGATCATCCCTTGCTGGTAAAACAGCACTgg                                                          >  1:872113/1‑149 (MQ=255)
                                                                                             aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                 <  2:317451/151‑1 (MQ=255)
                                                                                              ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGCGCCGGAAAATTTACACACGGGGGGCCCCACATCAGAAAGCGTGGGGCGCCCACCCGGCTGCGTTGTGAACATCCCTTGCTGATAAACCAACACGTGcacccc                                                    >  2:41272/1‑145 (MQ=255)
                                                                                                    ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGGGTCGCGTGGAAAAAATCAAAAAACGCAGCCCCCAACACAGAACGCGGGGGGCGCACACCCAAACGGGTAGCGATAATCACTATCTGATAAACAAGCGCTGGCCGCACTTCGCgg                                             >  1:969705/1‑148 (MQ=255)
                                                                                                     gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACACCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAt                                         <  2:1183041/151‑1 (MQ=255)
                                                                                                     gaggaAAGAAAAGAGCGACCATTTTGTGGCGGGGGTGCGCGGGAAAATTACAACAACGGAGTCCGCACCTCAGAAACCGCGGGGCGCCCAACCAGCTCGGGAGCCACCCTACCCTTCTGATTAACCACCACCgggcccccgttg                                                >  1:309369/1‑134 (MQ=255)
                                                                                                      aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGGAGCCCCCCCCTCAGAAAACGTGGGGGGCCCACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACCTCGCggg                                            >  2:633173/1‑147 (MQ=255)
                                                                                                      aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGAt                                         >  2:474612/1‑150 (MQ=255)
                                                                                                        gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCCCATCAGAAAGCGTGGGGGGCCCACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGAtgg                                       >  2:82424/1‑150 (MQ=255)
                                                                                                        gaaagaaaAGAGCGACCATTTTGTGGCGGAGGCGGGCGGGAAAAATTACCAACGGTACCCCCCCCCTCCAAAAGCGTGGGGCGCACACCCGGCTGCGCATCGCTCAACCCCTGCTGATAAAACAACACAGGATCCCCCTCGCggg                                            >  1:1227787/1‑145 (MQ=255)
                                                                                                         aaagaaaagaGCGACCATTTTGTGGCGGTTGCGCGCGGGAAAAATCAAACACGGGAGCCCCCACAACAGAAAACGGGGGGCGCACACCCAGCTGCGGAGCGATCATTCCTTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGAtgg                                       >  2:432775/1‑149 (MQ=255)
                                                                                                         aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACGGTAGCCCCCCCATCAGAAAGCGGGGTGCGCACACCGAGCCGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCCCCTCGCGGGTAtgg                                       >  1:1189511/1‑149 (MQ=255)
                                                                                                           agaaaagaGCGACCATTTTGTGGAGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGc                                   <  2:680372/151‑1 (MQ=255)
                                                                                                            gaaaagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAAAGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAACCAGCACCGGCACCACTTCGCGGGGAttgggggt                                  >  2:836434/1‑144 (MQ=255)
                                                                                                             aaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCt                                  <  1:465472/150‑1 (MQ=255)
                                                                                                               aagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTACACACGGTAGCCCCCACATCAAAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACCTCGCGGGGATGGGGGCtgt                                >  2:625669/1‑150 (MQ=255)
                                                                                                                  agCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGGAGCCCCCCCCTCAGAAAGCGTGGTGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCTGCCCCACCTTGCGGGGATGGTGGGGGTGAcc                            >  1:70779/1‑151 (MQ=255)
                                                                                                                  agCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAAACCACACTGGCACCACCTCGCGGGGATGGGGGGTGTGAcc                            >  2:388253/1‑151 (MQ=255)
                                                                                                                    cGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCaa                          <  2:377605/151‑1 (MQ=255)
                                                                                                                           tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                    <  1:94268/150‑1 (MQ=255)
                                                                                                                              gtggcggtggcgCGCGGGAAAATTCACACACGGTAGCCCCCCCATCAGAAAGCGTGGGGCGCCCACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGATGGGGGCTGTGACCAATGGTCgg                   >  1:135578/1‑148 (MQ=255)
                                                                                                                                 gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt             <  1:313645/151‑1 (MQ=255)
                                                                                                                                 gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt             <  1:569985/151‑1 (MQ=255)
                                                                                                                                 gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt             <  2:222189/151‑1 (MQ=255)
                                                                                                                                 gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt             <  2:888611/151‑1 (MQ=255)
                                                                                                                                 gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt             >  2:1250395/1‑151 (MQ=255)
                                                                                                                                  cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg            <  1:1064485/151‑1 (MQ=255)
                                                                                                                                      ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt        >  2:745128/1‑151 (MQ=255)
                                                                                                                                        cgcgcgGAAAAATTCACACACAGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg       >  2:868942/1‑150 (MQ=255)
                                                                                                                                         gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt     >  2:788429/1‑151 (MQ=255)
                                                                                                                                         gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt     >  2:568336/1‑151 (MQ=255)
                                                                                                                                            cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg  >  2:664413/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCG  >  NC_000913/3508514‑3508804

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