Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 858,425 G→T 5.2% L213M (CTG→ATG) ‡ ybiU ← DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913858,4250GT 5.2% 326.0 / 4.3 104L213M (CTG→ATG) ‡ybiUDUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base G (50/48);  new base T (2/4);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

GCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGAT  >  NC_000913/858276‑858574
                                                                                                                                                     |                                                                                                                                                     
gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcag                                                                                                                                                       <  2:523338/150‑1 (MQ=255)
 cccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcagc                                                                                                                                                      >  2:385622/1‑150 (MQ=255)
 cccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcagc                                                                                                                                                      >  2:20712/1‑150 (MQ=255)
 cccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcagc                                                                                                                                                      >  2:465984/1‑150 (MQ=255)
     ggAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGcagccagc                                                                                                                                                  >  1:581482/1‑150 (MQ=255)
         tGTCCGAAACACTTAACATATGGTCCTGTTGTCAACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGCCTCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGAAACGCTGGAAGCAGCCAGAGTTc                                                                                                                                              >  2:1614086/1‑150 (MQ=255)
          gTCCTAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTcc                                                                                                                                             <  2:1175634/150‑1 (MQ=255)
          gTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTcc                                                                                                                                             >  2:1636771/1‑150 (MQ=255)
          gTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTcc                                                                                                                                             <  2:1559153/150‑1 (MQ=255)
            ccGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAg                                                                                                                                           <  1:276501/150‑1 (MQ=255)
            ccGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAg                                                                                                                                           <  1:1331698/150‑1 (MQ=255)
             cGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGt                                                                                                                                          >  2:759434/1‑150 (MQ=255)
             cGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGAGCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatcagc                                                                                                                                          <  2:1598999/150‑7 (MQ=255)
              gAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACTTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatcagcg                                                                                                                                         <  2:1420950/150‑8 (MQ=255)
              gAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatcagcg                                                                                                                                         <  2:226866/150‑8 (MQ=255)
              gAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatcagcg                                                                                                                                         <  1:1152952/150‑8 (MQ=255)
                aaCACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGcc                                                                                                                                       <  2:1044081/150‑1 (MQ=255)
                  cacGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGcccc                                                                                                                                     >  1:1485730/1‑150 (MQ=255)
                   ccGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGAGCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCg                                                                                                                                    <  1:1168514/149‑1 (MQ=255)
                      gAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAg                                                                                                                                 >  1:1469955/1‑150 (MQ=255)
                         cATTTGGTGGTGTTGTCCACAGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCg                                                                                                                              >  1:357448/1‑150 (MQ=255)
                          aTTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCgg                                                                                                                             >  1:994819/1‑150 (MQ=255)
                          aTTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCgg                                                                                                                             >  1:1183950/1‑150 (MQ=255)
                           tttGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGt                                                                                                                            <  1:59591/150‑1 (MQ=255)
                            ttGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTa                                                                                                                           <  1:1138310/150‑1 (MQ=255)
                            ttGGTGGTGTTGTCCACCGTGTACTCTTAAACTTCCGTACGATGTGCCGCATGACAGGGATCATTTTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGATGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCTGAGTCGGTa                                                                                                                           <  1:1261654/150‑1 (MQ=255)
                             tggtggTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTAt                                                                                                                          >  1:1067686/1‑150 (MQ=255)
                              ggtggtGTTGTCGACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATg                                                                                                                         <  1:1137841/150‑1 (MQ=255)
                                tggtgTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATgcg                                                                                                                       >  2:1027619/1‑150 (MQ=255)
                                      tgtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCaa                                                                                                                 >  2:1244306/1‑150 (MQ=255)
                                       gtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAg                                                                                                                >  1:454049/1‑150 (MQ=255)
                                         ccACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAc                                                                                                              >  2:222221/1‑150 (MQ=255)
                                          cACCGTGTACTGTTCAACTTCCGTACGATTTGCCGCATGCCAGGGATCATATTGAGCCAGATTGCCATTAAAGACTTTGGCGAAAACGCGCTGATTCGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAcc                                                                                                             <  1:1737672/150‑1 (MQ=255)
                                          cACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAcc                                                                                                             <  2:207784/150‑1 (MQ=255)
                                            ccGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCtt                                                                                                           >  2:1499304/1‑150 (MQ=255)
                                               tgtACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTgg                                                                                                        <  1:385622/150‑1 (MQ=255)
                                                gtACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGa                                                                                                       <  1:486530/150‑1 (MQ=255)
                                                 tACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAg                                                                                                      <  1:1155692/150‑1 (MQ=255)
                                                     ctTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCg                                                                                                  >  2:1257859/1‑150 (MQ=255)
                                                      ttCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTACGGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGt                                                                                                 <  1:924529/150‑1 (MQ=255)
                                                          aCTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCg                                                                                             >  1:1428672/1‑150 (MQ=255)
                                                          aCTTCCGTACGATGTGCCGCATGCCAGGGATCATAATGCGCCAGATTGCAATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTTCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCg                                                                                             <  2:1379255/150‑1 (MQ=255)
                                                           cTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCTTTCCgg                                                                                            >  1:166479/1‑150 (MQ=255)
                                                            ttCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCggg                                                                                           >  1:1734094/1‑150 (MQ=255)
                                                            ttCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCggg                                                                                           <  1:1551812/150‑1 (MQ=255)
                                                                gTACGATGTGCCGCTTGGGAGGTATCATATTGCGCCATATTGCCATTAAAGTCGTTGGCGAAATCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAGCTTTGGAGGTCGTTCCGGGCGGa                                                                                       <  2:1558399/150‑1 (MQ=255)
                                                                  aCGATGTGCCGCCTGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGAcg                                                                                     <  2:1701947/150‑1 (MQ=255)
                                                                  aCGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGAcg                                                                                     >  2:1212155/1‑150 (MQ=255)
                                                                    gATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGGCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGAcggc                                                                                   <  2:1095829/150‑1 (MQ=255)
                                                                        tgCCGCATGCCAGGGATCATATTGCGCCAGATTGACATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGAcggcatcg                                                                               <  2:907376/150‑5 (MQ=255)
                                                                        ggCCGCACGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGAcggcggcg                                                                               <  2:320507/149‑1 (MQ=255)
                                                                          ccGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGATGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGa                                                                             >  1:945627/1‑150 (MQ=255)
                                                                          ccGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGa                                                                             >  1:799307/1‑150 (MQ=255)
                                                                              aTGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACg                                                                         <  2:1433659/150‑1 (MQ=255)
                                                                                 ccAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTc                                                                      >  2:1553644/1‑150 (MQ=255)
                                                                                  cAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCa                                                                     >  2:392904/1‑150 (MQ=255)
                                                                                    gggATCATATTGCGCCAGATTGCCCTTAAAGACGTTGGCGAACACGCGCTGATACGCTGGAAGCAGCCAGCGGTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAgg                                                                   <  1:234037/150‑1 (MQ=255)
                                                                                     ggATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAggg                                                                  >  2:752572/1‑150 (MQ=255)
                                                                                      gATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGt                                                                 <  1:384907/150‑1 (MQ=255)
                                                                                       aTCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCTGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTa                                                                <  2:933673/150‑1 (MQ=255)
                                                                                          atatTGCGACAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGCAAGAt                                                             >  2:1312742/1‑150 (MQ=255)
                                                                                            atTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGa                                                           >  2:1251704/1‑150 (MQ=255)
                                                                                            atTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGa                                                           >  2:1480581/1‑150 (MQ=255)
                                                                                                cgcCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCt                                                       >  2:422922/1‑150 (MQ=255)
                                                                                                 gcCAGATTGCCATTAAAGACGTTGGCGACAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTc                                                      <  2:241148/150‑1 (MQ=255)
                                                                                                 gcCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCTGACGGCGGCGGATACGGTCAGGGTAGATGACGCTc                                                      >  1:1626453/1‑150 (MQ=255)
                                                                                                  ccAGATTGCCATTAACGACGTGGGCGAAAACGCGCGGATACGCTGGAAGCAGCCAGCGTTCCAGGGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCGca                                                     <  1:1146456/150‑1 (MQ=255)
                                                                                                  ccAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTca                                                     >  1:1769110/1‑150 (MQ=255)
                                                                                                   aaGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCGCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGGCAGGGTAGATGACGCTcac                                                    <  2:865694/149‑1 (MQ=255)
                                                                                                       ttGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGCTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATcc                                                >  1:464046/1‑150 (MQ=255)
                                                                                                          ccATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCggg                                             >  2:609628/1‑150 (MQ=255)
                                                                                                          ccATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCggg                                             <  1:427610/150‑1 (MQ=255)
                                                                                                           cATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGt                                            >  2:1035517/1‑150 (MQ=255)
                                                                                                            aTTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGtt                                           >  1:1222405/1‑150 (MQ=255)
                                                                                                             ttAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTa                                          <  1:759434/150‑1 (MQ=255)
                                                                                                               aaaGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGATaaa                                        >  1:697145/1‑150 (MQ=255)
                                                                                                                aaGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAAc                                       >  2:1380382/1‑150 (MQ=255)
                                                                                                                     gTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCCTTCCAGTGCCCCGGAGTCGGTATGCACTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTg                                  >  2:795418/1‑150 (MQ=255)
                                                                                                                      ttGGCGAAAACGCGCTGATACGCGGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGc                                 <  1:1678478/150‑1 (MQ=255)
                                                                                                                        ggCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCtt                               >  1:123456/1‑150 (MQ=255)
                                                                                                                         gCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCttt                              >  1:823018/1‑150 (MQ=255)
                                                                                                                          cGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTc                             >  1:1114659/1‑150 (MQ=255)
                                                                                                                          cGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTc                             <  2:1183950/150‑1 (MQ=255)
                                                                                                                            aaaaCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCa                           >  1:1761730/1‑150 (MQ=255)
                                                                                                                             aaaCGCGCAGATACGCTGGAAGCAGCCAGCGGTCCAGGGCCCCGGCGGCGGGATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCAt                          <  1:318850/150‑1 (MQ=255)
                                                                                                                               aCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCa                        <  1:1416127/150‑1 (MQ=255)
                                                                                                                                cgcgCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCAc                       <  2:1439246/150‑1 (MQ=255)
                                                                                                                                cgcgCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCAc                       >  2:1518475/1‑150 (MQ=255)
                                                                                                                                  cgcTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACtt                     <  2:594855/150‑1 (MQ=255)
                                                                                                                                      ggaaCGCTGGCTGCTTCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCaa                 >  1:1296014/4‑150 (MQ=255)
                                                                                                                                       aTACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGGTCCGGGCGGCCGGCGGCGGAAACGGTCAGGGTAGATGACGCTCACATCCGGGTGACACCTTTGCTTTCCATCACTTTCaaa                >  1:709491/1‑150 (MQ=255)
                                                                                                                                         aCGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATg              <  1:1212155/150‑1 (MQ=255)
                                                                                                                                           gCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTc            >  1:1709710/1‑150 (MQ=255)
                                                                                                                                            cTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTcc           <  2:1818472/150‑1 (MQ=255)
                                                                                                                                              ggAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCcac         >  1:1712902/1‑150 (MQ=255)
                                                                                                                                              ggAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCAATCACTTTCAAATGTcccc         >  2:221875/1‑148 (MQ=255)
                                                                                                                                                   cagccagcGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACg    >  1:1258615/1‑150 (MQ=255)
                                                                                                                                                   cagccagcGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACg    <  2:419959/150‑1 (MQ=255)
                                                                                                                                                    agccagcGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGa   >  2:778029/1‑150 (MQ=255)
                                                                                                                                                     gccagcGTTCCAGTGCCCCGGAGGCGGTATGCGCTCCAAGACCGTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACAGCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGAt  >  2:1070748/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGAT  >  NC_000913/858276‑858574

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: