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.

GACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAA  >  NC_000913/2481622‑2481920
                                                                                                                                                      |                                                                                                                                                     
gACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGTAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCgg                                                                                                                                                       >  2:4717406/1‑150 (MQ=255)
gACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCgg                                                                                                                                                       >  2:5408643/1‑150 (MQ=255)
gACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCgg                                                                                                                                                       >  1:4727336/1‑150 (MQ=255)
 aCGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGa                                                                                                                                                      <  1:3378065/150‑1 (MQ=255)
 aCGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGa                                                                                                                                                      <  2:3058872/150‑1 (MQ=255)
 aCGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGa                                                                                                                                                      >  1:4501401/1‑150 (MQ=255)
   gATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATa                                                                                                                                                    <  1:5244071/150‑1 (MQ=255)
        tCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTcc                                                                                                                                               >  1:4549498/1‑150 (MQ=255)
         ccccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCt                                                                                                                                              <  2:698458/150‑1 (MQ=255)
          cccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTTCTGGCGGATAATTCCt                                                                                                                                              >  2:1707144/1‑150 (MQ=255)
          cccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCtt                                                                                                                                             <  1:65891/150‑1 (MQ=255)
           ccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGACAGAGCAAATGTTC‑TTTTTTCTGGCGTATAATTCCTta                                                                                                                                            >  2:824465/1‑150 (MQ=255)
            cATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTtaa                                                                                                                                           <  1:897153/150‑1 (MQ=255)
            cATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTtaa                                                                                                                                           <  2:4679786/150‑1 (MQ=255)
             aTAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTtaat                                                                                                                                          <  1:3791530/150‑1 (MQ=255)
              tAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTtaata                                                                                                                                         <  1:1018854/150‑1 (MQ=255)
                  gACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTtaataataa                                                                                                                                     <  2:1381984/150‑1 (MQ=255)
                      ttGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAActc                                                                                                                                 >  1:2566781/1‑150 (MQ=255)
                       tGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAActct                                                                                                                                <  2:5143953/150‑1 (MQ=255)
                         aTTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCtgt                                                                                                                              >  2:4492971/1‑150 (MQ=255)
                         aTTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCtgt                                                                                                                              <  2:2612499/150‑1 (MQ=255)
                           tAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGa                                                                                                                            >  2:4095163/1‑150 (MQ=255)
                           tAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGa                                                                                                                            <  2:1161682/150‑1 (MQ=255)
                              aaaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAg                                                                                                                         >  2:4917045/1‑150 (MQ=255)
                              aaaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAg                                                                                                                         >  2:191407/1‑150 (MQ=255)
                               aaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGTGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGt                                                                                                                        >  1:3882100/1‑150 (MQ=255)
                               aaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGt                                                                                                                        <  1:4402904/150‑1 (MQ=255)
                               aaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGt                                                                                                                        <  1:1456538/150‑1 (MQ=255)
                                  atCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGa                                                                                                                     >  1:2409978/1‑150 (MQ=255)
                                  atCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGACCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGa                                                                                                                     <  2:2831301/150‑1 (MQ=255)
                                   tCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGaa                                                                                                                    >  2:12245/1‑150 (MQ=255)
                                    ccAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAAt                                                                                                                   >  2:2866160/1‑150 (MQ=255)
                                        ccccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCATATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAAc                                                                                                               <  1:3705791/150‑1 (MQ=255)
                                        ccccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAAc                                                                                                               <  1:838001/150‑1 (MQ=255)
                                            aGCTAAAGTTATCACAAATATAACCACCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGcc                                                                                                           >  2:2924028/1‑150 (MQ=255)
                                             gCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGccc                                                                                                          >  1:306232/1‑150 (MQ=255)
                                              cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGcccc                                                                                                         >  2:4357563/1‑150 (MQ=255)
                                              cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGcccc                                                                                                         >  1:1764499/1‑150 (MQ=255)
                                                    ttATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCAtt                                                                                                   >  2:5377529/1‑150 (MQ=255)
                                                     tATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTa                                                                                                  >  2:4843516/1‑150 (MQ=255)
                                                      aTCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACAATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAACCTCTGTGACAGTGGAATTAACGGCCCCGCCATTaa                                                                                                 <  1:5470612/150‑1 (MQ=255)
                                                         acaAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAccc                                                                                              >  2:4219559/1‑150 (MQ=255)
                                                         acaAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAccc                                                                                              >  1:4868851/1‑150 (MQ=255)
                                                         acaAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAccc                                                                                              <  1:3968648/150‑1 (MQ=255)
                                                         acaAATAAACCCGCCCAAAAACGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCGGGCGGATAATTCCTTAATAATAAACACTGTGACAGTGGAATTAACGGCCCCGCCATTAAccc                                                                                              <  1:1256310/150‑1 (MQ=255)
                                                           aaaTATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCt                                                                                            >  2:857937/1‑150 (MQ=255)
                                                           aaaTATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCt                                                                                            <  2:4312102/150‑1 (MQ=255)
                                                              tataACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGaa                                                                                         <  1:3463081/150‑1 (MQ=255)
                                                               ataAGCGCCCAATATCGGCCCACATTTCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGGCCCGCCATTAACCCCTGAAc                                                                                        <  2:5132141/150‑1 (MQ=255)
                                                                    cGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGActc                                                                                   >  2:55147/1‑150 (MQ=255)
                                                                    cGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGActc                                                                                   >  1:3219761/1‑150 (MQ=255)
                                                                        cAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                               >  1:5437284/1‑150 (MQ=255)
                                                                        cAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                               >  1:5196923/1‑150 (MQ=255)
                                                                        cAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                               >  1:3659156/1‑150 (MQ=255)
                                                                        cAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                               >  2:5443805/1‑150 (MQ=255)
                                                                                gCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCa                                                                       >  2:254270/1‑150 (MQ=255)
                                                                                 cccAAATATCGGAGCTATAATCACGTTAATTTACAATATTGCAAGAGAAATTTTTA‑TTTTTTCTGGGGTATATTTCCTTAATAATAAACTCTGTGACTGTGGATTTAACGGCCCCGCCATTAACCCCTTAACGACTCTAAAAAATATCAg                                                                      <  2:3858272/150‑1 (MQ=255)
                                                                                    acaTATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcac                                                                   >  2:4525301/1‑150 (MQ=255)
                                                                                    acaTATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcac                                                                   >  1:2079177/1‑150 (MQ=255)
                                                                                         tCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCaa                                                              >  1:4164777/1‑150 (MQ=255)
                                                                                           ggAGCTATAATCACGTAATTTTACCATAATGCATGATAAATTTTTC‑TTTTTGCTGGCTTATAATTCATTAATAATAAACTCTTTGACTGTGTAATTAACGGCCCCGCAATTAACCCCACAACGACTCTAAAAAATATCAGCACATCAAGa                                                            <  1:762413/150‑1 (MQ=255)
                                                                                             agataTAATCAGGGTCTTTTACAATACTGCCAGTTAAATTTTTT‑TTTTTTGTGGGTTATATTTCTTTAATAATAAACTTTGTGAGTTTTTATTTAAGGGCCCCGCCATTACCCCCTTCACGACTCTAAAAAATATCAGCACATCAAGAtt                                                          <  1:1631033/146‑1 (MQ=255)
                                                                                             aGCGATAATCACGGTCATTGACAATAATGCCATAGAAAATTTTC‑TTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTGTTACTGTGGATTTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGAtt                                                          <  1:3117810/150‑1 (MQ=255)
                                                                                               cGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                        >  2:2345871/1‑150 (MQ=255)
                                                                                               cGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                        >  2:2801342/1‑150 (MQ=255)
                                                                                                gATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGTATATTTCCTTAATAATACACTCTGTGACAGTGGAATTAACGGCCCCGCATTTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGt                                                       <  2:3422072/150‑1 (MQ=255)
                                                                                                gATAATCACGGTAATTGACAATAATGCCAGAGAAAATGTTC‑TTTTTTCTGGCGTATATTTCCTTAATAATAAACTCTGTGACAGTGTATTTAACGGCCCGGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGt                                                       <  2:2056217/150‑1 (MQ=255)
                                                                                                  tAACCACTGTCATTGACAATAAGGCCAGAGAAAATGTTA‑TTTTTTGTGGCGGAAAATTCATTAATAATAAACTCTGTGACAGTGTAATTAACGGCCCCGCAATTAACCCCTCAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCg                                                     <  1:2893836/150‑1 (MQ=255)
                                                                                                     tCACGGTCTTTTAAAATAATGCCATATAATTTTTTT‑TTTTTTCGGGCGTATAATTCATTAATATTAAACTGTGTGACTGTGTATTTAAGGGCCCCGCATTTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATa                                                  <  2:4258856/150‑1 (MQ=255)
                                                                                                     tCACGGTCTTTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCGGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAAGGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATa                                                  <  2:2177661/150‑1 (MQ=255)
                                                                                                      cACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAg                                                 >  1:3619693/1‑150 (MQ=255)
                                                                                                        cGGTCATTGACAATAATGCCAGAGAAAATGTTC‑TTTTTTGGGGCGTAAAATCCATTAATAATAAACTCTGTGACTGTGGAATTAACGGCCCCGCCATTAACCCCTGCACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGc                                               <  1:1163553/150‑1 (MQ=255)
                                                                                                        cGGACATTGACCATAATGCCAGAGAAATTTTTC‑TTTTTTCTGGCGTATAATTCCTTAATAATAAACTATGTGACAGTGTAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGc                                               <  2:3947446/150‑1 (MQ=255)
                                                                                                         ggTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCt                                              >  2:524509/1‑150 (MQ=255)
                                                                                                          gTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTa                                             >  1:2801376/1‑150 (MQ=255)
                                                                                                           tatttGACCATAATGCCATATCAATTGTTC‑TTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTTTGACTGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTac                                            <  1:2683669/147‑1 (MQ=255)
                                                                                                                gACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                       >  2:4047568/1‑150 (MQ=255)
                                                                                                                  ccATAATCCAAGACAAAATTTTC‑TTTTTTCTGGCGTATAATACATAAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCAATTAACCCCTCAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTaa                                     <  1:407935/150‑1 (MQ=255)
                                                                                                                   cATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCGACACATTAAt                                    >  2:498277/1‑150 (MQ=255)
                                                                                                                      acggCCCGAGAAATTGTTC‑TTTTTTGGGGCGGATAATTCCTTAATAATAAACTCTGTGACTGGGGAATTAACGGCCCCGCCTTTAACCCCGGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGaa                                 <  1:2040521/147‑1 (MQ=255)
                                                                                                                      aaggCCACAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGAAAGTGGAATTAACGGCCCCGCCATTAACCCCTTAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGaa                                 <  2:3135522/147‑1 (MQ=255)
                                                                                                                      aaTGCCTGAGAAATTGTTT‑TTTTTTCTGGCGTATAAGTCCTTAATATTAAACGAGGGGACAGTGTAATTAACGGCCCCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACCTCAAGATTGGTCGATAGGCTACACATTAATGaa                                 <  2:1517072/150‑1 (MQ=255)
                                                                                                                         gCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGAc                              >  1:2172072/1‑150 (MQ=255)
                                                                                                                          cctgagCAATTGTTC‑TTTTTTCTGGGGGATATTTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACa                             <  2:361142/147‑1 (MQ=255)
                                                                                                                           cAGAGCAATTGTTC‑TTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAg                            <  2:3219850/150‑1 (MQ=255)
                                                                                                                                     gTTA‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaaaa                  <  2:1717369/150‑1 (MQ=255)
                                                                                                                                      ttC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAg                 >  2:2847397/1‑150 (MQ=255)
                                                                                                                                      ttC‑TTTTTTCTAGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaaaat                 >  2:1811406/1‑149 (MQ=255)
                                                                                                                                       tC‑TTTTTTCTGGCTGATATTTCCTTAATAATAAACTTTGTGAATGTGTATTTAACGGCCCCGCATTTAACCCGTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGt                <  2:3693515/150‑1 (MQ=255)
                                                                                                                                        c‑TTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTTTTACTGTGTAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGC‑AAAAAAAAGTGa              <  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)
                                                                                                                                                      |                                                                                                                                                     
GACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTC‑TTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAA  >  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: ‑