Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,522 G→T 5.5% 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 5.5% 403.1 / 9.9 136A26A (GCC→GCA) yadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (62/65);  new base T (3/5);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

Only 100 of 137 total aligned reads displayed.

CTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCC  >  NC_000913/151377‑151671
                                                                                                                                                 |                                                                                                                                                     
cTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                                                                                                   <  2:1183940/150‑1 (MQ=255)
   gTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGAc                                                                                                                                                <  2:1536710/150‑1 (MQ=255)
     tGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAATATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                                                                              <  1:1642579/150‑1 (MQ=255)
     tGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                                                                              <  1:450714/150‑1 (MQ=255)
     tGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAg                                                                                                                                              <  1:1349065/150‑1 (MQ=255)
      gATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgg                                                                                                                                             <  1:1910079/150‑1 (MQ=255)
          aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                                                                         >  1:1499331/1‑150 (MQ=255)
          aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTATGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                                                                         >  1:1874722/1‑150 (MQ=255)
           cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                        <  1:54617/150‑1 (MQ=255)
           cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                        >  2:1877283/1‑150 (MQ=255)
           cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                        >  2:1376222/1‑150 (MQ=255)
            tttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCTATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGc                                                                                                                                       <  1:263887/150‑1 (MQ=255)
              tCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaa                                                                                                                                     <  2:939066/150‑1 (MQ=255)
              tCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaa                                                                                                                                     >  2:608724/1‑150 (MQ=255)
                aGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCGGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaaa                                                                                                                                   >  2:1958867/1‑150 (MQ=255)
                 ggAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAAc                                                                                                                                  <  1:10287/150‑1 (MQ=255)
                 ggAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCAAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGAGCTGTTTGTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaataa                                                                                                                                  >  1:1456028/1‑147 (MQ=255)
                    gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGt                                                                                                                               >  1:1847829/1‑150 (MQ=255)
                    gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGt                                                                                                                               <  2:1750990/150‑1 (MQ=255)
                      aaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                                                                                                             >  1:119109/1‑150 (MQ=255)
                        cATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATg                                                                                                                           <  2:1254415/150‑1 (MQ=255)
                              cTGAATAATATCCGACCCCAACGTCGCCTACATCAATATTTGGGTGGTTACTGTTTCGGCTTTTGGTGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATGCTGTGAAGCAATGGCTGGTCTGGCGATGCAAATCAGTAATGtcagta                                                                                                                     <  1:783543/150‑3 (MQ=255)
                              cTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATTCAAAACAGTAATGacagaa                                                                                                                     <  1:1468005/150‑1 (MQ=255)
                               tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGATTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                                                                    >  2:1780818/1‑150 (MQ=255)
                                 aaaaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaaca                                                                                                                  <  2:870170/150‑1 (MQ=255)
                                 aaaaaTATCCCACCCCAACGACGACAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTTTTTTTTACATTTGAAACAAGGTCAATATCCTGTCCAGGTATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaaca                                                                                                                  <  1:7043/150‑1 (MQ=255)
                                  aaaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGCCAGGCGATGCAAAACAGTAATGacagaaacag                                                                                                                 >  2:655781/1‑150 (MQ=255)
                                  aaaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacag                                                                                                                 <  2:1116165/150‑1 (MQ=255)
                                    aaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                                                                               <  2:1837770/150‑1 (MQ=255)
                                    aaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                                                                               >  1:806132/1‑150 (MQ=255)
                                    aaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                                                                               <  2:1392749/150‑1 (MQ=255)
                                       atCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGCCAGAAACAGAatta                                                                                                            >  2:1295663/1‑150 (MQ=255)
                                         cccACCCCAACGACGCAAAGATCAATATTACCCTGGTTAATGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAAGAAAAACAGTAATGACAGAAACATAattatt                                                                                                          >  2:271490/1‑150 (MQ=255)
                                          ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtc                                                                                                         >  2:490071/1‑150 (MQ=255)
                                          ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtc                                                                                                         <  2:433520/150‑1 (MQ=255)
                                          ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCAAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGAAAGAAACAGAattatta                                                                                                         >  2:457413/1‑149 (MQ=255)
                                           cACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtct                                                                                                        >  1:1580274/1‑150 (MQ=255)
                                             ccccAACGCCGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                                                                                      >  2:1047392/1‑150 (MQ=255)
                                              cccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCtt                                                                                                     <  2:1472197/150‑1 (MQ=255)
                                              cccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCtt                                                                                                     >  1:1660279/1‑150 (MQ=255)
                                               ccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTCTACAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                                                                                                    >  1:330963/1‑150 (MQ=255)
                                                     acgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGc                                                                                              <  2:1164622/150‑1 (MQ=255)
                                                      cgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCt                                                                                             >  2:957176/1‑150 (MQ=255)
                                                        ccAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTa                                                                                           >  2:1948895/1‑150 (MQ=255)
                                                          aaGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTAAATATCCTGTCCTGCAATGGCTGGACAGGCGATGAAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACg                                                                                         >  2:303197/1‑150 (MQ=255)
                                                            gATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt                                                                                       >  1:636489/1‑150 (MQ=255)
                                                                 atatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                                                                                  >  1:1174567/1‑150 (MQ=255)
                                                                 atatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                                                                                  <  2:2045790/150‑1 (MQ=255)
                                                                  tatTACCCTGGGTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTgg                                                                                 <  1:1545628/150‑1 (MQ=255)
                                                                    ttACCCTGGTTACTGATTCCGCTTTTGCAGGTGCAGTTTTTTACATTGGCAACAAGGTCAAAATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                               <  2:91111/150‑1 (MQ=255)
                                                                       cccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGc                                                                            >  2:762515/1‑150 (MQ=255)
                                                                       cccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGc                                                                            >  2:1610620/1‑150 (MQ=255)
                                                                         cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGctt                                                                          >  2:121994/1‑148 (MQ=255)
                                                                          tGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCata                                                                         <  1:1717923/150‑1 (MQ=255)
                                                                          tGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCata                                                                         >  2:1046012/1‑150 (MQ=255)
                                                                          tGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGAAATGACAGAAACAGAATAATTCTCTTCATCAGCTTACGTTGAGTGGGATGCata                                                                         >  1:237415/1‑150 (MQ=255)
                                                                              tACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCAtatata                                                                     >  1:180873/1‑150 (MQ=255)
                                                                                    ttCCGCTTTTTCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCt                                                               <  2:288423/150‑1 (MQ=255)
                                                                                       cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa                                                            >  2:2038165/1‑150 (MQ=255)
                                                                                       cGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGaa                                                            >  2:673684/1‑150 (MQ=255)
                                                                                         cTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATc                                                          <  1:1877284/150‑1 (MQ=255)
                                                                                           tttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTACAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTATCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTACTGAATCtt                                                        >  2:1453773/1‑149 (MQ=255)
                                                                                            ttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGt                                                       <  2:1989524/150‑1 (MQ=255)
                                                                                             tGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                      >  2:1124579/1‑150 (MQ=255)
                                                                                             tGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                      <  2:1847828/150‑1 (MQ=255)
                                                                                               cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAg                                                    <  1:608724/150‑1 (MQ=255)
                                                                                                aGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGt                                                   >  2:1069594/1‑150 (MQ=255)
                                                                                                 ggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCGTGTCCTGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTc                                                  <  1:1778482/150‑1 (MQ=255)
                                                                                                 ggTGCTGTTATTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGta                                                  >  1:1707388/1‑149 (MQ=255)
                                                                                                   tGCTGTTTTTTACATTGGCAACAACGTCAATATCCTGTCCAGCAATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTg                                                <  2:1110585/150‑1 (MQ=255)
                                                                                                    gCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATCGCTGGACAGGCGATGAAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATAAGCTTACCTTGAGTGGAATGCATATATAGTTCATGAATCTGTAAGTCTga                                               >  1:283650/1‑149 (MQ=255)
                                                                                                     cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtt                                              <  2:1729548/150‑1 (MQ=255)
                                                                                                     cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtt                                              >  1:465049/1‑150 (MQ=255)
                                                                                                      tGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtta                                             <  1:1318679/150‑1 (MQ=255)
                                                                                                           tttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtat                                        <  1:1893314/150‑1 (MQ=255)
                                                                                                             tACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAg                                      >  2:510801/1‑150 (MQ=255)
                                                                                                                aTTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGtt                                   >  1:433883/1‑150 (MQ=255)
                                                                                                                aTTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGAATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGtt                                   >  2:1812453/1‑150 (MQ=255)
                                                                                                                   ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACt                                >  1:1548974/1‑150 (MQ=255)
                                                                                                                    gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGTCAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCCGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTg                               <  1:1047393/150‑1 (MQ=255)
                                                                                                                    gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTTTTAGTTATAGGTTaatt                               >  2:583084/1‑147 (MQ=255)
                                                                                                                    gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTg                               <  1:490071/150‑1 (MQ=255)
                                                                                                                    gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTg                               <  2:546035/150‑1 (MQ=255)
                                                                                                                      aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTa                             <  2:1874721/150‑1 (MQ=255)
                                                                                                                         aaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAt                          <  1:1756591/150‑1 (MQ=255)
                                                                                                                              aaaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGCCAGAAACAGAATTAATCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCAGTTAGTTATAGGTTACTGTAAATGTGGc                     <  1:1124577/149‑1 (MQ=255)
                                                                                                                                   tCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCGGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTg                >  1:718634/1‑150 (MQ=255)
                                                                                                                                   tCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTACTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTg                >  2:831288/1‑150 (MQ=255)
                                                                                                                                    ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGc               >  1:950190/1‑150 (MQ=255)
                                                                                                                                    ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAGGTGGCTGGTGc               <  1:1905587/150‑1 (MQ=255)
                                                                                                                                     cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCa              >  1:663594/1‑150 (MQ=255)
                                                                                                                                     cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGTTGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTTTAGGTTACTGTAAATGTGGCTGGTGCa              >  2:1657694/1‑150 (MQ=255)
                                                                                                                                      tGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAg             <  1:1610621/150‑1 (MQ=255)
                                                                                                                                         ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGACTCTGTTAGTCTGTTAGTTATAGGTGACTGGTAATGTGGCTGGTGCAGaaa          >  2:1437574/1‑150 (MQ=255)
                                                                                                                                         ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaa          >  1:347714/1‑150 (MQ=255)
                                                                                                                                         ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaa          >  2:611922/1‑150 (MQ=255)
                                                                                                                                           aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc        <  1:914146/150‑1 (MQ=255)
                                                                                                                                           aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACACAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc        <  1:938255/150‑1 (MQ=255)
                                                                                                                                            gCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACTGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACt       <  1:1069595/150‑1 (MQ=255)
                                                                                                                                              aaTGGCTGGACAGGCGATGCAATACAGTAATTACAGAAAAAGAATTATTCTCTTCATCAGCTTACGTTGAGTGTTATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTAATGTAAATGTGGCTGGTGCAGAAAACTTa     <  1:744964/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCC  >  NC_000913/151377‑151671

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