Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,822,659 G→A R173H (CGT→CAT)  spoU → putative RNA methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,822,6590GA75.4% 650.9 / 182.1 317R173H (CGT→CAT) spoUputative RNA methylase
Reads supporting (aligned to +/- strand):  ref base G (42/36);  new base A (125/114);  total (167/150)
Fisher's exact test for biased strand distribution p-value = 8.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 319 total aligned reads displayed.

CGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCC  >  NC_000913/3822567‑3822751
                                                                                             |                                                                                            
  gCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATg                                                                                            >  2:713727/1‑93 (MQ=255)
  gCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATg                                                                                            >  2:6714837/1‑93 (MQ=255)
   cATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGa                                                                                           <  1:3499858/93‑1 (MQ=255)
   cATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGa                                                                                           <  2:5689553/93‑1 (MQ=255)
    aTGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGaa                                                                                          <  2:6241985/93‑1 (MQ=255)
     tGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGaaa                                                                                         <  1:4537521/93‑1 (MQ=255)
      ggTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGGTACCTGCATGaaa                                                                                         >  1:2416083/1‑93 (MQ=255)
      ggTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGaaaa                                                                                        <  1:2665823/93‑1 (MQ=255)
           aGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCa                                                                                   <  1:7325286/93‑1 (MQ=255)
            gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCAt                                                                                  <  2:3837843/93‑1 (MQ=255)
              cGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGt                                                                                <  1:1026504/93‑1 (MQ=255)
               gCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGtt                                                                               <  2:5516790/93‑1 (MQ=255)
                cTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTg                                                                              >  1:6899770/1‑93 (MQ=255)
                 tGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGc                                                                             >  2:6006245/1‑93 (MQ=255)
                 tGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGc                                                                             <  2:3012711/93‑1 (MQ=255)
                    aTGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCgg                                                                          >  1:151599/1‑93 (MQ=255)
                      gTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGaa                                                                        <  2:4437169/93‑1 (MQ=255)
                      gTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGaa                                                                        <  1:1591346/93‑1 (MQ=255)
                            gTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCagag                                                                  >  2:1444412/1‑93 (MQ=255)
                             ttGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGc                                                                 >  2:314082/1‑93 (MQ=255)
                              tGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGca                                                                >  1:3127375/1‑93 (MQ=255)
                              tGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGca                                                                >  2:1897831/1‑93 (MQ=255)
                               gCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGcaa                                                               <  2:3495376/93‑1 (MQ=255)
                                ccTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGcaac                                                              <  2:1833352/93‑1 (MQ=255)
                                 cTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaaca                                                             >  1:5082393/1‑93 (MQ=255)
                                 cTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGcaaca                                                             >  1:2291976/1‑93 (MQ=255)
                                 cTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGcaaca                                                             >  1:2004757/1‑93 (MQ=255)
                                  tCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                            >  1:1172289/1‑93 (MQ=255)
                                  tCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                            <  1:4431851/93‑1 (MQ=255)
                                     gCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGc                                                         >  1:1107237/1‑93 (MQ=255)
                                     gCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGc                                                         >  2:4888058/1‑93 (MQ=255)
                                     gCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGc                                                         >  2:5657946/1‑93 (MQ=255)
                                      cccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGcc                                                        <  2:3609050/93‑1 (MQ=255)
                                       ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCt                                                       >  2:584866/1‑93 (MQ=255)
                                        cTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtg                                                      >  2:2786066/1‑93 (MQ=255)
                                        cTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtg                                                      >  1:3032497/1‑93 (MQ=255)
                                        cTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtg                                                      >  1:4207276/1‑93 (MQ=255)
                                        cTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtg                                                      >  1:3331693/1‑93 (MQ=255)
                                         tCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgt                                                     <  2:2534567/93‑1 (MQ=255)
                                         tCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgt                                                     <  1:1317374/93‑1 (MQ=255)
                                          cATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgtt                                                    >  1:3775008/1‑93 (MQ=255)
                                           aTTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttg                                                   >  1:875534/1‑93 (MQ=255)
                                           aTTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttg                                                   >  1:6785340/1‑93 (MQ=255)
                                           aTTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttg                                                   >  1:7136952/1‑93 (MQ=255)
                                            ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttgt                                                  >  1:2577306/1‑93 (MQ=255)
                                            ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttgt                                                  <  2:2257290/93‑1 (MQ=255)
                                             tCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTgttgtt                                                 >  2:2063094/1‑93 (MQ=255)
                                              cTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGttt                                                <  1:4737770/93‑1 (MQ=255)
                                               tttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTg                                               >  1:7353052/1‑93 (MQ=255)
                                               tttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTg                                               >  2:7245109/1‑93 (MQ=255)
                                               tttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTg                                               >  2:5271096/1‑93 (MQ=255)
                                                 tACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGaa                                             <  1:6413542/93‑1 (MQ=255)
                                                 tACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGaa                                             <  1:4471989/93‑1 (MQ=255)
                                                   cGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAgg                                           >  2:1016875/1‑93 (MQ=255)
                                                   cGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAgg                                           >  1:1779664/1‑93 (MQ=255)
                                                    gAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggc                                          >  1:5917927/1‑93 (MQ=255)
                                                    gAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggc                                          >  1:2630090/1‑93 (MQ=255)
                                                     aaGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcg                                         >  1:5650030/1‑93 (MQ=255)
                                                      aGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcgg                                        >  2:4661619/1‑93 (MQ=255)
                                                       gCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcggc                                       <  1:534090/93‑1 (MQ=255)
                                                       gCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcggc                                       >  2:778544/1‑93 (MQ=255)
                                                        cccAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCt                                      >  2:4098870/1‑93 (MQ=255)
                                                        cccAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCt                                      >  1:7243839/1‑93 (MQ=255)
                                                         ccAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTa                                     <  1:7458470/93‑1 (MQ=255)
                                                         ccAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTa                                     >  2:1541969/1‑93 (MQ=255)
                                                          cAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTAt                                    >  2:4535666/1‑93 (MQ=255)
                                                           aGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATc                                   >  1:6011413/1‑93 (MQ=255)
                                                           aGCGTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATc                                   <  1:5477811/93‑1 (MQ=255)
                                                            gCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATcc                                  >  2:3340139/1‑93 (MQ=255)
                                                              gTCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCgg                                >  2:6835962/1‑93 (MQ=255)
                                                               tCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                               <  1:684941/93‑1 (MQ=255)
                                                               tCAGCGGCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                               <  1:2509290/93‑1 (MQ=255)
                                                                    ggCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTgg                          <  2:117399/93‑1 (MQ=255)
                                                                     gCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGc                         >  2:2379848/1‑93 (MQ=255)
                                                                     gCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGc                         >  1:4801749/1‑93 (MQ=255)
                                                                     gCAAAATGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGc                         >  1:4063359/1‑93 (MQ=255)
                                                                       aaaaTGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGa                       >  1:6267639/1‑93 (MQ=255)
                                                                       aaaaTGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGa                       >  1:3742229/1‑93 (MQ=255)
                                                                        aaaTGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGaa                      >  1:906202/1‑93 (MQ=255)
                                                                        aaaTGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGaa                      >  2:2656098/1‑93 (MQ=255)
                                                                         aaTGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGaaa                     <  1:6651746/93‑1 (MQ=255)
                                                                          aTGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAg                    >  2:7064938/1‑93 (MQ=255)
                                                                          aTGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAg                    <  1:4114383/93‑1 (MQ=255)
                                                                           tGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGt                   <  2:5145014/93‑1 (MQ=255)
                                                                           tGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGt                   <  2:5930079/93‑1 (MQ=255)
                                                                           tGCAGGCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGt                   >  1:2847175/1‑93 (MQ=255)
                                                                              aGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGc                >  1:5281560/1‑93 (MQ=255)
                                                                                gCATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaa              <  1:1208794/93‑1 (MQ=255)
                                                                                 cATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaa             <  2:4248030/93‑1 (MQ=255)
                                                                                 cATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaa             >  1:3599960/1‑93 (MQ=255)
                                                                                 cATGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaa             >  2:4568627/1‑93 (MQ=255)
                                                                                   tGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAAc           >  2:3165751/1‑93 (MQ=255)
                                                                                   tGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAAc           <  1:584453/93‑1 (MQ=255)
                                                                                   tGTACCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAAc           >  1:6726874/1‑93 (MQ=255)
                                                                                       aCCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTCAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAACGCa        <  1:6208097/93‑1 (MQ=255)
                                                                                       aCCTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAACGCa        <  2:1361133/93‑1 (MQ=255)
                                                                                         cTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAACGCaaa      <  1:2088533/93‑1 (MQ=255)
                                                                                         cTGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAACGCaaa      >  2:3265016/1‑93 (MQ=255)
                                                                                          tGCATGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAg     >  2:5100170/1‑93 (MQ=255)
                                                                                             gTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGgcc  >  2:164166/1‑93 (MQ=255)
                                                                                             |                                                                                            
CGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCC  >  NC_000913/3822567‑3822751

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