Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,522 G→T 9.0% 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 9.0% 320.2 / 15.6 116A26A (GCC→GCAyadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (49/56);  new base T (4/7);  total (53/63)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 119 total aligned reads displayed.

CCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCC  >  NC_000913/151373‑151671
                                                                                                                                                     |                                                                                                                                                     
   ccTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTAACCTGGTTACTGATTCCGCTTTTGCAGGTGCAGTTTTTTACATTGGAAACAAGGTCAATATACTGTCAAGCAATggat                                                                                                                                                    >  1:838339/1‑148 (MQ=255)
    cTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                                                                                                   <  2:1131277/150‑1 (MQ=255)
    cTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                                                                                                   >  2:1482223/1‑150 (MQ=255)
     tGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTgg                                                                                                                                                  >  1:992247/1‑150 (MQ=255)
      ggTTGATAACTTTCAGGAGTACCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGa                                                                                                                                                 >  2:1183063/1‑150 (MQ=255)
        ttGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCAAAGATCAATATAACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                                                                               >  1:1170050/1‑150 (MQ=255)
         tGATAACTTTCAGGAGTAACATTACGTGAATAATATGCCACCCCAACGACGCCAAGATTAATATTCCCCTGGTTACTGATTCCGCTTTTGCAGTTGCTGTTTTTTTCATTGGCAACATGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                                                                              <  1:305059/150‑1 (MQ=255)
         tGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGCTTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGAAATGGCTGGACAg                                                                                                                                              >  2:335006/1‑150 (MQ=255)
         tGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                                                                              <  2:609090/150‑1 (MQ=255)
         tGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAg                                                                                                                                              >  2:515809/1‑150 (MQ=255)
          taTAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgg                                                                                                                                             <  1:1565746/149‑1 (MQ=255)
          gATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgg                                                                                                                                             <  1:823720/150‑1 (MQ=255)
             aaCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGa                                                                                                                                          >  1:117254/1‑150 (MQ=255)
              aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                                                                         >  1:1711414/1‑150 (MQ=255)
              aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                                                                         >  2:190449/1‑150 (MQ=255)
              aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                                                                         <  2:456077/150‑1 (MQ=255)
               cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                        <  1:1059883/150‑1 (MQ=255)
                   cAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                                                                                                    >  1:1094519/1‑150 (MQ=255)
                    aGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaaa                                                                                                                                   >  2:1173657/1‑150 (MQ=255)
                     ggAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAAc                                                                                                                                  >  2:1442215/1‑150 (MQ=255)
                     ggAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAAc                                                                                                                                  >  1:1515625/1‑150 (MQ=255)
                     ggAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAAc                                                                                                                                  >  2:958529/1‑150 (MQ=255)
                      gAGTAACATTACCTGAAAAATATCCCACCCCAACTACGCCAAGATCAAAATTACCATTGTTACTAATTCCGCTTTTGAAGGTGCTGTTTTTTACATTTGCTACAAGGTCAATATCATGTACAGCAATTGCGGGACAGGCGATGCAAAACa                                                                                                                                 >  1:1068138/1‑150 (MQ=255)
                       aGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                                                                                >  1:566704/1‑150 (MQ=255)
                       aGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTACAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAg                                                                                                                                <  2:660045/150‑1 (MQ=255)
                        gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGt                                                                                                                               >  2:1253618/1‑150 (MQ=255)
                        gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGt                                                                                                                               >  2:1260176/1‑150 (MQ=255)
                        gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACACTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGt                                                                                                                               <  2:836829/150‑1 (MQ=255)
                         taaccatACCTGAAAAAAATCCCTCCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTATGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCAGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTa                                                                                                                              <  1:283140/144‑1 (MQ=255)
                          aaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                                                                                                             >  1:885184/1‑150 (MQ=255)
                           aCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAAt                                                                                                                            >  2:774368/1‑150 (MQ=255)
                           aCAATACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAAt                                                                                                                            <  2:1155122/150‑1 (MQ=255)
                            cATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATg                                                                                                                           <  2:678245/150‑1 (MQ=255)
                               tACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGaca                                                                                                                        <  1:1260176/150‑1 (MQ=255)
                               tACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGaca                                                                                                                        <  1:1253618/150‑1 (MQ=255)
                                 ccTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacaga                                                                                                                      <  2:967201/150‑1 (MQ=255)
                                   tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                                                                    >  1:385282/1‑150 (MQ=255)
                                   tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                                                                    >  1:1437797/1‑150 (MQ=255)
                                    gAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaac                                                                                                                   <  2:1718630/150‑1 (MQ=255)
                                    gAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaac                                                                                                                   <  2:686496/150‑1 (MQ=255)
                                            tCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattat                                                                                                           >  2:1172357/1‑150 (MQ=255)
                                            tCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattat                                                                                                           >  2:960415/1‑150 (MQ=255)
                                             cccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatt                                                                                                          >  2:1117805/1‑150 (MQ=255)
                                             cccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatt                                                                                                          >  1:566439/1‑150 (MQ=255)
                                              ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtc                                                                                                         >  2:827743/1‑150 (MQ=255)
                                              caaCCCCAACGACTCCAAGATCAATATTACCCTGGTTACTGATTCAGCTTTTGCATGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtc                                                                                                         <  2:330311/148‑1 (MQ=255)
                                                aCCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtctc                                                                                                       <  1:1173657/150‑1 (MQ=255)
                                                   ccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGCCAGCAACAGAATTATTCTCTtc                                                                                                    >  2:328348/1‑150 (MQ=255)
                                                   ccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                                                                                                    >  2:1029394/1‑150 (MQ=255)
                                                    cAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                                   <  2:941135/150‑1 (MQ=255)
                                                    cAACGACGCCAAGATCAATATTACCCTGGTTACTGATACCGCTTTTGCAGGAGCTCTTATTTACATTGGCATCATGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATAATTCTCTtca                                                                                                   >  1:1159530/1‑150 (MQ=255)
                                                             aaaGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTTTTTTTTTCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACATAATTATTCTCTTCATCAGCTTAc                                                                                          <  1:1278722/149‑1 (MQ=255)
                                                              aaGATCAATATTAACCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACg                                                                                         <  2:817732/150‑1 (MQ=255)
                                                               aGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGt                                                                                        <  2:1094519/150‑1 (MQ=255)
                                                                gATCTATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGTAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt                                                                                       >  2:705948/1‑150 (MQ=255)
                                                                  tCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                                                                     <  2:934758/150‑1 (MQ=255)
                                                                      tatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTgg                                                                                 <  1:1634386/150‑1 (MQ=255)
                                                                          aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATg                                                                             <  1:515809/150‑1 (MQ=255)
                                                                           cccTGGTTACTGATTCCTCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGGCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGc                                                                            <  2:853790/150‑1 (MQ=255)
                                                                            ccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                                                                           <  2:603947/150‑1 (MQ=255)
                                                                              tGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCata                                                                         <  1:958529/150‑1 (MQ=255)
                                                                                    cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGt                                                                   >  1:119242/1‑150 (MQ=255)
                                                                                    cTGAGTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTACAGCAATGGCTGGACAGGCTATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATAAGCTTACGTTGATTTTGATGCATATATAGt                                                                   >  2:158718/1‑150 (MQ=255)
                                                                                      gATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCATTATCCTGTGCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGTCAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtta                                                                 <  2:694023/150‑2 (MQ=255)
                                                                                      gATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc                                                                 >  1:1002689/1‑150 (MQ=255)
                                                                                      gATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTACAGCAATGTCTGGACAGGCGATGCAAAAAAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc                                                                 >  1:727677/1‑150 (MQ=255)
                                                                                       aTTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGATTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc                                                                <  1:190449/150‑1 (MQ=255)
                                                                                        ttCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCt                                                               >  1:1114096/1‑150 (MQ=255)
                                                                                            gCTTTTTCAGGTTCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATTGCTGGACAGGCGATGCAAAACAGTAATGACTGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAt                                                           <  1:600065/150‑1 (MQ=255)
                                                                                              ttttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCt                                                         <  2:178558/150‑1 (MQ=255)
                                                                                               tttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTg                                                        <  1:1414496/150‑1 (MQ=255)
                                                                                               tttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTg                                                        <  1:74073/150‑1 (MQ=255)
                                                                                                 tGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                      >  1:1590015/1‑150 (MQ=255)
                                                                                                    gggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGt                                                   >  2:477390/2‑150 (MQ=255)
                                                                                                     ggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATTCAAAACATTAATGACAGACACAGAAATATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATATGTAAGta                                                  >  2:120846/1‑149 (MQ=255)
                                                                                                     ggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTc                                                  <  2:1711414/150‑1 (MQ=255)
                                                                                                      gTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGATGGACAGGCGATGCAAAACAGTAATGACAGATACAGAAATATTCTATTCATCAGCTTACGTTGATTGGGATGCATATATATTTCCTGAATCTGTAAGTCt                                                 >  2:1136848/1‑150 (MQ=255)
                                                                                                       tGCTGTTTTTTACTTTGTCAACATGGTCAATTTCCTGTCCAGCAATGGCTGGACAGGGGATGCAATACAGTATTGAGTGTAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTg                                                <  2:1159530/150‑1 (MQ=255)
                                                                                                         cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtt                                              <  1:960415/150‑1 (MQ=255)
                                                                                                         cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtt                                              <  2:1053642/150‑1 (MQ=255)
                                                                                                         cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtt                                              >  1:611781/1‑150 (MQ=255)
                                                                                                               tttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtat                                        >  1:591394/1‑150 (MQ=255)
                                                                                                               tttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAAAAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGGtat                                        >  2:56236/1‑150 (MQ=255)
                                                                                                                 tACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAg                                      <  2:385282/150‑1 (MQ=255)
                                                                                                                  aCATTGGCAACAAGGTCAATATCCTGACCAGAAATGGCTTGACAGGCGATGCAAAACAGGAATGACAGAAACATAATTATTCTCTACATCAGCTTACGTTGAGTGGGATGCATATTTAGTTCCTGAATCTGTAAGTCTGTTAGTTATAgg                                     <  2:116988/150‑1 (MQ=255)
                                                                                                                   cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGt                                    <  2:207057/150‑1 (MQ=255)
                                                                                                                      tGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTAc                                 <  1:1482223/150‑1 (MQ=255)
                                                                                                                      tGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTAc                                 <  2:1333413/150‑1 (MQ=255)
                                                                                                                          aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTa                             <  1:1034925/150‑1 (MQ=255)
                                                                                                                           acaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTaa                            <  2:1170050/150‑1 (MQ=255)
                                                                                                                               ggTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtgt                        >  1:1673227/1‑150 (MQ=255)
                                                                                                                                  aaaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGTGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                     <  1:1656193/149‑1 (MQ=255)
                                                                                                                                   catcgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCt                    >  1:1451825/6‑150 (MQ=255)
                                                                                                                                   aaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCt                    <  2:1515625/150‑1 (MQ=255)
                                                                                                                                    atcgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTACTGAATCTGTAAGTCTGTTAGTTAGAGGTTAGTGTAAATGTGGCTg                   >  2:1293569/5‑150 (MQ=255)
                                                                                                                                     tatCCTGTCAAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTgg                  <  2:1013086/150‑1 (MQ=255)
                                                                                                                                       tCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTg                >  1:1565403/1‑150 (MQ=255)
                                                                                                                                       tCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTg                >  2:603823/1‑150 (MQ=255)
                                                                                                                                                  aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTa     <  1:477390/150‑1 (MQ=255)
                                                                                                                                                    tGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTAcc   <  1:705948/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCC  >  NC_000913/151373‑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: