Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,519 A→C 7.8% I27M (ATT→ATG yadK ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,5190AC 7.8% 323.8 / 15.2 116I27M (ATT→ATGyadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (52/53);  new base C (7/2);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 120 total aligned reads displayed.

TTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTA  >  NC_000913/151370‑151668
                                                                                                                                                     |                                                                                                                                                     
tttCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCa                                                                                                                                                       >  1:1516472/1‑150 (MQ=255)
   ccACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATg                                                                                                                                                    <  1:1467716/150‑1 (MQ=255)
           ttGATAACTATCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                                                                            <  1:573/150‑1 (MQ=255)
            tGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTTCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAg                                                                                                                                           >  1:1915951/1‑150 (MQ=255)
             gATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgg                                                                                                                                          <  1:1205776/150‑1 (MQ=255)
              aTAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                                                                         >  2:631702/1‑150 (MQ=255)
              aTAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                                                                         <  2:1723733/150‑1 (MQ=255)
              aTAACTTTCAGGAGTAACATTACATGAAAAATATCCCACCCAAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGTTGCTGTTTTTTACATTGGAAACAAGGTCAATATCCTGTCCAGCAATGGCTGGAAAGGc                                                                                                                                         >  1:1169873/1‑150 (MQ=255)
                 aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCAAGCAATGGCTGGACAGGCGAt                                                                                                                                      >  1:1815411/1‑150 (MQ=255)
                  cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                     >  2:1275790/1‑150 (MQ=255)
                  cTTTCAGGAGTAACATTACATGAAAAATATCCCACCCCAACGCCGCCAAGATCAATATTACCCTGGTTACAGATTCCGCTTTTGCAGGAGCTGTTTTTTACAGTGGCAACAAGGTCAATATCCTGTCCAGCTATGGCTGGACACCCGATg                                                                                                                                     >  2:277319/1‑150 (MQ=255)
                   tttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGc                                                                                                                                    <  2:1711494/150‑1 (MQ=255)
                    ttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCa                                                                                                                                   >  1:583416/1‑150 (MQ=255)
                      cAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGCCGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTATTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                                                                                                 >  2:1423830/1‑150 (MQ=255)
                      cAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTAGCAGGTGCTGTTTTTTACATTGGCAACAACGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                                                                                                 <  1:1272528/150‑1 (MQ=255)
                        ggAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACAGATTCCGCTTTTGCAGGTGGTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCTATGGCTGGACAGGCGATGCAAAAc                                                                                                                               >  2:114755/1‑150 (MQ=255)
                                  tACCTGAAAAATATACCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGACAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGaca                                                                                                                     >  2:477957/1‑150 (MQ=255)
                                            atatCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAat                                                                                                           <  2:1517481/150‑1 (MQ=255)
                                               tCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGAAAGAAACAGAattat                                                                                                        >  2:1285919/1‑150 (MQ=255)
                                                 ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtc                                                                                                      >  2:1066949/1‑150 (MQ=255)
                                                  cACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGCAACAGAATTATtct                                                                                                     >  1:392867/1‑150 (MQ=255)
                                                    ccccAAAGAAGCCAAGATAAATATTACACTTGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                                                                                   <  1:827613/150‑1 (MQ=255)
                                                          cgacgCCAAGATCAATATTACACTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAAAAAGGTCAATATCATGTCCAGCAATGGCTGGACAGGCGATTCAAAACAGAAATGACAGAAACAGAATTAATATCTtcaaca                                                                                             >  1:1370878/1‑150 (MQ=255)
                                                            ccgCCACCATCAATATTACCCAGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTATTACATTGCCAACAAGGTCAATATCCTGGCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGc                                                                                           <  1:118020/149‑1 (MQ=255)
                                                             cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTATTGACAGAAACAGAATTATTCTCTTCATCAGCt                                                                                          >  1:742939/1‑150 (MQ=255)
                                                               ccAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTa                                                                                        >  1:1261119/1‑150 (MQ=255)
                                                                cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTACAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTAc                                                                                       >  1:1231829/1‑150 (MQ=255)
                                                                  aGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGt                                                                                     <  1:1794539/150‑1 (MQ=255)
                                                                   gATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt                                                                                    >  1:1771395/1‑150 (MQ=255)
                                                                     tCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGTTGCAAAACAGGAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGATGa                                                                                  >  2:328670/1‑150 (MQ=255)
                                                                     tCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                                                                  <  2:1263585/150‑1 (MQ=255)
                                                                       aaTATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                <  1:658383/150‑1 (MQ=255)
                                                                        atatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGAAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                                                                               <  2:1617833/150‑1 (MQ=255)
                                                                          atTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTggg                                                                             <  1:721440/150‑1 (MQ=255)
                                                                             aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATg                                                                          <  2:1148725/150‑1 (MQ=255)
                                                                             aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATg                                                                          <  2:851040/150‑1 (MQ=255)
                                                                             aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGAGGGATg                                                                          <  2:50377/150‑1 (MQ=255)
                                                                                cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATTACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                                                                       <  2:111646/150‑1 (MQ=255)
                                                                                cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                                                                       <  1:629123/150‑1 (MQ=255)
                                                                                cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                                                                       >  1:1263589/1‑150 (MQ=255)
                                                                                cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                                                                       <  1:1553800/150‑1 (MQ=255)
                                                                                      aCTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGAAAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAg                                                                 <  2:402868/150‑1 (MQ=255)
                                                                                       cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGt                                                                >  1:126821/1‑150 (MQ=255)
                                                                                        tGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGTAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                                                               >  2:1086036/1‑150 (MQ=255)
                                                                                         gactCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc                                                              >  1:1920928/4‑150 (MQ=255)
                                                                                          aTTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc                                                             >  2:1278427/1‑150 (MQ=255)
                                                                                           ttCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCt                                                            <  1:913365/150‑1 (MQ=255)
                                                                                            tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTg                                                           <  1:1120211/150‑1 (MQ=255)
                                                                                            tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTg                                                           <  1:867014/150‑1 (MQ=255)
                                                                                             ccGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa                                                          >  1:1417939/1‑150 (MQ=255)
                                                                                              cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa                                                         <  2:1698137/150‑1 (MQ=255)
                                                                                              cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa                                                         >  1:886043/1‑150 (MQ=255)
                                                                                              cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa                                                         <  1:713119/150‑1 (MQ=255)
                                                                                              cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATACTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa                                                         >  1:86665/1‑150 (MQ=255)
                                                                                              cGCTTTTGCAGGTGCTGTGTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTAACGTTGAGGGGGATGCATATATAGTTCCTGaa                                                         <  1:202465/150‑1 (MQ=255)
                                                                                                cTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATc                                                       >  1:1690015/1‑150 (MQ=255)
                                                                                                 ttttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCt                                                      <  2:1279837/150‑1 (MQ=255)
                                                                                                    tGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                   >  2:1445490/1‑150 (MQ=255)
                                                                                                    tGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTATCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                   <  1:1368480/150‑1 (MQ=255)
                                                                                                             tGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttt                                          >  1:600893/1‑149 (MQ=255)
                                                                                                              gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGTACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttag                                         <  2:141521/150‑1 (MQ=255)
                                                                                                                tttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttagtt                                       >  1:1612382/1‑150 (MQ=255)
                                                                                                                 ttttACATTGGCAACCAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGATTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttagata                                      >  1:438845/1‑150 (MQ=255)
                                                                                                                      cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTTTAGGt                                 <  1:238848/150‑1 (MQ=255)
                                                                                                                      cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGt                                 >  2:1144892/1‑150 (MQ=255)
                                                                                                                         tGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTAc                              <  2:405081/150‑1 (MQ=255)
                                                                                                                         tGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTAc                              <  2:811001/150‑1 (MQ=255)
                                                                                                                          ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACt                             >  1:1712876/1‑150 (MQ=255)
                                                                                                                          ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACt                             <  2:1023008/150‑1 (MQ=255)
                                                                                                                           gCAACAAGGTCAATATCCTTTCCAGCAATGGCTGCACAGGCGATGCAAAACAGTAATGACAGAACCAGAATTATTCTCTTCATCTGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTg                            <  2:290852/150‑1 (MQ=255)
                                                                                                                           gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTg                            >  1:1366004/1‑150 (MQ=255)
                                                                                                                           gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTg                            <  2:854527/150‑1 (MQ=255)
                                                                                                                           gCAACAAGGTCAATATCCTGTCCAGAAATGGCTGGGCAGGCGATGCTACACAGTAATGTCAGAAAGAGAAATATTCTCTTCTTCAGCTGACCTTGAGGGGGATGCAGTTATAGTTCCTGAATCTGTAAGTCTGATAGTGATAGGTTACTg                            >  2:370689/1‑150 (MQ=255)
                                                                                                                             aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTTTTCTCTTCATCAGCTTACGTTGCGTGGGATTCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTa                          <  2:1580852/150‑1 (MQ=255)
                                                                                                                                  ggTCAATATCCTGTCCAGCAATGGCTTGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtgt                     <  1:631702/150‑1 (MQ=255)
                                                                                                                                     cAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAATCAGATTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                  >  1:1818078/1‑150 (MQ=255)
                                                                                                                                     cAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                  >  2:540159/1‑150 (MQ=255)
                                                                                                                                     aaaTATCCTGTCCAGCAATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                  <  2:616665/149‑1 (MQ=255)
                                                                                                                                     aaaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                  <  1:1275790/149‑1 (MQ=255)
                                                                                                                                     aaaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                  <  1:504661/149‑1 (MQ=255)
                                                                                                                                      aaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCt                 <  1:744811/150‑1 (MQ=255)
                                                                                                                                      aaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCt                 >  1:1643737/1‑150 (MQ=255)
                                                                                                                                         atCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGt              >  1:785410/1‑150 (MQ=255)
                                                                                                                                           ccTGTCTAGAAATGGCTGGACAGGCGATGCAAAACAGAAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCTTATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATTGGTTACTGTAAATGTCGCTGGTGc            <  2:94921/150‑1 (MQ=255)
                                                                                                                                           ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTgt            >  2:85787/1‑149 (MQ=255)
                                                                                                                                           ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGc            >  1:254508/1‑150 (MQ=255)
                                                                                                                                           ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGc            >  2:1440792/1‑150 (MQ=255)
                                                                                                                                           ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGc            >  2:522371/1‑150 (MQ=255)
                                                                                                                                           ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGATGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTCTAAATGTGGCTGGTGc            >  1:751182/1‑150 (MQ=255)
                                                                                                                                            cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCa           >  1:591852/1‑150 (MQ=255)
                                                                                                                                            cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCa           >  1:1712805/1‑150 (MQ=255)
                                                                                                                                            cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCa           >  1:860905/1‑150 (MQ=255)
                                                                                                                                             tGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAg          >  2:200423/1‑150 (MQ=255)
                                                                                                                                             tGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAg          >  1:661172/1‑150 (MQ=255)
                                                                                                                                               tCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaa        >  1:1545299/1‑150 (MQ=255)
                                                                                                                                                 cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCGGTAAGTCTGGTAGTTATAGGTTACTGTATATGTGGCTGGTGCAGacaa      >  2:1191053/1‑150 (MQ=255)
                                                                                                                                                  aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc     <  2:1062920/150‑1 (MQ=255)
                                                                                                                                                  aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGATCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc     <  1:271869/150‑1 (MQ=255)
                                                                                                                                                     aaTGGCTGGACAGTCTATTCAATACATTAATGACAGAATCATAATTATTCTCTTCATCAGCTTACGTTGAGTGTGATGCATATATAGTTCCTGAATCTTTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGTTGCAGAAAACTTa  <  2:203345/150‑1 (MQ=255)
                                                                                                                                                     aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTa  <  2:1958950/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTA  >  NC_000913/151370‑151668

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