Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,481,771 G→T 8.6% 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,7710GT 8.6% 1698.0 / 29.6 514P136T (CCG→ACG) emrYputative multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (260/189);  new base T (2/61);  total (262/252)
Fisher's exact test for biased strand distribution p-value = 1.55e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 527 total aligned reads displayed.

GACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAA  >  NC_000913/2481622‑2481920
                                                                                                                                                     |                                                                                                                                                     
gACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCgg                                                                                                                                                       >  1:2317703/1‑150 (MQ=255)
   gATAATCCCCATAGTGACATTGATTAAAATTATCAAACCCCAGCTAAAGATATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAAACACGGACAAAGACAATAATTCCACAGCAAATGTTCATTTATATCGCGGATa                                                                                                                                                    >  1:6262442/1‑150 (MQ=255)
     tAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAAt                                                                                                                                                  >  2:8135120/1‑150 (MQ=255)
     tAAACCCCATAGGGACATTGATTAAAAATATCCAACCCCATCTAAAGTTATCACAATCATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAAt                                                                                                                                                  <  2:8555321/150‑1 (MQ=255)
      aaTCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAAtt                                                                                                                                                 <  1:1019389/150‑1 (MQ=255)
      aaTCCCCATACCCACATTGAATAAAAATAACCAACCCCAGCTAAAGTAAACACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAAtt                                                                                                                                                 <  1:6788980/150‑1 (MQ=255)
         ccccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCt                                                                                                                                              <  1:1515282/150‑1 (MQ=255)
         ccccATAGGGACATTGATTAAAAATATACAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCt                                                                                                                                              >  2:7163213/1‑150 (MQ=255)
           ccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTta                                                                                                                                            <  1:5871726/150‑1 (MQ=255)
           ccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTta                                                                                                                                            >  2:167474/1‑150 (MQ=255)
           ccATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTta                                                                                                                                            <  2:2573761/150‑1 (MQ=255)
            cATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaa                                                                                                                                           <  1:2043762/150‑1 (MQ=255)
            cATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaa                                                                                                                                           <  2:4552363/150‑1 (MQ=255)
            cATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaa                                                                                                                                           <  2:1019951/150‑1 (MQ=255)
                gggACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaataat                                                                                                                                       <  1:370711/150‑1 (MQ=255)
                 ggACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaataata                                                                                                                                      >  1:1271080/1‑150 (MQ=255)
                   aCATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATaaa                                                                                                                                    >  2:3529565/1‑150 (MQ=255)
                   aCATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATaaa                                                                                                                                    <  2:1280355/150‑1 (MQ=255)
                       tGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCAAGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAActct                                                                                                                                >  1:1355568/1‑150 (MQ=255)
                           tAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGa                                                                                                                            >  1:5741603/1‑150 (MQ=255)
                             aaaaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCTCGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACa                                                                                                                          <  1:851069/150‑1 (MQ=255)
                             aaaaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACa                                                                                                                          <  1:3938703/150‑1 (MQ=255)
                             aaaaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACa                                                                                                                          >  2:336066/1‑150 (MQ=255)
                              aaaTATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAg                                                                                                                         >  2:411239/1‑150 (MQ=255)
                                atatCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTg                                                                                                                       <  1:6457743/150‑1 (MQ=255)
                                  atCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGa                                                                                                                     <  1:7364860/150‑1 (MQ=255)
                                     cAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAAtt                                                                                                                  >  1:8152688/1‑150 (MQ=255)
                                     cAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTCTGACAGTGGAAta                                                                                                                  >  1:5307378/1‑149 (MQ=255)
                                              cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGcccc                                                                                                         >  1:4541351/1‑150 (MQ=255)
                                              cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGcccc                                                                                                         >  1:7220163/1‑150 (MQ=255)
                                               tAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCg                                                                                                        >  2:5204703/1‑150 (MQ=255)
                                               tAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCg                                                                                                        >  1:5771293/1‑150 (MQ=255)
                                                aaaGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCANGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGNAATTAACGGCCCCGc                                                                                                       <  2:6758841/150‑1 (MQ=255)
                                                  aGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCa                                                                                                     >  2:3078178/1‑150 (MQ=255)
                                                      aTCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTaa                                                                                                 >  2:8274137/1‑150 (MQ=255)
                                                       tCACAAATATACCCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAc                                                                                                <  2:8380525/150‑1 (MQ=255)
                                                        cacaAATATACCCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCTTTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAcc                                                                                               <  1:103715/150‑1 (MQ=255)
                                                         acacATATACCCGCCCAAGCTCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAccc                                                                                              <  2:7037317/150‑1 (MQ=255)
                                                          caAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGTAATTAACGGCCCCGCCATTAAcccc                                                                                             <  1:661418/150‑1 (MQ=255)
                                                          caAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAcccc                                                                                             <  2:2449615/150‑1 (MQ=255)
                                                           aaaTATACCCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCt                                                                                            <  2:8082371/150‑1 (MQ=255)
                                                           aaaTATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCt                                                                                            >  1:6956628/1‑150 (MQ=255)
                                                           aaaTATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCt                                                                                            >  2:848687/1‑150 (MQ=255)
                                                              tataACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGATCATAATGCCATAGAAAATGTTCTTTTTTCGGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGaa                                                                                         <  1:5158043/150‑1 (MQ=255)
                                                              tataACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGaa                                                                                         <  2:610963/150‑1 (MQ=255)
                                                              tataACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGaa                                                                                         <  1:4951451/150‑1 (MQ=255)
                                                                  aCCGCCCAATATCGCCCCACATATCGGAGCTATAATCACGTTCTTTGACCATAATGCCATAGAAATTTTTATTTTTTCTGGCGTATATTTCCTTAATAATAAACTCTGTGACAGTGTAATTAACGGCCCCGCCATTAACCCCTGAACGAc                                                                                     <  1:2712561/150‑1 (MQ=255)
                                                                    cGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGActc                                                                                   >  1:8370334/1‑150 (MQ=255)
                                                                    cGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGActc                                                                                   >  1:2924274/1‑150 (MQ=255)
                                                                      cccAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATACTGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGACTTAACGGCCCCGCCATTAACCCCTGAACGACTCTa                                                                                 >  2:2373288/1‑150 (MQ=255)
                                                                          atatCGGCCCACATATCGTAGCGATATTCACGGTCATTTACCATAATGCCAGAGAAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaaaa                                                                             <  2:4921483/150‑1 (MQ=255)
                                                                          atatCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaaaa                                                                             >  2:5852461/1‑150 (MQ=255)
                                                                           tatCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATTTTCTTTTTTCTGGCGGATAATTCATTATTAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCTTTAACCCCTGAACGACTCTaaaaaa                                                                            <  2:3577471/150‑1 (MQ=255)
                                                                             tCGGCCCACTTATCGGAGCTATAATCACGGTCATTGACAATAATGCCATAGAAAATGTTCTTTTTTCTGGGGTATAATTCCTTAATAATAAACTCTGTGAATGTGGATTTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAata                                                                          <  2:188613/150‑1 (MQ=255)
                                                                               ggCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATc                                                                        >  2:4307420/1‑150 (MQ=255)
                                                                                   cacaTATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGca                                                                    <  2:4851575/150‑1 (MQ=255)
                                                                                    acaTATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcac                                                                   >  1:905037/1‑150 (MQ=255)
                                                                                     caTATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcaca                                                                  >  2:7236850/1‑150 (MQ=255)
                                                                                        atCGGTGCTATAATCACGTTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCa                                                               <  2:5073155/150‑1 (MQ=255)
                                                                                            gAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTACTTAATAATAAACTCTGTGACAGTGGAATTAAAGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATAAGCACAACAAGAt                                                           >  1:755864/1‑150 (MQ=255)
                                                                                              gCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTg                                                         >  2:574504/1‑150 (MQ=255)
                                                                                                taTAATCACGTTATTGGACAATATTGCAATAGAAATTTTTATTTTTTCTGGGGTATAATTCCTTAATAATAAACTCTGTTACTGTGTAATTAAGGGCCCCGCCATTAACCCCTGCACGACTCTAAAAAATATCAGCACATCAAGATTGGt                                                       <  1:4300095/149‑1 (MQ=255)
                                                                                                 aTAATCTCGTTCTTTGACCATATTGCCAGTGAATTTTTTTTTTTTTGTGGGGTTTATTTCCTTATTAATAAACTGTGTGACTTGGTAATTAAGGGCCCCGCCATTACCCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTc                                                      <  1:2433300/150‑1 (MQ=255)
                                                                                                     tCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATATTTCCTTAATAATAAACTCTGTGACAGTGGAATTAAGGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATa                                                  <  2:5564233/150‑1 (MQ=255)
                                                                                                     tCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTta                                                                                                                                            <  1:4825801/60‑1 (MQ=255)
                                                                                                     tCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATa                                                  >  1:5548005/1‑150 (MQ=255)
                                                                                                      cACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAg                                                 >  1:17769/1‑150 (MQ=255)
                                                                                                        cGGTCATTTACAATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGc                                               <  2:6286248/150‑1 (MQ=255)
                                                                                                            cATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTaca                                           >  1:408725/1‑150 (MQ=255)
                                                                                                            aaTTGACAATAATGCCAGAGCAATTGTTATTTTTTCGGGCGTATAATTCCTTAAAAATAAACCCTGTGACTTTGGAATTAACGGCCCCGCCATTAACCCCTCAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTaca                                           <  1:2161926/149‑1 (MQ=255)
                                                                                                               tGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAt                                        >  1:4598863/1‑150 (MQ=255)
                                                                                                                 aCCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTGTGACTGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTa                                      <  2:2827554/150‑1 (MQ=255)
                                                                                                                 aCCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTtaa                                                                                                                                           >  2:1704091/1‑49 (MQ=255)
                                                                                                                 aCCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACTGTGTAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTa                                      <  1:5001440/150‑1 (MQ=255)
                                                                                                                   cATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTGTGACTGTGGAATTAACGGCCCCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAAt                                    <  2:7871513/150‑1 (MQ=255)
                                                                                                                   cATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAAt                                    >  2:4361369/1‑150 (MQ=255)
                                                                                                                   cATAATCCCAGAGAAATTGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAAt                                    <  1:7712441/150‑1 (MQ=255)
                                                                                                                   aaTAATGCCAGAGAAAATGTTATTTTTTCTGGCGTATAATTCATTAATAATAAACTCTGTGACAGTGTATTTAACGGCCCCGCAATTAACCCCTTAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAAt                                    <  1:5649286/149‑1 (MQ=255)
                                                                                                                      aaTGCAAGAGCAAATGTTCTTTTTTCTGGGGTATAATTCCTTAATAATAAACTCTGTGACAGTGGATTTAACGGCCCCGCCTTTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGaa                                 <  2:8654618/150‑1 (MQ=255)
                                                                                                                          ccAGAGCAACTGTTCTTTTTTCTGGCGGTTAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACa                             <  2:1944486/150‑1 (MQ=255)
                                                                                                                          caagatCAATTGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGTAATTAACGGCCCCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACa                             <  2:5638226/148‑1 (MQ=255)
                                                                                                                            atagaaaaTTTTATTTTTTCTGGCGTATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCCCAATTAACCCCTCAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGc                           <  1:4710348/145‑1 (MQ=255)
                                                                                                                            agagCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGt                                                                                                                        >  2:4644728/1‑57 (MQ=255)
                                                                                                                             gagCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCg                          >  2:1373963/1‑150 (MQ=255)
                                                                                                                               gCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaa                        >  1:3340005/1‑150 (MQ=255)
                                                                                                                                cAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaa                       >  1:795523/1‑150 (MQ=255)
                                                                                                                                     tttattttttttGGGGTATAATTCATTAATAATAAACTCTGTGACTGTGTAATTAAGGGCCCCCCCATTAACCCCTTAACGAATATAAAAAATATCAGCACATCAAGTTTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaaaa                  <  1:3072654/146‑1 (MQ=255)
                                                                                                                                     tttCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaaaa                  <  2:2179285/149‑1 (MQ=255)
                                                                                                                                         ttttttCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGa              >  2:3981462/1‑150 (MQ=255)
                                                                                                                                         ttttttCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGa              >  2:7416493/1‑150 (MQ=255)
                                                                                                                                          tttttCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGAc             >  2:6484775/1‑150 (MQ=255)
                                                                                                                                          tttttCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGAc             >  2:5413192/1‑150 (MQ=255)
                                                                                                                                          tttttCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGAc             >  2:2057970/1‑150 (MQ=255)
                                                                                                                                          tttttCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGAc             >  1:1225305/1‑150 (MQ=255)
                                                                                                                                                tGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAg       >  2:3700381/1‑150 (MQ=255)
                                                                                                                                                 ggCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGt      >  2:7707712/1‑150 (MQ=255)
                                                                                                                                                   cGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaa    >  2:7546692/1‑150 (MQ=255)
                                                                                                                                                   cGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaa    >  1:6045665/1‑150 (MQ=255)
                                                                                                                                                    ggATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGGCAGCGAAAAAAAAGTGACTGAAAGTaaa   >  1:3508796/1‑150 (MQ=255)
                                                                                                                                                     gATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaaaa  >  1:1632989/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: