Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,180,699 C→T 5.1% H93H (CAC→CAT) ‡ ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,180,6990CT 5.1% 432.7 / 5.6 134H93H (CAC→CAT) ‡ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base C (61/66);  new base T (2/5);  total (63/71)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

AGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACGCCGGTAA  >  NC_000913/3180551‑3180848
                                                                                                                                                    |                                                                                                                                                     
aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGTCAGCCATGACTGGCGCAca                                                                                                                                                      <  2:1272637/150‑2 (MQ=255)
aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGAAcc                                                                                                                                                      >  1:1040050/1‑150 (MQ=255)
   ttGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGc                                                                                                                                                   <  1:813721/150‑1 (MQ=255)
     gAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGCTGACCAAATTCCGCATTCCAAAACACACATGGAGTTTTGTGCGCCAGTCCTGGCTGACGCACCAGCCa                                                                                                                                                 >  1:232476/1‑150 (MQ=255)
      aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGCTCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCag                                                                                                                                                >  1:152557/1‑149 (MQ=255)
      aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCAt                                                                                                                                                >  1:458060/1‑150 (MQ=255)
      aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCAt                                                                                                                                                >  2:1708303/1‑150 (MQ=255)
       aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAtg                                                                                                                                               <  2:1583995/150‑2 (MQ=255)
       aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATc                                                                                                                                               <  2:1198700/150‑1 (MQ=255)
         aGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAtgac                                                                                                                                             <  2:1566463/150‑4 (MQ=255)
         aGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGc                                                                                                                                             <  1:1695381/150‑1 (MQ=255)
           cTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCtt                                                                                                                                           <  1:342144/150‑1 (MQ=255)
            tGGAAGAAGTCACCGCCGAACTGCACCAGATTTGCCTGAACGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTACGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCttt                                                                                                                                          <  1:941836/150‑1 (MQ=255)
            tGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCttt                                                                                                                                          <  2:1445972/150‑1 (MQ=255)
             ggAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTa                                                                                                                                         >  1:683085/1‑150 (MQ=255)
             ggAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGCAGTTTTGTGCGCCAGTCATGGATTACGCACCAGCCATCGCTTTa                                                                                                                                         >  1:1125780/1‑150 (MQ=255)
               aagaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTAtt                                                                                                                                       <  2:565879/150‑1 (MQ=255)
               aagaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTAtt                                                                                                                                       <  2:1131972/150‑1 (MQ=255)
               aagaagTCACCGCCGAACACCACCCGATGAGCCTTAAAGTGGTGGAAAAAGTGATCGCCAGCGACGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGAGTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTAtt                                                                                                                                       <  1:864579/150‑1 (MQ=255)
                 gaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGCGGAAAAAGTGATCGCAAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcg                                                                                                                                     <  2:1600789/150‑1 (MQ=255)
                 gaagGCACCGCCGAAATGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTTATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcg                                                                                                                                     <  2:833072/150‑1 (MQ=255)
                  aagGCACCGCCGAACTGCACCTGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAACAACTCCTGGAGTTATGTGCGCCAGTCATGGCTGTCTCACCAGACATCGCTTTATTcgc                                                                                                                                    >  2:887440/1‑150 (MQ=255)
                   agTCACCGCCGAACTGCACGAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcgcg                                                                                                                                   <  2:1340987/150‑1 (MQ=255)
                   agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcgcg                                                                                                                                   <  1:793040/150‑1 (MQ=255)
                     tCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTc                                                                                                                                 <  1:742993/150‑1 (MQ=255)
                      cACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCt                                                                                                                                <  2:402851/150‑1 (MQ=255)
                        ccgccgAACTGCACCATATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAACCATCGCTTTATTCGCGTCTTg                                                                                                                              <  2:596239/150‑1 (MQ=255)
                        ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                              <  2:31202/150‑1 (MQ=255)
                        ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                              <  2:411449/150‑1 (MQ=255)
                        ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGGTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                              <  2:325101/150‑1 (MQ=255)
                         cgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGTTCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGa                                                                                                                             >  2:501654/1‑150 (MQ=255)
                         cgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGa                                                                                                                             >  2:1651773/1‑150 (MQ=255)
                          gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGACTGACGCACCAGCCATCGCtt                                                                                                                                           <  1:1698737/135‑1 (MQ=255)
                           ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                           <  1:679400/150‑1 (MQ=255)
                           ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                           >  2:1656406/1‑150 (MQ=255)
                           ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                           >  1:1010320/1‑150 (MQ=255)
                           ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                           <  2:260690/150‑1 (MQ=255)
                                cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                                      <  2:976714/150‑1 (MQ=255)
                                cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                                      <  1:179522/150‑1 (MQ=255)
                                 tGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGt                                                                                                                     >  2:421062/1‑150 (MQ=255)
                                 tGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGt                                                                                                                     >  1:837131/1‑150 (MQ=255)
                                   cACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTTATCTGGCGTgg                                                                                                                   <  1:981617/150‑1 (MQ=255)
                                    aCCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTggg                                                                                                                  <  1:933931/150‑1 (MQ=255)
                                    aCCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTggg                                                                                                                  >  1:1311907/1‑150 (MQ=255)
                                       aGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGCCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATg                                                                                                               >  2:892061/1‑150 (MQ=255)
                                        gATGTGTCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATgg                                                                                                              >  1:47726/1‑150 (MQ=255)
                                        gATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATgg                                                                                                              <  1:224853/150‑1 (MQ=255)
                                             gCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTg                                                                                                         <  1:272265/150‑1 (MQ=255)
                                               cTGATAGTGGTGGAAAATGTGATCGCCAGCGATGAGCTGATGACCAAATTGAGCATTCCAAAACACACCTGGAGGTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATGGATTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                                                                       <  1:264742/150‑1 (MQ=255)
                                               cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCTAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                                                                       >  2:897896/1‑150 (MQ=255)
                                               cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                                                                       >  2:987270/1‑150 (MQ=255)
                                                 gAAAGTGGTGGAAAAAGTGAACGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGa                                                                                                     <  2:1190738/150‑1 (MQ=255)
                                                    aGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAAcc                                                                                                  <  1:1656406/150‑1 (MQ=255)
                                                         tggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGAGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAAc                                                                                             >  1:194634/1‑150 (MQ=255)
                                                          ggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACt                                                                                            <  2:654383/150‑1 (MQ=255)
                                                            aaaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTc                                                                                          <  2:906081/150‑1 (MQ=255)
                                                             aaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTTATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCt                                                                                         <  2:850586/150‑1 (MQ=255)
                                                             aaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCt                                                                                         <  2:806871/150‑1 (MQ=255)
                                                               aaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTgg                                                                                       <  2:1671284/150‑1 (MQ=255)
                                                                 gTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAACCTTCTGGaa                                                                                     >  1:368947/1‑150 (MQ=255)
                                                                  tGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaaa                                                                                    >  1:923605/1‑150 (MQ=255)
                                                                   gATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaaaa                                                                                   <  2:773569/150‑1 (MQ=255)
                                                                     tCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATa                                                                                 <  1:1637808/150‑1 (MQ=255)
                                                                        ccAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGGTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACg                                                                              <  1:1155798/150‑1 (MQ=255)
                                                                          aGCGATGAGCTTATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTTTGCGCCATTCATGGCTGACTCACCAGCCATCGCTTTATTCTCGTCTTGATCTGGCGTGGGATGTCACTGGTGAACCTAAACTTCTGGAAAATAACGcc                                                                            <  1:1613359/150‑1 (MQ=255)
                                                                              aTGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGGGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATa                                                                        >  1:143713/1‑150 (MQ=255)
                                                                                  gCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGcc                                                                    <  2:187030/150‑1 (MQ=255)
                                                                                        gaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTc                                                              >  1:1599092/1‑150 (MQ=255)
                                                                                          ccAAATTCCGCATTCCAACACACAGCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGGGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                            <  1:604626/150‑1 (MQ=255)
                                                                                          ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                            >  2:1469220/1‑150 (MQ=255)
                                                                                            aaaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACta                                                          >  1:629368/1‑150 (MQ=255)
                                                                                             aaTTCCGGATTCCAAAACACACCTGGAGTTTTGTGCGACAGTCAGGTCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGTCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtat                                                         <  1:981263/150‑1 (MQ=255)
                                                                                             aaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtat                                                         >  2:1442594/1‑150 (MQ=255)
                                                                                              aTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtata                                                        <  1:743090/150‑1 (MQ=255)
                                                                                                tCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACg                                                      <  1:531141/150‑1 (MQ=255)
                                                                                                 ccGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGa                                                     >  1:863612/1‑150 (MQ=255)
                                                                                                    cATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAggc                                                  >  1:510650/1‑150 (MQ=255)
                                                                                                     aTTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAggcg                                                 >  1:1267664/1‑150 (MQ=255)
                                                                                                         cAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGt                                             <  1:1603528/150‑1 (MQ=255)
                                                                                                          aaaaCACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGtt                                            >  2:1553664/1‑150 (MQ=255)
                                                                                                               acacCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCAACGCTTTATTATCTTCTTTATCTGGCGTGGGATGTCAATGTTGAACATAAACTTCTGGAAAATAACGCCGATACGCCAACGTAAATATACGAGGCGGAGTTCTTTc                                       >  2:1111775/1‑150 (MQ=255)
                                                                                                               acacCTGGAGTTGTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTc                                       <  2:908810/150‑1 (MQ=255)
                                                                                                                  ccTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGGCGTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCGTTCAGt                                    >  2:1040383/1‑150 (MQ=255)
                                                                                                                        gTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATct                              >  1:1727090/1‑150 (MQ=255)
                                                                                                                        gTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATct                              >  1:1647728/1‑150 (MQ=255)
                                                                                                                           tggtgCGTCAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggc                           >  2:413913/3‑150 (MQ=255)
                                                                                                                            tgtgCGCCAGTCAGGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggct                          <  2:193374/150‑1 (MQ=255)
                                                                                                                            ggtgCGTCAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggct                          >  1:1279397/2‑150 (MQ=255)
                                                                                                                             gtgCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggctg                         >  2:742032/1‑150 (MQ=255)
                                                                                                                               gcgcCAGGCATGGCTGACGTACCAGCCATCGCTTTATTCGCGTCTTGATATGGCGTGGGATGGCACTGTTGATCCTAAACTTCTGGAATATATCGCGGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggctggt                       <  2:232476/150‑2 (MQ=255)
                                                                                                                                 taaagcgaTGGCTGGTGCGTCAGCCATCGCTTTTTTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAg                     >  1:113565/8‑150 (MQ=255)
                                                                                                                                  tcAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGa                    >  1:787973/2‑150 (MQ=255)
                                                                                                                                  ccAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGa                    >  1:616668/1‑150 (MQ=255)
                                                                                                                                      cgaTGGCTGGTGCGTCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAg                >  1:757078/3‑150 (MQ=255)
                                                                                                                                       cATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGc               <  1:1614785/150‑1 (MQ=255)
                                                                                                                                        aTGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCt              >  1:845087/1‑150 (MQ=255)
                                                                                                                                          ggCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTa            >  1:399884/1‑150 (MQ=255)
                                                                                                                                                cGCACCAGCAATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACGCCg      <  2:1279397/150‑1 (MQ=255)
                                                                                                                                                 gCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGACCCTAAACTTCTGGAAAATAACGCCGATATGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACGCCgg     >  1:1232347/1‑150 (MQ=255)
                                                                                                                                                    ccAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACGCCGGTaa  >  1:490460/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACGCCGGTAA  >  NC_000913/3180551‑3180848

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: