Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,143,747 A→C 5.4% I207M (ATT→ATG rne ← ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,143,7470AC 5.4% 296.4 / 7.4 93I207M (ATT→ATGrneribonuclease E
Reads supporting (aligned to +/- strand):  ref base A (52/36);  new base C (3/2);  total (55/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCTTCCGGCAGTTCA  >  NC_000913/1143601‑1143894
                                                                                                                                                  |                                                                                                                                                   
aTTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATc                                                                                                                                                  <  2:189731/150‑1 (MQ=255)
  tttGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAg                                                                                                                                                <  1:908958/150‑1 (MQ=255)
  tttGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAg                                                                                                                                                <  2:586607/150‑1 (MQ=255)
       tgctgAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATTCATCAGGagag                                                                                                                                           <  2:1447442/150‑4 (MQ=255)
       tgctgAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACg                                                                                                                                           >  1:1276357/1‑150 (MQ=255)
        gctgAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATTCATCAGGagagc                                                                                                                                          <  2:627360/150‑5 (MQ=255)
           gAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGc                                                                                                                                       >  2:1444017/1‑150 (MQ=255)
             aaTCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCg                                                                                                                                     <  1:543955/150‑1 (MQ=255)
             aaTCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCg                                                                                                                                     >  2:1197566/1‑150 (MQ=255)
              aTCCGGGCGCCCTAAAGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGAGGAATCAGGAACGGGGCCgg                                                                                                                                    >  2:885748/1‑150 (MQ=255)
                ccGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGggc                                                                                                                                  >  2:1610409/1‑150 (MQ=255)
                 cGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGggcg                                                                                                                                 <  1:1444017/150‑1 (MQ=255)
                 cGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGggcg                                                                                                                                 <  2:1990949/150‑1 (MQ=255)
                  gggCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGggcgg                                                                                                                                >  2:1242231/1‑150 (MQ=255)
                    gCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCt                                                                                                                              <  2:1265254/150‑1 (MQ=255)
                     cGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCtt                                                                                                                             <  2:271842/150‑1 (MQ=255)
                     cGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCtt                                                                                                                             <  2:1859076/150‑1 (MQ=255)
                       aCCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTc                                                                                                                           >  2:589453/1‑150 (MQ=255)
                          tAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGc                                                                                                                        >  2:1592462/1‑150 (MQ=255)
                           aaTGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCg                                                                                                                       <  1:1373895/150‑1 (MQ=255)
                           aaTGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCg                                                                                                                       >  1:1310344/1‑150 (MQ=255)
                              gCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCt                                                                                                                    >  2:41825/1‑150 (MQ=255)
                                    atatGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTg                                                                                                              <  1:549553/150‑1 (MQ=255)
                                     tatGCTGACGTGCCAGTTCGAGCTCTTTCTTGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCTGAATGAGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGa                                                                                                             <  2:1351425/150‑1 (MQ=255)
                                       tGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCtt                                                                                                                             >  2:793864/1‑132 (MQ=255)
                                       tGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCtt                                                                                                                             <  1:793864/132‑1 (MQ=255)
                                       tGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATg                                                                                                           <  2:1348378/150‑1 (MQ=255)
                                       tGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATg                                                                                                           <  1:1490522/150‑1 (MQ=255)
                                        gCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATgg                                                                                                          >  2:7681/1‑150 (MQ=255)
                                                gCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATc                                                                                                                                                  <  2:430389/102‑1 (MQ=255)
                                                gCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATc                                                                                                                                                  >  1:430389/1‑102 (MQ=255)
                                                gCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAg                                                                                                  >  2:1378809/1‑150 (MQ=255)
                                                gCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAg                                                                                                  >  2:99/1‑150 (MQ=255)
                                                gCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAg                                                                                                  >  1:513874/1‑150 (MQ=255)
                                                  cAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAgtg                                                                                                >  2:414544/1‑150 (MQ=255)
                                                  cAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAgtg                                                                                                >  2:1163542/1‑150 (MQ=255)
                                                   aGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAgtgt                                                                                               >  2:1670587/1‑150 (MQ=255)
                                                       cGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCa                                                                                           <  2:1992953/150‑1 (MQ=255)
                                                        gAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAg                                                                                          <  1:896551/150‑1 (MQ=255)
                                                          gCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGAc                                                                                        >  1:562082/1‑150 (MQ=255)
                                                                  gggTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCtt                                                                                >  1:908671/1‑150 (MQ=255)
                                                                   ggTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTa                                                                               <  1:1592462/150‑1 (MQ=255)
                                                                    gTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTaa                                                                              >  1:682698/1‑150 (MQ=255)
                                                                           aTAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCAt                                                                       >  2:812511/1‑150 (MQ=255)
                                                                            tAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCAtt                                                                      >  1:133810/1‑150 (MQ=255)
                                                                                 aTTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAgc                                                                 >  2:626773/1‑150 (MQ=255)
                                                                                    tCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGcc                                                              >  2:414532/1‑150 (MQ=255)
                                                                                     cGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTATCAGACGGAAGCTTAAATCCCATTGCAGAGCCt                                                             >  1:1872902/1‑150 (MQ=255)
                                                                                       ccGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTca                                                           >  1:1994518/1‑150 (MQ=255)
                                                                                             tCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGAt                                                     >  1:1041702/1‑150 (MQ=255)
                                                                                             tCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGAt                                                     <  2:1994518/150‑1 (MQ=255)
                                                                                             tCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGAt                                                     >  1:275336/1‑150 (MQ=255)
                                                                                             tCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGAt                                                     >  1:1048507/1‑150 (MQ=255)
                                                                                              ccTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAAt                                                                            <  2:1363800/126‑1 (MQ=255)
                                                                                              ccTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAAt                                                                            >  1:1363800/1‑126 (MQ=255)
                                                                                              ccTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGcc                                                              >  2:1741/1‑140 (MQ=255)
                                                                                              ccTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGcc                                                              <  1:1741/140‑1 (MQ=255)
                                                                                              ccTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGAtt                                                    >  2:1702097/1‑150 (MQ=255)
                                                                                               cTGACGTAAGTAATAGCGGAATGAGCGAACTATCAAGTTGCTCTACTGATGAATAAGGAACGGGGCCGGGAGGCTTTCAGCGGCTTTTTTGATGGCTTCAAAGTGTTTCAGAAGGAAGCTTAAATAAAATTGCAGCGCCTCAGCAGAttt                                                   >  2:1742176/1‑150 (MQ=255)
                                                                                                 gACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGc                                                 <  1:361475/150‑1 (MQ=255)
                                                                                                 gACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGc                                                 >  2:217302/1‑150 (MQ=255)
                                                                                                  aCGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGcc                                                <  2:1510455/150‑1 (MQ=255)
                                                                                                   cGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTATATCGCATTGCAGCGCCTCAGCAGATTTGcgg                                               >  2:365802/1‑148 (MQ=255)
                                                                                                   cGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCg                                               <  2:1750587/150‑1 (MQ=255)
                                                                                                    gtaagtaaTCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGa                                              >  2:388400/1‑150 (MQ=255)
                                                                                                      aagtaaTCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACg                                            <  1:142751/150‑1 (MQ=255)
                                                                                                        gtaaTCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGcc                                          >  1:581329/1‑150 (MQ=255)
                                                                                                         taaTCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCa                                         >  1:411371/1‑150 (MQ=255)
                                                                                                          aaTCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAg                                        >  1:1898414/1‑150 (MQ=255)
                                                                                                           aTCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGc                                       >  2:53693/1‑150 (MQ=255)
                                                                                                           aTCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGc                                       >  1:233538/1‑150 (MQ=255)
                                                                                                                 gAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTgcg                                 <  2:520304/150‑1 (MQ=255)
                                                                                                                  aaTGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTgcgc                                <  1:1610409/150‑1 (MQ=255)
                                                                                                                  aaTGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTgcgc                                <  2:223726/150‑1 (MQ=255)
                                                                                                                   aTGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCa                               <  2:1815533/150‑1 (MQ=255)
                                                                                                                   aTGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCa                               <  1:1197566/150‑1 (MQ=255)
                                                                                                                    tGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCAc                              >  1:975182/1‑150 (MQ=255)
                                                                                                                     gcgcGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGCCGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACg                             >  1:1052904/1‑150 (MQ=255)
                                                                                                                     gcgcGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCTCACg                             <  2:562082/150‑1 (MQ=255)
                                                                                                                      cgcgAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGa                            >  2:897579/1‑150 (MQ=255)
                                                                                                                       gcgAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCCGCGGTGCGCACGAt                           >  2:537394/1‑150 (MQ=255)
                                                                                                                        cgAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAGATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATa                          >  1:1934894/1‑150 (MQ=255)
                                                                                                                        cgAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATa                          <  1:1163542/150‑1 (MQ=255)
                                                                                                                            cGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGcc                      <  2:275336/150‑1 (MQ=255)
                                                                                                                               ttACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCa                   >  2:1441188/1‑150 (MQ=255)
                                                                                                                                  cGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTTAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGc                >  2:1973626/1‑150 (MQ=255)
                                                                                                                                  cGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCACCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGc                >  2:1802542/1‑150 (MQ=255)
                                                                                                                                   gTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGcc               <  1:388400/150‑1 (MQ=255)
                                                                                                                                     tGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATAGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCtt             >  1:859818/1‑150 (MQ=255)
                                                                                                                                     cccgttcCTGATTCATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCtt             >  1:751203/6‑150 (MQ=255)
                                                                                                                                      ccgttcCTGATTCATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCTTc            >  1:445369/5‑150 (MQ=255)
                                                                                                                                      ccgttcCTGATTCATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCTTc            >  1:565920/5‑150 (MQ=255)
                                                                                                                                          tcCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCTTCCGGc        <  1:1378809/150‑1 (MQ=255)
                                                                                                                                              gatgaATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCTTCCGGCAGtt    <  2:1898414/150‑1 (MQ=255)
                                                                                                                                                tgaATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCTTCCGGCAGTTCa  >  1:1290240/1‑150 (MQ=255)
                                                                                                                                                tgaATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCTTCCGGCAGTTCa  >  2:1398096/1‑150 (MQ=255)
                                                                                                                                                tgaATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCTTCCGGCAGTTCa  <  2:233538/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTGATGGCTTCCCAGTGTTTCAGACGGAAGCTTAAATCCCATTGCAGCGCCTCAGCAGATTTGCCGACGCCAGCGGTGCGCACGATAAGCCCCATGCCTTCCGGCAGTTCA  >  NC_000913/1143601‑1143894

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