Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,652 G→C 6.6% A578A (GCG→GCC rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,6520GC 6.6% 221.5 / 3.4 81A578A (GCG→GCCrsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base G (21/54);  new base C (3/3);  total (24/57)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCG  >  NC_000913/1708506‑1708788
                                                                                                                                                  |                                                                                                                                        
gCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGaa                                                                                                                                       >  1:213781/1‑150 (MQ=255)
 cccGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                                                                      <  1:655644/150‑1 (MQ=255)
 cccGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                                                                      <  1:1729914/150‑1 (MQ=255)
  ccGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAcc                                                                                                                                     >  2:2009877/1‑150 (MQ=255)
    gTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTACTGCGGAACCAg                                                                                                                                   >  1:1989452/1‑150 (MQ=255)
     tAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGa                                                                                                                                  >  1:1202915/1‑150 (MQ=255)
        aGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGaaca                                                                                                                               <  2:825715/150‑1 (MQ=255)
           gcAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAacaaca                                                                                                                            <  2:1565234/150‑1 (MQ=255)
            cgcgCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAg                                                                                                                           <  2:917551/147‑1 (MQ=255)
            cAAGCCAGAGCGAATCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGACGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAAGAGGCTAATGCGGAACCAGAACAACAg                                                                                                                           <  1:175103/150‑1 (MQ=255)
               gCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAACCTGGAACCGCAACAGGCTAATGCGGAACCAGAACAACAGGTc                                                                                                                        >  1:1676623/1‑150 (MQ=255)
               gCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGAGCCAAAGCGCGCAAGCTGGCACAGCAACAGGCTAATACGGAACCAGAAAAACAGGTc                                                                                                                        >  2:1863947/1‑150 (MQ=255)
                 cAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAAAGCCGAACCAGAACACCAGGTCgc                                                                                                                      >  1:844875/1‑149 (MQ=255)
                 cAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGa                                                                                                                      <  1:1112621/150‑1 (MQ=255)
                   gagCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATc                                                                                                                    >  2:1683550/1‑150 (MQ=255)
                    agCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATcc                                                                                                                   >  2:1018487/1‑150 (MQ=255)
                       gAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCcgcg                                                                                                                <  2:396932/150‑1 (MQ=255)
                         aaCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAACAACAGGTCGATCCGCGCa                                                                                                              >  1:1949082/1‑150 (MQ=255)
                           cAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCaaa                                                                                                            <  2:1536047/150‑1 (MQ=255)
                            aGGCAGATCTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAg                                                                                                           <  1:826049/150‑1 (MQ=255)
                            aGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAg                                                                                                           <  2:1710947/150‑1 (MQ=255)
                             ggCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAgc                                                                                                          >  1:172937/1‑150 (MQ=255)
                               cAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGACGAGCCGCGCAAAgccc                                                                                                        >  1:1629227/1‑149 (MQ=255)
                               aaGATCTGCAGGAAAATCCGGACGCAGCAAACGTTGCTGATCCTCGTAAAACTGCCGTTGAAGCAGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccg                                                                                                        <  2:172937/149‑1 (MQ=17)
                                 caaCTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAgccgcc                                                                                                      <  2:1287610/149‑1 (MQ=255)
                                   aCTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTAGAAGCAGCTATCGCCCGCGCCAACGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGccgccgg                                                                                                    >  2:1495427/1‑149 (MQ=255)
                                    cTGCAGCAAACTAACGACGCAGCAACCGTTTCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGTAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                                   <  1:1564008/150‑1 (MQ=255)
                                    cTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                                   <  1:1067016/150‑1 (MQ=255)
                                    cTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGACGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                                   <  1:91767/150‑1 (MQ=255)
                                     tGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCg                                                                                                  <  2:1111382/150‑1 (MQ=255)
                                     tGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTAGATCCGCGCAAAGCAGCCGtac                                                                                                  >  2:301514/1‑148 (MQ=255)
                                      gcagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCAAGAACAACAGCTCGATCCGCGCACAGACGCCGTCgc                                                                                                 >  1:190412/1‑149 (MQ=25)
                                       cagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGCAGCAGCCACCTCGTAAAGCCCGACCAGAACAGCAGGTCGATCCGCGCAACTCCGCCGTAGaa                                                                                                >  1:450066/1‑150 (MQ=255)
                                       cagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGATCCGGTCGAACCGCGCCAGGCCGCCGTGGac                                                                                                >  1:467469/1‑149 (MQ=21)
                                        tgcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                                                                                               <  2:452587/149‑1 (MQ=255)
                                           aaaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCg                                                                                            <  2:391675/150‑1 (MQ=255)
                                           aaaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTAGAAGCAGCTATAGCCCGCGCCAAAGCGCGCAAGATGGAACACCAACAGGCTAATGAGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCg                                                                                            >  2:145782/1‑150 (MQ=21)
                                            aaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                                           <  1:1018487/150‑1 (MQ=255)
                                            aaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                                           <  1:1600986/150‑1 (MQ=255)
                                            aaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                                           <  1:1815806/150‑1 (MQ=255)
                                            aaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                                           <  1:184057/150‑1 (MQ=255)
                                             aCTAACGACGCATCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGAAACAGGCTAATTAGGAACCACAACAACACTTCTATCCGCGCAAAGCCGCCGTCGAAGCCGcg                                                                                          >  2:425524/1‑149 (MQ=21)
                                              cTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAAGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                                                         >  1:676746/1‑150 (MQ=255)
                                                 ccgacgCAGCCCCCCTTGTTGATCCTCGCAAACCAGCCGGCGAAGCAGCTATAGCCCGCGCCACAGCGCGCAAGCGGGAACAGCAACAGGCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg                                                                                      <  1:145782/149‑1 (MQ=11)
                                                 acgacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg                                                                                      <  2:1756503/150‑1 (MQ=255)
                                                   gacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                                                    <  2:426255/150‑1 (MQ=255)
                                                    acgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGTTCGATCAGCGCATAGCCGCCGTCGAAGCCGCTAtagacc                                                                                   >  1:641008/1‑145 (MQ=18)
                                                    acgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGccc                                                                                   <  1:1565111/150‑1 (MQ=255)
                                                    acgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGGTGAAGCAGCTATCTCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGccc                                                                                   <  2:169852/150‑1 (MQ=255)
                                                       cagcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg                                                                                <  2:1180997/150‑1 (MQ=255)
                                                         gcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                              <  2:718867/150‑1 (MQ=255)
                                                         gcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                              <  2:343755/150‑1 (MQ=255)
                                                         gcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                              <  1:1569696/150‑1 (MQ=255)
                                                         gcaACCGAGGATTATCCACGAAAAACTGCCGGAGACGCAGCTATCTCCCGCGACAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                              <  1:262865/150‑1 (MQ=14)
                                                             ccGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAg                                                                          <  2:526369/150‑1 (MQ=255)
                                                                 tGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                      <  2:975531/150‑1 (MQ=255)
                                                                 tGCTGATCCACGAAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                      <  2:1825170/150‑1 (MQ=255)
                                                                 ggCTGATCCCCGTAAAACGGCCGTTGTAGCAGCTATCGCCCGCGCCAAAGCGAGCAAGATTGAACAGCAACAGGCTAATGCGTAACCAGAACAACATGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                      <  2:338954/149‑1 (MQ=17)
                                                                   cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCa                                                                    <  1:404530/150‑1 (MQ=255)
                                                                         cACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAATCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTgg                                                              <  1:856128/150‑1 (MQ=255)
                                                                          aCGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                             <  1:2009877/150‑1 (MQ=255)
                                                                               aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                                        <  1:578370/150‑1 (MQ=255)
                                                                               aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                                        <  2:745636/150‑1 (MQ=255)
                                                                                  cTGCCGTTGTAGCAGCTATAGCCCGCGCCAAATCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAGCAGAACAACAGGTCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAAc                                                     <  2:35338/150‑1 (MQ=14)
                                                                                   ttccGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACa                                                    <  1:1675221/148‑1 (MQ=255)
                                                                                   tGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACa                                                    <  1:1713121/150‑1 (MQ=255)
                                                                                          gAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg                                                                                                                <  2:421967/83‑1 (MQ=255)
                                                                                          gAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg                                                                                                                >  1:421967/1‑83 (MQ=255)
                                                                                            agcagcTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGc                                           <  1:896406/150‑1 (MQ=255)
                                                                                              cagcTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTAGATCCGCGCAAAGCAGCCGTCGAAAACGCTATTGCCAGAGACAAAGAGAGAAATCTGTAACATCAACATGCTAATGctt                                         >  2:1801902/1‑148 (MQ=14)
                                                                                              cagcTATCGACCCCCCACAATCCCGCAAGCAGGAACATCAACAGCCTACTGAGGAACCAGCAACACATGCCGATCCTCTCAAATCCGCCGTCGAAGCCGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCgg                                         <  1:100959/150‑1 (MQ=2)
                                                                                                     atcccGCGACTAAGCGCGCAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCaga                                  <  2:714559/148‑1 (MQ=255)
                                                                                                       acctggcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATTCGGAACCAGAACAACAGGTCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaag                                <  2:1120148/145‑1 (MQ=255)
                                                                                                             ccAAAGCGCGCAAGCTGGCAGAGCAACAGGCTAATGCGGAACCAGAACACGTGGGCGATCCTCGGAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                          <  2:1989452/150‑1 (MQ=14)
                                                                                                                   gccgcAATCTGGAAGATCTTCAGGCTAATTAGGAACCAGAAAACCATGCGGATCCTCGCAAACCAGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                    <  2:1097586/148‑1 (MQ=14)
                                                                                                                     gggaaGCTGGAACACCACCAGGCAATTGAGGACCCAGCACCCGAGGTGGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcg                  <  2:28775/147‑1 (MQ=14)
                                                                                                                                aCAGCAACAGCCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccgc       <  1:1700375/150‑1 (MQ=255)
                                                                                                                                    cAACAGCATACTGAGGACCCAGAACACGAGGTCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTc   <  1:1148508/150‑1 (MQ=14)
                                                                                                                                     atcAGTCTAATGCCGAAGCAGAAACACAGGTCGATCCACGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCACGCAAAGCCGCCGTCg  <  2:1831145/148‑1 (MQ=14)
                                                                                                                                     aaCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCgag                                                                                <  1:1495299/72‑3 (MQ=21)
                                                                                                                                     aaCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCgag                                                                                >  2:1495299/1‑70 (MQ=21)
                                                                                                                                                  |                                                                                                                                        
GCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCG  >  NC_000913/1708506‑1708788

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