Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,522 G→T 8.2% A26A (GCC→GCA yadK ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,5220GT 8.2% 299.0 / 16.9 118A26A (GCC→GCAyadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (54/54);  new base T (7/3);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.01e-01

Only 100 of 119 total aligned reads displayed.

CCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTA  >  NC_000913/151373‑151668
                                                                                                                                                     |                                                                                                                                                  
ccACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATg                                                                                                                                                    <  1:1467716/150‑1 (MQ=255)
   ccTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCt                                                                                                                                                 <  1:1294132/150‑1 (MQ=255)
        ttGATAACTATCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                                                                            <  1:573/150‑1 (MQ=255)
          gATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgg                                                                                                                                          <  1:1205776/150‑1 (MQ=255)
           aTAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                                                                         >  2:631702/1‑150 (MQ=255)
           aTAACTTTCAGGAGTAACATTACATGAAAAATATCCCACCCAAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGTTGCTGTTTTTTACATTGGAAACAAGGTCAATATCCTGTCCAGCAATGGCTGGAAAGGc                                                                                                                                         >  1:1169873/1‑150 (MQ=255)
              aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                                                                      <  2:1333533/150‑1 (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)
                   cAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGCCGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTATTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                                                                                                 >  2:1423830/1‑150 (MQ=255)
                   cAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTAGCAGGTGCTGTTTTTTACATTGGCAACAACGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                                                                                                 <  1:1272528/150‑1 (MQ=255)
                     ggAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACAGATTCCGCTTTTGCAGGTGGTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCTATGGCTGGACAGGCGATGCAAAAc                                                                                                                               >  2:114755/1‑150 (MQ=255)
                       aGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                                                                             <  1:90372/150‑1 (MQ=255)
                          aaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                                                                                                          >  1:796932/1‑150 (MQ=255)
                               tACCTGAAAAATATACCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGACAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGaca                                                                                                                     >  2:477957/1‑150 (MQ=255)
                                      aaaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacag                                                                                                              >  2:1090246/1‑150 (MQ=255)
                                             cccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatt                                                                                                       <  1:624269/150‑1 (MQ=255)
                                               cACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGCAACAGAATTATtct                                                                                                     >  1:392867/1‑150 (MQ=255)
                                                 ccccAAAGAAGCCAAGATAAATATTACACTTGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                                                                                   <  1:827613/150‑1 (MQ=255)
                                                 caccAACGACGCCTAGATCTATATTTCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                                                                                   <  2:1814056/148‑1 (MQ=255)
                                                   ccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                                                                                                 >  1:1659790/1‑150 (MQ=255)
                                                       cgacgCCAAGATCAATATTACACTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAAAAAGGTCAATATCATGTCCAGCAATGGCTGGACAGGCGATTCAAAACAGAAATGACAGAAACAGAATTAATATCTtcaaca                                                                                             >  1:1370878/1‑150 (MQ=255)
                                                         ccgCCACCATCAATATTACCCAGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTATTACATTGCCAACAAGGTCAATATCCTGGCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGc                                                                                           <  1:118020/149‑1 (MQ=255)
                                                          cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTATTGACAGAAACAGAATTATTCTCTTCATCAGCt                                                                                          >  1:742939/1‑150 (MQ=255)
                                                             cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTACAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTAc                                                                                       >  1:1231829/1‑150 (MQ=255)
                                                              aaGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACg                                                                                      <  2:1520118/150‑1 (MQ=255)
                                                              aaGATCAATATTACCCTGGTTACTGATTCCGCTCATGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACg                                                                                      <  1:1501067/150‑1 (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)
                                                                           cccTGGTTACTGATTCCGCTTTTCCAGGTGCTGTTTTTTACAAAGGCAACAAGGTAAATATCCTGTCCAGCAATGGCTGGACACGCGATGAAAAACAGTAATGACAGAAACAGAATTATTATCTTAATCACCTTACGTTGAGTGGGATga                                                                         >  1:203345/1‑149 (MQ=255)
                                                                             cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATTACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                                                                       <  2:111646/150‑1 (MQ=255)
                                                                             cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                                                                       <  1:629123/150‑1 (MQ=255)
                                                                             cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                                                                       <  1:1553800/150‑1 (MQ=255)
                                                                             cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCat                                                                       >  1:1263589/1‑150 (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)
                                                                                          ccGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa                                                          >  1:1417939/1‑150 (MQ=255)
                                                                                           cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa                                                         <  1:452854/150‑1 (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                             <  2:1023008/150‑1 (MQ=255)
                                                                                                                        gCAACAAGGTCAATATCCTTTCCAGCAATGGCTGCACAGGCGATGCAAAACAGTAATGACAGAACCAGAATTATTCTCTTCATCTGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTg                            <  2:290852/150‑1 (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)
                                                                                                                               ggTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtgt                     >  1:716493/1‑150 (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:504661/149‑1 (MQ=255)
                                                                                                                                  aaaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                  <  1:1275790/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)
                                                                                                                                         cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCa           >  1:591852/1‑150 (MQ=255)
                                                                                                                                         cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCa           >  1:1712805/1‑150 (MQ=255)
                                                                                                                                         cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGCAACAGAATTATGCTCTTCATCAGCTGACGTTGAGGGGGATGCATATATAGTTCCTGAAGCTGTAAGTCTGTTAGTTCTAGGTTACTGTAAATGTGGCTGGTGCa           <  1:1799133/150‑1 (MQ=255)
                                                                                                                                         cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCa           >  1:1246796/1‑150 (MQ=255)
                                                                                                                                         cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCa           >  1:860905/1‑150 (MQ=255)
                                                                                                                                          tGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAg          >  2:200423/1‑150 (MQ=255)
                                                                                                                                               aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc     <  2:1062920/150‑1 (MQ=255)
                                                                                                                                                  aaTGGCTGGACAGTCTATTCAATACATTAATGACAGAATCATAATTATTCTCTTCATCAGCTTACGTTGAGTGTGATGCATATATAGTTCCTGAATCTTTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGTTGCAGAAAACTTa  <  2:203345/150‑1 (MQ=255)
                                                                                                                                                  aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTa  <  2:1958950/150‑1 (MQ=255)
                                                                                                                                                  aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTa  <  2:1261119/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTA  >  NC_000913/151373‑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: