Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,481,771 G→T 12.3% P136T (CCG→ACG)  emrY ← putative multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,481,7710GT12.3% 898.8 / 19.4 283P136T (CCG→ACG) emrYputative multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (139/99);  new base T (2/43);  total (141/142)
Fisher's exact test for biased strand distribution p-value = 1.88e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

GACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAA  >  NC_000913/2481622‑2481920
                                                                                                                                                      |                                                                                                                                                     
gACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGTAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCgg                                                                                                                                                       >  2:4717406/1‑150 (MQ=255)
gACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCgg                                                                                                                                                       >  2:5408643/1‑150 (MQ=255)
gACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCgg                                                                                                                                                       >  1:4727336/1‑150 (MQ=255)
 aCGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGa                                                                                                                                                      <  1:3378065/150‑1 (MQ=255)
 aCGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGa                                                                                                                                                      <  2:3058872/150‑1 (MQ=255)
 aCGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGa                                                                                                                                                      >  1:4501401/1‑150 (MQ=255)
   gATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATa                                                                                                                                                    <  1:5244071/150‑1 (MQ=255)
        tCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTcc                                                                                                                                               >  1:4549498/1‑150 (MQ=255)
         ccccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCt                                                                                                                                              <  2:698458/150‑1 (MQ=255)
          cccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTTCTGGCGGATAATTCCt                                                                                                                                              >  2:1707144/1‑150 (MQ=255)
          cccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCtt                                                                                                                                             <  1:65891/150‑1 (MQ=255)
           ccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGACAGAGCAAATGTTCTTTTTTCTGGCGTATAATTCCTta                                                                                                                                            >  2:824465/1‑150 (MQ=255)
            cATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaa                                                                                                                                           <  1:897153/150‑1 (MQ=255)
            cATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaa                                                                                                                                           <  2:4679786/150‑1 (MQ=255)
             aTAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaat                                                                                                                                          <  1:3791530/150‑1 (MQ=255)
              tAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaata                                                                                                                                         <  1:1018854/150‑1 (MQ=255)
                  gACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaataataa                                                                                                                                     <  2:1381984/150‑1 (MQ=255)
                      ttGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAActc                                                                                                                                 >  1:2566781/1‑150 (MQ=255)
                       tGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAActct                                                                                                                                <  2:5143953/150‑1 (MQ=255)
                         aTTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCtgt                                                                                                                              >  2:4492971/1‑150 (MQ=255)
                         aTTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCtgt                                                                                                                              <  2:2612499/150‑1 (MQ=255)
                           tAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGa                                                                                                                            >  2:4095163/1‑150 (MQ=255)
                           tAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGa                                                                                                                            <  2:1161682/150‑1 (MQ=255)
                              aaaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAg                                                                                                                         >  2:4917045/1‑150 (MQ=255)
                              aaaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAg                                                                                                                         >  2:191407/1‑150 (MQ=255)
                               aaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGTGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGt                                                                                                                        >  1:3882100/1‑150 (MQ=255)
                               aaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGt                                                                                                                        <  1:4402904/150‑1 (MQ=255)
                               aaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGt                                                                                                                        <  1:1456538/150‑1 (MQ=255)
                                  atCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGa                                                                                                                     >  1:2409978/1‑150 (MQ=255)
                                  atCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGACCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGa                                                                                                                     <  2:2831301/150‑1 (MQ=255)
                                   tCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGaa                                                                                                                    >  2:12245/1‑150 (MQ=255)
                                    ccAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAAt                                                                                                                   >  2:2866160/1‑150 (MQ=255)
                                        ccccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCATATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAAc                                                                                                               <  1:3705791/150‑1 (MQ=255)
                                        ccccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAAc                                                                                                               <  1:838001/150‑1 (MQ=255)
                                            aGCTAAAGTTATCACAAATATAACCACCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGcc                                                                                                           >  2:2924028/1‑150 (MQ=255)
                                             gCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGccc                                                                                                          >  1:306232/1‑150 (MQ=255)
                                              cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGcccc                                                                                                         >  2:4357563/1‑150 (MQ=255)
                                              cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGcccc                                                                                                         >  1:1764499/1‑150 (MQ=255)
                                                    ttATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCAtt                                                                                                   >  2:5377529/1‑150 (MQ=255)
                                                     tATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTa                                                                                                  >  2:4843516/1‑150 (MQ=255)
                                                      aTCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACAATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAACCTCTGTGACAGTGGAATTAACGGCCCCGCCATTaa                                                                                                 <  1:5470612/150‑1 (MQ=255)
                                                         acaAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAccc                                                                                              >  2:4219559/1‑150 (MQ=255)
                                                         acaAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAccc                                                                                              >  1:4868851/1‑150 (MQ=255)
                                                         acaAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAccc                                                                                              <  1:3968648/150‑1 (MQ=255)
                                                         acaAATAAACCCGCCCAAAAACGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCGGGCGGATAATTCCTTAATAATAAACACTGTGACAGTGGAATTAACGGCCCCGCCATTAAccc                                                                                              <  1:1256310/150‑1 (MQ=255)
                                                           aaaTATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCt                                                                                            >  2:857937/1‑150 (MQ=255)
                                                           aaaTATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCt                                                                                            <  2:4312102/150‑1 (MQ=255)
                                                              tataACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGaa                                                                                         <  1:3463081/150‑1 (MQ=255)
                                                               ataAGCGCCCAATATCGGCCCACATTTCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGGCCCGCCATTAACCCCTGAAc                                                                                        <  2:5132141/150‑1 (MQ=255)
                                                                    cGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGActc                                                                                   >  2:55147/1‑150 (MQ=255)
                                                                    cGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGActc                                                                                   >  1:3219761/1‑150 (MQ=255)
                                                                        cAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                               >  1:5437284/1‑150 (MQ=255)
                                                                        cAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                               >  1:5196923/1‑150 (MQ=255)
                                                                        cAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                               >  1:3659156/1‑150 (MQ=255)
                                                                        cAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                               >  2:5443805/1‑150 (MQ=255)
                                                                                gCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCa                                                                       >  2:254270/1‑150 (MQ=255)
                                                                                 cccAAATATCGGAGCTATAATCACGTTAATTTACAATATTGCAAGAGAAATTTTTATTTTTTCTGGGGTATATTTCCTTAATAATAAACTCTGTGACTGTGGATTTAACGGCCCCGCCATTAACCCCTTAACGACTCTAAAAAATATCAg                                                                      <  2:3858272/150‑1 (MQ=255)
                                                                                    acaTATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcac                                                                   >  2:4525301/1‑150 (MQ=255)
                                                                                    acaTATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcac                                                                   >  1:2079177/1‑150 (MQ=255)
                                                                                         tCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCaa                                                              >  1:4164777/1‑150 (MQ=255)
                                                                                           ggAGCTATAATCACGTAATTTTACCATAATGCATGATAAATTTTTCTTTTTGCTGGCTTATAATTCATTAATAATAAACTCTTTGACTGTGTAATTAACGGCCCCGCAATTAACCCCACAACGACTCTAAAAAATATCAGCACATCAAGa                                                            <  1:762413/150‑1 (MQ=255)
                                                                                             agataTAATCAGGGTCTTTTACAATACTGCCAGTTAAATTTTTTTTTTTTGTGGGTTATATTTCTTTAATAATAAACTTTGTGAGTTTTTATTTAAGGGCCCCGCCATTACCCCCTTCACGACTCTAAAAAATATCAGCACATCAAGAtt                                                          <  1:1631033/146‑1 (MQ=255)
                                                                                             aGCGATAATCACGGTCATTGACAATAATGCCATAGAAAATTTTCTTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTGTTACTGTGGATTTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGAtt                                                          <  1:3117810/150‑1 (MQ=255)
                                                                                               cGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                        >  2:2345871/1‑150 (MQ=255)
                                                                                               cGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                        >  2:2801342/1‑150 (MQ=255)
                                                                                                gATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGTATATTTCCTTAATAATACACTCTGTGACAGTGGAATTAACGGCCCCGCATTTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGt                                                       <  2:3422072/150‑1 (MQ=255)
                                                                                                gATAATCACGGTAATTGACAATAATGCCAGAGAAAATGTTCTTTTTTCTGGCGTATATTTCCTTAATAATAAACTCTGTGACAGTGTATTTAACGGCCCGGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGt                                                       <  2:2056217/150‑1 (MQ=255)
                                                                                                  tAACCACTGTCATTGACAATAAGGCCAGAGAAAATGTTATTTTTTGTGGCGGAAAATTCATTAATAATAAACTCTGTGACAGTGTAATTAACGGCCCCGCAATTAACCCCTCAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCg                                                     <  1:2893836/150‑1 (MQ=255)
                                                                                                     tCACGGTCTTTTAAAATAATGCCATATAATTTTTTTTTTTTTCGGGCGTATAATTCATTAATATTAAACTGTGTGACTGTGTATTTAAGGGCCCCGCATTTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATa                                                  <  2:4258856/150‑1 (MQ=255)
                                                                                                     tCACGGTCTTTGACCATAATGCCAGAGCAAATGTTCTTTTTTCGGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAAGGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATa                                                  <  2:2177661/150‑1 (MQ=255)
                                                                                                      cACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAg                                                 >  1:3619693/1‑150 (MQ=255)
                                                                                                        cGGTCATTGACAATAATGCCAGAGAAAATGTTCTTTTTTGGGGCGTAAAATCCATTAATAATAAACTCTGTGACTGTGGAATTAACGGCCCCGCCATTAACCCCTGCACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGc                                               <  1:1163553/150‑1 (MQ=255)
                                                                                                        cGGACATTGACCATAATGCCAGAGAAATTTTTCTTTTTTCTGGCGTATAATTCCTTAATAATAAACTATGTGACAGTGTAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGc                                               <  2:3947446/150‑1 (MQ=255)
                                                                                                         ggTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCt                                              >  2:524509/1‑150 (MQ=255)
                                                                                                          gTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTa                                             >  1:2801376/1‑150 (MQ=255)
                                                                                                           tatttGACCATAATGCCATATCAATTGTTCTTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTTTGACTGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTac                                            <  1:2683669/147‑1 (MQ=255)
                                                                                                                gACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                       >  2:4047568/1‑150 (MQ=255)
                                                                                                                  ccATAATCCAAGACAAAATTTTCTTTTTTCTGGCGTATAATACATAAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCAATTAACCCCTCAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTaa                                     <  1:407935/150‑1 (MQ=255)
                                                                                                                   cATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCGACACATTAAt                                    >  2:498277/1‑150 (MQ=255)
                                                                                                                      acggCCCGAGAAATTGTTCTTTTTTGGGGCGGATAATTCCTTAATAATAAACTCTGTGACTGGGGAATTAACGGCCCCGCCTTTAACCCCGGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGaa                                 <  1:2040521/147‑1 (MQ=255)
                                                                                                                      aaggCCACAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGAAAGTGGAATTAACGGCCCCGCCATTAACCCCTTAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGaa                                 <  2:3135522/147‑1 (MQ=255)
                                                                                                                      aaTGCCTGAGAAATTGTTTTTTTTTCTGGCGTATAAGTCCTTAATATTAAACGAGGGGACAGTGTAATTAACGGCCCCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACCTCAAGATTGGTCGATAGGCTACACATTAATGaa                                 <  2:1517072/150‑1 (MQ=255)
                                                                                                                         gCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGAc                              >  1:2172072/1‑150 (MQ=255)
                                                                                                                          cctgagCAATTGTTCTTTTTTCTGGGGGATATTTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACa                             <  2:361142/147‑1 (MQ=255)
                                                                                                                           cAGAGCAATTGTTCTTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAg                            <  2:3219850/150‑1 (MQ=255)
                                                                                                                                     gTTATTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaaaa                  <  2:1717369/150‑1 (MQ=255)
                                                                                                                                      ttCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAg                 >  2:2847397/1‑150 (MQ=255)
                                                                                                                                      ttCTTTTTTCTAGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaaaat                 >  2:1811406/1‑149 (MQ=255)
                                                                                                                                       tCTTTTTTCTGGCTGATATTTCCTTAATAATAAACTTTGTGAATGTGTATTTAACGGCCCCGCATTTAACCCGTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGt                <  2:3693515/150‑1 (MQ=255)
                                                                                                                                        cTTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTTTTACTGTGTAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCAAAAAAAAGTGa              <  2:4874602/150‑1 (MQ=255)
                                                                                                                                          ttttttCTGGGGGATAATTCCTTAATAATAAACTCTGTTAATGTTGAATTAACGGCCCCGACATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGa              <  2:2410097/150‑1 (MQ=255)
                                                                                                                                          ttttctCTGGCGGATAATTCATTAAAAATAAAATCTTTGACAGTGTAATTAACGGCCCGGCAATTACCCCATGCACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGa              <  2:41341/150‑1 (MQ=255)
                                                                                                                                           tttttCTGGGGTATAATTCCTTAATAATAAACTCTGTGACTGTGTAATTAACGGCCCCGCATTTAACCCCTGAACGACTCTAAAAAATATCAGCAAATCAAGTTTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGAc             <  1:2956065/150‑1 (MQ=255)
                                                                                                                                           tttttCTGGGGTATAATTCATTAATAATAAACTCTGTGACAGGGGAATTAAGGGCCCCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGAc             <  1:4997144/150‑1 (MQ=255)
                                                                                                                                           tttttCTGGCGTATAATTCCTTAATAATAAACTCTGTGACAGTGTAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGAc             <  1:191407/150‑1 (MQ=255)
                                                                                                                                               ccTGGCGTATAATACCTAAATAAAAAACTCTGTGACTGTGTAATTAACGGCCCCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGaa         <  1:1931637/149‑1 (MQ=255)
                                                                                                                                                cgggCGGATAAGTCATTAAAAATAAACTCTGTGACAGTGTAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAAAATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGaaa        <  1:5313911/148‑1 (MQ=255)
                                                                                                                                                 tGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAg       >  2:2275797/1‑150 (MQ=255)
                                                                                                                                                   gCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTa     >  1:2356165/1‑150 (MQ=255)
                                                                                                                                                   gCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTa     >  1:1395145/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAA  >  NC_000913/2481622‑2481920

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: