Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,822,652 T→C Y171H (TAC→CAC)  spoU → putative RNA methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,822,6520TC100.0% 944.0 / NA 282Y171H (TAC→CAC) spoUputative RNA methylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (145/137);  total (145/137)

Only 100 of 282 total aligned reads displayed.

TGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAA  >  NC_000913/3822563‑3822741
                                                                                         |                                                                                         
tgatCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCGTGCAcc                                                                                        <  1:5558206/93‑1 (MQ=255)
tgatCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCAcc                                                                                        <  2:1545187/93‑1 (MQ=255)
tgatCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCAcc                                                                                        <  2:6023684/93‑1 (MQ=255)
 gatCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCt                                                                                       <  2:2717075/93‑1 (MQ=255)
  atCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTg                                                                                      >  1:389307/1‑93 (MQ=255)
   tCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGc                                                                                     <  1:5617961/93‑1 (MQ=255)
     ggCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGt                                                                                   >  1:1823005/1‑93 (MQ=255)
      gCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTg                                                                                  >  2:3185035/1‑93 (MQ=255)
      gCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTg                                                                                  >  1:4326027/1‑93 (MQ=255)
      gCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTg                                                                                  >  2:6657538/1‑93 (MQ=255)
       cATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGa                                                                                 >  1:399599/1‑93 (MQ=255)
          ggTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGaaaa                                                                              <  2:1978918/93‑1 (MQ=255)
              cAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGc                                                                          <  2:1262878/93‑1 (MQ=255)
              cAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGc                                                                          <  2:5247215/93‑1 (MQ=255)
              cAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGc                                                                          <  1:3232004/93‑1 (MQ=255)
               aGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCa                                                                         <  1:1616303/93‑1 (MQ=255)
                gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCAt                                                                        >  1:4243104/1‑93 (MQ=255)
                gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCAt                                                                        >  2:4860421/1‑93 (MQ=255)
                   gCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGtt                                                                     <  1:3801286/93‑1 (MQ=255)
                     tGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGttac                                                                   <  1:3040193/93‑3 (MQ=255)
                     tGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGc                                                                   >  2:257782/1‑93 (MQ=255)
                      gAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGcc                                                                  <  2:2276849/93‑1 (MQ=255)
                       aaTGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCg                                                                 >  2:3077369/1‑93 (MQ=255)
                         tGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGa                                                               >  2:5095007/1‑93 (MQ=255)
                         tGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGa                                                               >  2:1865315/1‑93 (MQ=255)
                             tCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCa                                                           <  1:3726699/93‑1 (MQ=255)
                              ccGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCag                                                          <  2:590195/93‑1 (MQ=255)
                               cGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCaga                                                         <  1:1278102/93‑1 (MQ=255)
                                gTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCagag                                                        <  1:6289485/93‑1 (MQ=255)
                                 ttGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGc                                                       >  2:3817696/1‑93 (MQ=255)
                                    ccTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaac                                                    <  2:713852/93‑1 (MQ=255)
                                     cTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaaca                                                   >  2:8912/1‑93 (MQ=255)
                                      tCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                  <  1:5711128/93‑1 (MQ=255)
                                      tCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                  >  1:4193339/1‑93 (MQ=255)
                                       cAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCaacaac                                                 <  1:5087120/93‑1 (MQ=255)
                                       cAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCaacaac                                                 <  2:5612959/93‑1 (MQ=255)
                                        aGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACg                                                >  2:6085120/1‑93 (MQ=255)
                                        aGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACg                                                >  1:6440231/1‑93 (MQ=255)
                                         gCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGc                                               >  1:1003733/1‑93 (MQ=255)
                                         gCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGc                                               >  2:63867/1‑93 (MQ=255)
                                         gCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGc                                               >  1:4621322/1‑93 (MQ=255)
                                          cccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGcc                                              <  1:4813439/93‑1 (MQ=255)
                                          cccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGcc                                              >  1:6723979/1‑93 (MQ=255)
                                          cccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGcc                                              >  2:1699389/1‑93 (MQ=255)
                                           ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCt                                             >  2:5551848/1‑93 (MQ=255)
                                           ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCt                                             >  1:2676621/1‑93 (MQ=255)
                                           ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCt                                             <  2:6513593/93‑1 (MQ=255)
                                           ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCt                                             <  2:3405142/93‑1 (MQ=255)
                                            cTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtg                                            >  1:4980584/1‑93 (MQ=255)
                                            cTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtg                                            >  2:3299151/1‑93 (MQ=255)
                                             tCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgt                                           <  2:1192164/93‑1 (MQ=255)
                                             tCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgt                                           >  1:2138445/1‑93 (MQ=255)
                                              cATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgtt                                          >  2:6230959/1‑93 (MQ=255)
                                               aTTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttg                                         <  2:2744823/93‑1 (MQ=255)
                                               aTTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttg                                         >  1:831198/1‑93 (MQ=255)
                                               aTTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttg                                         >  1:4734095/1‑93 (MQ=255)
                                                ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttgt                                        >  1:58880/1‑93 (MQ=255)
                                                 tCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTgttgtt                                       <  1:1207235/93‑1 (MQ=255)
                                                   tttACGAAGCCCAGCGTCAGCTGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTg                                     <  2:5286203/93‑1 (MQ=255)
                                                     tACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGaa                                   <  1:3409233/93‑1 (MQ=255)
                                                      aCGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAg                                  <  2:220614/93‑1 (MQ=255)
                                                        gAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggc                                <  1:6693590/93‑1 (MQ=255)
                                                         aaGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcg                               >  2:5420458/1‑93 (MQ=255)
                                                         aaGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcg                               >  1:6024140/1‑93 (MQ=255)
                                                          aGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcgg                              <  2:2890036/93‑1 (MQ=255)
                                                           gCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcggc                             <  1:972118/93‑1 (MQ=255)
                                                            cccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCt                            <  2:3532834/93‑1 (MQ=255)
                                                             ccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTa                           >  2:1369364/1‑93 (MQ=255)
                                                               aGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATc                         >  2:2152566/1‑93 (MQ=255)
                                                                 cGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCg                       <  2:4436477/93‑1 (MQ=255)
                                                                  gTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCgg                      >  2:4642975/1‑93 (MQ=255)
                                                                  gTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCgg                      >  2:3651476/1‑93 (MQ=255)
                                                                   tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCggg                     >  1:582574/1‑92 (MQ=255)
                                                                   tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                     >  1:179214/1‑93 (MQ=255)
                                                                   tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                     >  2:1836998/1‑93 (MQ=255)
                                                                   tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                     >  2:2237007/1‑93 (MQ=255)
                                                                   tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                     <  1:5388436/93‑1 (MQ=255)
                                                                     aGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGc                   <  1:529460/93‑1 (MQ=255)
                                                                       cGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTg                 <  1:4114020/93‑1 (MQ=255)
                                                                        ggCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAATCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTgg                <  2:4325981/93‑1 (MQ=255)
                                                                        ggCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTgg                >  2:499564/1‑93 (MQ=255)
                                                                         gCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGc               <  1:137680/93‑1 (MQ=255)
                                                                         gCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGc               <  1:487842/93‑1 (MQ=255)
                                                                         gCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGc               <  1:2359/93‑1 (MQ=255)
                                                                         gCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGc               <  1:2481404/93‑1 (MQ=255)
                                                                           aaaaTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGa             >  2:1078076/1‑93 (MQ=255)
                                                                           aaaaTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGa             >  1:1994395/1‑93 (MQ=255)
                                                                           aaaaTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGGTGTTTGAAGGCGGCTATCCGGTGGTGGCGa             >  1:2830504/1‑93 (MQ=255)
                                                                            aaaTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGaa            <  2:4004953/93‑1 (MQ=255)
                                                                              aTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAg          >  2:2440097/1‑93 (MQ=255)
                                                                                gCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTc        <  1:6502538/93‑1 (MQ=255)
                                                                                gCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTc        <  2:4720474/93‑1 (MQ=255)
                                                                                 cAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCg       >  2:1841395/1‑93 (MQ=255)
                                                                                 cAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCg       <  1:5921937/93‑1 (MQ=255)
                                                                                  aGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGc      <  1:796174/93‑1 (MQ=255)
                                                                                  aGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGc      <  2:5736082/93‑1 (MQ=255)
                                                                                  aGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGc      <  2:4773781/93‑1 (MQ=255)
                                                                                   ggCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCa     >  1:3323966/1‑93 (MQ=255)
                                                                                      aTGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaaa  >  1:343392/1‑93 (MQ=255)
                                                                                      aTGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaaa  >  2:820565/1‑93 (MQ=255)
                                                                                         |                                                                                         
TGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAA  >  NC_000913/3822563‑3822741

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: