Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,591,500 A→G 9.2% intergenic (‑221/‑157) yjiY ← / → tsr putative transporter/methyl‑accepting chemotaxis protein I, serine sensor receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,591,5000AG 9.2% 334.9 / 20.2 120intergenic (‑221/‑157)yjiY/tsrputative transporter/methyl‑accepting chemotaxis protein I, serine sensor receptor
Reads supporting (aligned to +/- strand):  ref base A (66/43);  new base G (3/8);  total (69/51)
Fisher's exact test for biased strand distribution p-value = 5.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

Only 100 of 126 total aligned reads displayed.

ATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAG  >  NC_000913/4591353‑4591649
                                                                                                                                                   |                                                                                                                                                     
aTTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGTGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTAata                                                                                                                                                     >  1:596896/1‑150 (MQ=255)
aTTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCata                                                                                                                                                     >  2:826619/1‑150 (MQ=255)
aTTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTGATTTTAGGGATGTAAGCGGGCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCata                                                                                                                                                     <  1:1141711/150‑1 (MQ=255)
       gCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACata                                                                                                                                              >  2:1778631/1‑150 (MQ=255)
         tctCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTCAGGTAAATAACTGAGGGGTTATTTTAGGGATGTAAGCGTGCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGGTTAAGTGACGATAAATATGTGATATGAATCACATAtt                                                                                                                                            <  1:1148917/148‑1 (MQ=255)
         tATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATAtt                                                                                                                                            <  1:1664005/150‑1 (MQ=255)
            ccTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTtat                                                                                                                                         <  1:1428241/150‑1 (MQ=255)
              tGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACTATAAATATGTGATTCATATCACATATTTatat                                                                                                                                       >  1:277070/1‑150 (MQ=255)
              tGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTtatcg                                                                                                                                       <  2:346018/150‑3 (MQ=255)
              tGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTtatcg                                                                                                                                       <  1:1048864/150‑3 (MQ=255)
              tGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGCTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTtatcg                                                                                                                                       <  1:752439/150‑3 (MQ=255)
                  gTGATGTTGCCAGATTTACGCGCTCGTAATGTAAGGTAAATAACTGAGTGGTTATTTTACGGATGTAATCGGTGAGTTTTGCGGATGAGCGGGAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATtgtg                                                                                                                                   <  2:942903/150‑1 (MQ=255)
                  gTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATtgtg                                                                                                                                   <  2:601860/150‑1 (MQ=255)
                  gTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATtgtg                                                                                                                                   <  2:200442/150‑1 (MQ=255)
                   tGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGa                                                                                                                                  <  2:1035061/150‑1 (MQ=255)
                   tGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTtatcgtcact                                                                                                                                  <  2:393556/150‑8 (MQ=39)
                    gATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTATAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGaa                                                                                                                                 >  1:832384/1‑150 (MQ=255)
                         tGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGAAGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAAtt                                                                                                                            >  2:739646/1‑150 (MQ=255)
                             agagTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCTGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATg                                                                                                                        >  2:79340/1‑150 (MQ=255)
                                gTTACGCGCGCGTAAAGTAAGGTATATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATGACATATTTATATTGTGAATAATGTATGTaa                                                                                                                     >  1:307010/1‑150 (MQ=255)
                                gTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTaa                                                                                                                     >  1:1024325/1‑150 (MQ=255)
                                    cgcgcgcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCTGTTGAGCGCCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATTTAaataa                                                                                                                 >  2:1226824/1‑150 (MQ=255)
                                    cgcgcgcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAaataa                                                                                                                 >  1:44378/1‑150 (MQ=255)
                                      cgcgcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATg                                                                                                               >  2:1599814/1‑150 (MQ=255)
                                          cgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTTAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATATTGGCtt                                                                                                           >  1:327282/1‑150 (MQ=255)
                                          cgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGAGAATAATTTATGTAAATAATGGCtt                                                                                                           >  1:1805076/1‑150 (MQ=255)
                                           gTAAAGTAAGGTAAATAACTGAGTGGTTATTTAAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAAATTAAGTAAATAATGGCttt                                                                                                          >  2:274447/1‑150 (MQ=255)
                                            tAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATAACATATTTATATTGTGAATAATTTATGTAAATAATGGCtttt                                                                                                         >  1:142504/1‑150 (MQ=255)
                                             aaaGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCCCATATTTATATTGTGAATAATTTATGTAAATAATGGCttttt                                                                                                        <  2:794158/150‑1 (MQ=255)
                                               aGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTaa                                                                                                      >  1:261974/1‑150 (MQ=255)
                                                ttAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTaaa                                                                                                     <  1:1372386/149‑1 (MQ=255)
                                                gTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTaaa                                                                                                     >  1:1055978/1‑150 (MQ=255)
                                                gTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGAGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTaaa                                                                                                     >  1:1188053/1‑150 (MQ=255)
                                                 tAAGGTAAATAAATGAGTGGATATTTTAGGGATGTAAGCGGTAAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAAAAATGGCTTTTTAAAt                                                                                                    <  1:710592/150‑1 (MQ=255)
                                                  aaGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAtt                                                                                                   <  2:395334/150‑1 (MQ=255)
                                                  aaGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAtt                                                                                                   >  2:1507837/1‑150 (MQ=255)
                                                   tggTAAATAACTGAGTGTTTATTTTATGGATGTAAGCGGTGAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCAGATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTc                                                                                                  <  2:1481228/149‑1 (MQ=255)
                                                   aGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCTGTCAGTTTTGCGGTTGAGCGGCAAATGCGTCGTTTAAGTTACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTc                                                                                                  <  2:1291575/150‑1 (MQ=255)
                                                   aGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAtta                                                                                                  >  1:776820/1‑149 (MQ=255)
                                                    ggTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                                                                                                 >  1:1311219/1‑150 (MQ=255)
                                                     gTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTTag                                                                                                <  2:499549/150‑1 (MQ=255)
                                                       aaaTAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCagag                                                                                              >  1:531939/1‑150 (MQ=255)
                                                         aTAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAgtg                                                                                            <  1:1746515/150‑1 (MQ=255)
                                                          tAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAgtgt                                                                                           >  2:405490/1‑150 (MQ=255)
                                                           aaCTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAgtgtt                                                                                          >  2:754556/1‑149 (MQ=255)
                                                           aaCTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGtgtg                                                                                          <  1:1517520/150‑1 (MQ=255)
                                                             cTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAATTAATGGCTTTTTAAATTCAGAGTGTGaa                                                                                        >  1:230635/1‑150 (MQ=255)
                                                             cTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGaa                                                                                        >  1:190614/1‑150 (MQ=255)
                                                             cTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCAAATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGaa                                                                                        >  1:1674289/1‑150 (MQ=255)
                                                              tGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAAt                                                                                       <  2:1268284/150‑1 (MQ=255)
                                                                   ggTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAAt                                                                                  <  1:1484201/150‑1 (MQ=255)
                                                                       aTTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCAc                                                                              <  1:846020/150‑1 (MQ=255)
                                                                         tttAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTc                                                                            >  1:795084/1‑150 (MQ=255)
                                                                          ttAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCg                                                                           >  1:772569/1‑150 (MQ=255)
                                                                           tAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCgg                                                                          <  2:1658947/150‑1 (MQ=255)
                                                                           tAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTCAAATTCAGAGTGTGAATAAAATTCACTCgg                                                                          <  1:721247/150‑1 (MQ=255)
                                                                           tAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAACGCGTCGTTTAAGTGACGATAAATATGTGATTCATAACACATATTTAAATTGTGAATAATTTATGTAAATAATTTCTTTTTAAATTCAGATACAGAATAAAAATCACTCgc                                                                          >  2:946695/1‑149 (MQ=255)
                                                                            aGGGATGTAAGCGTTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAAAAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGc                                                                         <  2:717033/150‑1 (MQ=255)
                                                                            aGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATTTTATTCAcactct                                                                         <  2:324035/150‑6 (MQ=255)
                                                                             gggATGTAAGCGGTCAGTTTTGCGGTTTAGCTGCAAAGTCGTCGTTTAAGAGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCg                                                                        <  1:339864/150‑1 (MQ=255)
                                                                                   tAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAAtct                                                                  >  2:839416/1‑150 (MQ=255)
                                                                                   tAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAAtct                                                                  >  2:627458/1‑150 (MQ=255)
                                                                                   tAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAAtct                                                                  >  2:681312/1‑150 (MQ=255)
                                                                                      gCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTAACTCGGCGTAATCTCcg                                                               >  2:1139423/1‑150 (MQ=255)
                                                                                        ggTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCcgcg                                                             >  2:1805208/1‑150 (MQ=255)
                                                                                          tCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCggg                                                           >  1:727618/1‑150 (MQ=255)
                                                                                             gTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGata                                                        >  2:1392796/1‑150 (MQ=255)
                                                                                               tttGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATAtt                                                      <  1:1769404/150‑1 (MQ=255)
                                                                                                   cGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATa                                                  <  1:1778631/150‑1 (MQ=255)
                                                                                                    ggTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATaa                                                 >  2:511151/1‑150 (MQ=255)
                                                                                                        tagaggCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGtt                                             <  1:1544945/146‑1 (MQ=255)
                                                                                                          gCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGtttt                                           >  2:967575/1‑150 (MQ=255)
                                                                                                          gCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGtttt                                           >  2:1680594/1‑150 (MQ=255)
                                                                                                          gCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGtttt                                           <  1:1543234/150‑1 (MQ=255)
                                                                                                          gCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGtttt                                           >  2:1181116/1‑150 (MQ=255)
                                                                                                            ggCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTtttc                                         >  2:6643/1‑150 (MQ=255)
                                                                                                             gCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTtttcc                                        <  1:1005618/150‑1 (MQ=255)
                                                                                                             gCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTtttcc                                        >  2:260745/1‑150 (MQ=255)
                                                                                                                   gCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCa                                  <  1:521211/150‑1 (MQ=255)
                                                                                                                    cgtcgtTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAg                                 <  1:989579/150‑1 (MQ=255)
                                                                                                                    cgtcgtTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTACTTTCCAg                                 >  2:1316791/1‑150 (MQ=255)
                                                                                                                         tttAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGa                            <  2:531939/150‑1 (MQ=255)
                                                                                                                            aaGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGaaaa                         >  1:1143037/1‑150 (MQ=255)
                                                                                                                             aGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAAt                        <  2:399435/150‑1 (MQ=255)
                                                                                                                               tGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCt                      >  1:423526/1‑150 (MQ=255)
                                                                                                                               tGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCt                      >  1:530433/1‑150 (MQ=255)
                                                                                                                                gACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTACAGGCCGAAAATCtt                     >  1:957318/1‑150 (MQ=255)
                                                                                                                                caataTAAATATGTGATATGAATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCtt                     >  1:1283295/5‑150 (MQ=255)
                                                                                                                                  cGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGc                   >  1:780122/1‑150 (MQ=255)
                                                                                                                                  cGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGGGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGc                   <  2:1171536/150‑1 (MQ=255)
                                                                                                                                   gATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCa                  >  1:971124/1‑150 (MQ=255)
                                                                                                                                    aTAAATATGTGATATGAATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCAt                 >  2:472237/1‑150 (MQ=255)
                                                                                                                                      aaaTATGTGATATGAATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCg               >  1:911028/1‑150 (MQ=255)
                                                                                                                                       aaTATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCTCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCgg              >  1:998947/1‑150 (MQ=255)
                                                                                                                                        atatGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGt             >  2:1500824/1‑150 (MQ=255)
                                                                                                                                         tatGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTc            <  1:1260355/150‑1 (MQ=255)
                                                                                                                                          atggtaatCATATAACATATTTATATTGTGTATAATTTATTTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAATCTTGCATCGGTCca          <  1:1501134/143‑1 (MQ=255)
                                                                                                                                          atGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTcc           >  2:692605/1‑150 (MQ=255)
                                                                                                                                             tgATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCcaca        >  1:1163926/1‑150 (MQ=255)
                                                                                                                                                   atatCACATATTTATATTGTGAATAATTTATGTAAATAATTGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCAGCGGGATATTCATAAAGTTATTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACATGAAag  >  2:749777/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAG  >  NC_000913/4591353‑4591649

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