Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,375,1230TG44.0% 149.3 / 17.4 85G85G (GGT→GGGdcuDputative transporter
Reads supporting (aligned to +/- strand):  ref base T (6/41);  new base G (37/0);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 6.66e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAGCGCCTC  >  NC_000913/3374976‑3375272
                                                                                                                                                   |                                                                                                                                                          
gatgTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTa                                                                                                                                                          <  1:753599/150‑1 (MQ=255)
   gtgCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTAt                                                                                                                                                         >  1:15837/1‑148 (MQ=255)
      cGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTGTGGCCg                                                                                                                                                    >  2:320421/1‑150 (MQ=255)
            cTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCg                                                                                                                                                    >  2:766977/1‑144 (MQ=255)
              ggTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACAt                                                                                                                                              <  2:272797/148‑1 (MQ=255)
              ggTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGGGGGTTTGGCCggttcaaggg                                                                                                                                           >  2:97252/1‑142 (MQ=255)
                     gTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGGGGGTATGCCCGCCACATTGGGCGCAt                                                                                                                                   >  2:44315/1‑150 (MQ=255)
                        cGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAggg                                                                                                                               <  1:536684/151‑1 (MQ=255)
                        cGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAggg                                                                                                                               <  1:131516/151‑1 (MQ=255)
                         ggtggtGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCTACATGGAGCGCATAgggg                                                                                                                              >  1:822193/1‑151 (MQ=255)
                             gtgtACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAg                                                                                                                          <  2:831305/151‑1 (MQ=255)
                               gtACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATTTCCCCCCCATGGAGGGGCTAGGGGGCAgg                                                                                                                         >  1:870931/1‑149 (MQ=255)
                                  cTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATTCCCGCTACATGGGGCGCGAGGGGGGCAGgcgc                                                                                                                      >  2:517918/1‑150 (MQ=255)
                                    cGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGCCCCCCCCCTGGGGCGCGAAAGGGGccgggggg                                                                                                                      >  2:536358/1‑143 (MQ=255)
                                      aTCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTTGCCGCCACATGGGGGGCCAGAGGGGccg                                                                                                                           >  2:537867/1‑141 (MQ=255)
                                       tCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGGCCCCCACATGGAGCGCCTAGGGGGCAggtccg                                                                                                                     >  2:542049/1‑141 (MQ=255)
                                        cccACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATgg                                                                                                              <  1:176179/151‑1 (MQ=255)
                                          caccaAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCCCCCCCTGGAGCGCATAGGGGGCCGGcgcgcggggg                                                                                                               >  1:906437/1‑145 (MQ=255)
                                            ccaAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGGGGGTTTGCCCGCTACATGGAGCGGATAGGGGGCAGGCGCGCGATgggg                                                                                                            >  1:849006/1‑147 (MQ=255)
                                             caAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGCGGGTATGTCCGCTACATGGAGAGCCTAGGGGCCAGACGCGCGGTGGCGAg                                                                                                          >  2:138613/1‑150 (MQ=255)
                                                aaGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTg                                                                                                      <  2:335454/151‑1 (MQ=255)
                                                  gcagcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTCCCCCCACATGGGGCGCCTAGGGGGGAGTCGCGCGGTGGGGGGgttgtt                                                                                                   >  1:63437/1‑151 (MQ=255)
                                                     gcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAg                                                                                                <  2:319924/151‑1 (MQ=255)
                                                     gcGGCTACTTGAGCGTCTATATTTATAATGAAATCCTGCGGATGCTGTCCAACCGCATTGCCGGATTGGGGCTCTCGATTAGGGCGGTGGGCGGTTAGGCACGCTATATGGAGCGCATAGGGGCTAGTCGCGCGATTGTGAGCTTGTTAAg                                                                                                <  2:111704/151‑1 (MQ=255)
                                                           aCTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAgccgc                                                                                           <  1:832355/149‑1 (MQ=255)
                                                           aCTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGccc                                                                                         <  1:175138/151‑1 (MQ=255)
                                                           aCTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGccc                                                                                         <  2:25667/151‑1 (MQ=255)
                                                            cTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTTCCGCTACCTGGGGCGCATAGGGGGCAGGCGCGCGGTGGGGGGGTTGTTAAAACGccc                                                                                         >  1:558622/1‑150 (MQ=255)
                                                              tGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGt                                                                                       <  1:615743/150‑1 (MQ=255)
                                                                  cgtcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCCCCCTGGAGGGCATAGGGGGCAGGCGCGCGATGGGGAGGTTGTTTAACCGCCCGTTtaa                                                                                   >  1:503619/1‑150 (MQ=255)
                                                                    tcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTTCCCGCCCCATGGAGGGCCTAGGGGGCCGGCGCGGGATGGGGGGGGTGTTTAGCCGcccccttaaa                                                                                  >  1:448965/1‑142 (MQ=255)
                                                                      gATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGCGGGCGGGTATTCGCGCCACAAGGGGAGCAAAGGGGGCAGGCGCGCGCTGGGGGAGGTGTTTAACCCCCCGGTTaaaaac                                                                                >  1:486120/1‑148 (MQ=255)
                                                                        tatTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCCCCACATGGGGCGCCTAGGGGGCAGTCGCGCGGTGGGGGGGTTGTTTAGCCGCCCGGTtaaa                                                                                  >  1:33925/1‑145 (MQ=255)
                                                                         atTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCTGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGGCCGCTACATGGAGGGCATAGGGGGCCGGCGCGCGACGGGGAGCGTGTTTAAGCGCCCCCTAAAAATCAtt                                                                            >  1:497084/1‑150 (MQ=255)
                                                                         atTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGGCCGCTACAAGGAGGGCCAAGGGGGCAGGCGCGGGGTGGGGAGGGTGTTTAGGCCCCCCGTTAAAATCATtc                                                                           >  2:201188/1‑151 (MQ=255)
                                                                           ttATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcg                                                                          <  1:61460/150‑1 (MQ=255)
                                                                             ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcgctc                                                                       <  1:163755/151‑1 (MQ=255)
                                                                             ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcgctc                                                                       <  1:369982/151‑1 (MQ=255)
                                                                             ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcgct                                                                        <  2:363209/150‑1 (MQ=255)
                                                                               aaTGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcgctcgc                                                                     <  1:96592/151‑1 (MQ=255)
                                                                                   aaaTCCTGCGCATGCTGTCCACCCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGta                                                                 <  1:579262/151‑1 (MQ=255)
                                                                                    aaTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGtat                                                                <  2:665274/151‑1 (MQ=255)
                                                                                     aTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGAGGGGGGGTATGTCCGCTCCCTGGGGCGCGTAGGGGGCCTTCGCGCGGTGGTGGGGGGGTTAGGACGCCCGTTCAAAAACATttgctc                                                                       >  2:224110/1‑143 (MQ=255)
                                                                                      tCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGtata                                                               <  2:394065/150‑1 (MQ=255)
                                                                                       ccTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTTCCCGCCCCATGGAGCGCATAGGGGGCAGGCGCGCGATGGGGGGCTTGTTTAGGCGCCCCCTAAAAAACATTTGCTCCCCGGAtatt                                                             >  1:451420/1‑151 (MQ=255)
                                                                                         tGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGTCCCCCACATGGAGCGCCTAGGGGGCAGTCGCGCGATGGGGAGCTTGTTTAGACGCCCGTTAAAAATCATTCGCTCGCCGTAtattat                                                           >  1:234493/1‑151 (MQ=255)
                                                                                          gcgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTAtattat                                                           <  2:259392/150‑1 (MQ=255)
                                                                                          gcgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTAtatta                                                            <  1:744684/149‑1 (MQ=255)
                                                                                          gcgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATattatt                                                          <  2:158325/151‑1 (MQ=255)
                                                                                           cgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGGGGGTATGCCCGCCACATTGGGGGCCTAGGGGGCAGGCGCGCGGTGGGGAGCTTGTTTAGCCGCCCCGTTAAAATCATTTGCTCGCCGTATattatt                                                          >  1:122597/1‑150 (MQ=255)
                                                                                            gcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCTACATGGGGCGCCTAGGGGGCCGTCGCGCGGCGGGGGGGTTGTTAAAACGCCCCCTTAAAAACCTTCGCTCGCCCGATattatttt                                                        >  1:567654/1‑149 (MQ=255)
                                                                                             cATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTCCCGCCACATGGAGGGCCTAGGGGGCCGGCGCGCGGTGGGGAGGTTGTTTAAGCGCCCGGTAAAAAACATTTGCTCGCCGGATattttttt                                                        >  1:74752/1‑148 (MQ=255)
                                                                                             cATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTg                                                       <  2:552668/151‑1 (MQ=255)
                                                                                              aTGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGt                                                      <  2:652120/151‑1 (MQ=255)
                                                                                                  tGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGTTTTGCCCCCCACATGGAGCGCCTAGGGGGCAGTCGCGCGATGGGGGGGTTTTTTAGCCGCCCCCTTAAACTCATTTGCTCGCCGTGTattttttt                                                        >  1:351660/1‑143 (MQ=255)
                                                                                                   gTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCgg                                                   <  1:881343/149‑1 (MQ=255)
                                                                                                      cAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACt                                              <  1:535679/151‑1 (MQ=255)
                                                                                                       aaCCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGCCTGCCCCCCCCCCGGCGCCCCTAGGGGCCCGCCCCGCCCGGGGGGGCCTGTTTAGCCCCCCCCCAAAACCCATTCCCCCCCCCCCTATTATTCTTTccg                                                   >  2:360075/1‑143 (MQ=255)
                                                                                                         cctcATTGCCGGATTGGGACTGTCGATTATGGCGGTAGGTGGATATGCCCGCTACATGGGGCGCATAGGGGCCACTCGCGCGATGGTGAGTTTGTTATGCCGCCAGTTAAAACTCATTCGATCGCCGTATATTATTCTGTCGGCAACTTAc                                           <  2:534360/148‑1 (MQ=255)
                                                                                                         ccGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCCCTACATGGAGCGCATAGGGGGCCGTCGCGCGGGGGGGAGCGTTTTTAGGCGCCCGTTAAAAATCATTCGCTCGCCGCATATTTTTTTTTCGGGAAcct                                             >  1:308267/1‑147 (MQ=255)
                                                                                                           gCATTGCCGGATTGGGGCTGTCGATTATGGCGGTCGGCGGGTATGACCGCCAGATGGAGCGCGTAGGGGGGAGGCGCGCGGGGGTGGGGTTGTTATGACGCCCCGTAAAAATCCTTCGCTCGCCGTCTATTTTTCTGTCTGcgaaatccg                                          >  2:497838/1‑142 (MQ=255)
                                                                                                                gCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGCCCGCCCCATGGGGGGCATAGGGGGCCGGCCCGCGATGGTGGGCTTTTTTAACCCCCCCGTAAAAATCATTTGCTCGCCGCATATTTTTTTTTCGGCCAACTACGTCa                                       >  1:536636/1‑148 (MQ=255)
                                                                                                                gCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGGCCCCCACCTGGGGGGCGCAGGGGGCCGGCGCGCGAGGGGGGGGTTGTTTAGTCGCCCCTTAAAAATCATTTGCTCGCCGGATATTTTTTTTTCTGCCCCTTACGGCAt                                      >  2:876373/1‑149 (MQ=255)
                                                                                                                gCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCg                                    <  2:586452/151‑1 (MQ=255)
                                                                                                                gCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGGCCGCTACATGGAGCGCATAGGGGGCCGTCGCGCGGAGGGGAGATTTTTTAACCGCCCGCTAAAAATCAATTGCTCGCCGGATATTTTTTTTTTCGGGACTTACGTCATCg                                    >  1:45100/1‑151 (MQ=255)
                                                                                                                 ccGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTTTTCCCCCCACATGGAGCGCATAGGGGGCAGGGGCCCCGTGGGGGGGTTGTTAAGCCGCCCCTTAAAAAACCTTCGCTCGCCCGATTTTTTTTTGTCGGCAACCTACGTCATCgg                                   >  1:373859/1‑151 (MQ=255)
                                                                                                                  cGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTCCCCCTACAAGGAGGGCCAAGGGGGCAGGCGCGCGATGGGGGGGTTGTTTAGCCGCCCGGTAAAAATCATTCTCTCGCCGGATATTTTTTTTTCGGCCACCTATGTCATCg                                    >  1:287459/1‑149 (MQ=255)
                                                                                                                       tGGGGTTGCCAATTATGGCGGTGGGCGGTTATGCCCGCTGCATGCAGCGCATAGGGGCCGCTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAGCCCCATTCGCGCGCCGTATATTATTCTGTCGTCGACTTACGTCATCGGCCaa                               <  2:13215/149‑1 (MQ=255)
                                                                                                                        ggggCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCGACATGGAGCGCATAGGGGCCAGTCGCGCGATGTTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAAt                             <  2:625233/150‑1 (MQ=255)
                                                                                                                        ggggCTGTCGATTATGGCGGGGGGGGGGGATGTCCGCCACATGGAGCGCGTAGGGGGCAGGCGCGCGATGGGGGGGTTGTTTAGCCGCCCCGTAAAAAACATTTGCTTCCCCTGTATTATTCTTTGGGCCACTTActgcctctgcccaa                              >  2:82068/1‑136 (MQ=255)
                                                                                                                          ggCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCAt                          <  1:11891/151‑1 (MQ=255)
                                                                                                                          ggCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCAt                          <  1:692290/151‑1 (MQ=255)
                                                                                                                          ggCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCAt                          <  1:792296/151‑1 (MQ=255)
                                                                                                                           gCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATg                         <  1:58703/151‑1 (MQ=255)
                                                                                                                               tCGATTAGGGCGGTGGGCGGTTATACCCGGTACATGGAGCCCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAGTCATGgcg                      <  2:740810/150‑1 (MQ=255)
                                                                                                                                     aTGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGttt                <  2:889446/150‑1 (MQ=255)
                                                                                                                                      tGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTttat              >  2:596109/1‑151 (MQ=255)
                                                                                                                                      tGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTttat              >  1:108914/1‑151 (MQ=255)
                                                                                                                                      tGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTttat              >  1:313178/1‑151 (MQ=255)
                                                                                                                                       ggCGGTGGGCGGTTATGCCCTCTACATGGGGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTTGCCGTATATTATTCTGTCGGCATCTTACGTCATCGGCCAAATCATGGCGCAGTttat              <  2:778055/150‑1 (MQ=255)
                                                                                                                                       ggCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTttatt             >  1:160571/1‑151 (MQ=255)
                                                                                                                                            tGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTAcc          >  1:327542/1‑149 (MQ=255)
                                                                                                                                              ggCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAgcg      <  1:694464/151‑1 (MQ=255)
                                                                                                                                              ggCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAgcg      <  1:479199/151‑1 (MQ=255)
                                                                                                                                                  gTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAGCGCCTc  <  2:767111/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                          
GATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAGCGCCTC  >  NC_000913/3374976‑3375272

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: