Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,814 A→G 6.8% A632A (GCA→GCG 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,8140AG 6.8% 192.0 / 2.6 75A632A (GCA→GCGrsxCSoxR iron‑sulfur cluster reduction factor component, putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base G (2/4);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGA  >  NC_000913/1708576‑1709044
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                        
aTCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAACTGGAACAGCAACAGGCTAATGAGGAAACAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCACAAGCa                                                                                                                                                                                                                                          >  1:500209/1‑239 (MQ=255)
  ccgcatAAAACTGCGGTAGAAGCAGCTATCGCACGCGCCAAAGCGCGCAAGCTGGAACAGCCACAGGCTAATGCGGAACCAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACg                                                                                                                                                                                                                                        <  1:418798/234‑1 (MQ=255)
  ccACGTATACCTGCCACGAAATCAGCTATCGCCCGCTCCAAAGCGAGCACGCTGGACCAGAACCTCAGTAATGCGGAACCAACACCAGAGGTCGATCCGCGTAAAGCCGCCGTCGAAGCCGCTATCGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACg                                                                                                                                                                                                                                        <  2:661028/239‑1 (MQ=14)
  ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCAc                                                                                                                                                                                                                                         >  1:1139762/1‑238 (MQ=255)
  ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAc                                                                                                                                                                                                                                         >  1:661028/1‑238 (MQ=255)
  ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACg                                                                                                                                                                                                                                        <  1:681044/239‑1 (MQ=255)
  ccACGAAAGACTGCGGGTGAGGCAGCCAGCGACCGCGCCAACGGGGGGAAGCTGGAACAGCAACAGCCAAATGCGGACCCAGAACAACTGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACg                                                                                                                                                                                                                                        <  2:1145367/239‑1 (MQ=14)
     cGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCaa                                                                                                                                                                                                                                     >  1:389008/1‑239 (MQ=255)
     cGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCACGTGCCAAAGCACGCAAGCTGGAACCGCAACCGCATACACCGGAAACGGAGCAAAAGATGGACTCGGGCAAAACCGCCGTGGAAGCCGCAACTGCAGATGACAAAACAAGCaa                                                                                                                                                                                                                                     >  2:308384/1‑239 (MQ=14)
          aaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                                                                                                                                                                                                                               >  1:876486/1‑239 (MQ=255)
               ccGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGCCTACTGCGGCACCAAAAGAACAGGCCGATCCGCCCAAAGCCGCCGTCGCAGCGGCCATTCCCCGGCCCAAACACCGCAAGCTGGAGCa                                                                                                                                                                                                                           >  2:1054818/1‑238 (MQ=14)
                 gtggAAGCAACTACCACCCGCCCACAAGAGCGCAAACTGGAACAGCGACAGGCTAATCCAGACGCAGAAGAACAGGTCGATCCGCGCACAGCCGCCGTTGAAGACGCTATCGCCGGCGCCAAAGCACGAAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAAAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCa                                                                                                                                                                                                                        <  2:1123289/236‑1 (MQ=14)
                   tGAAGCAGCTCGCGCCCGCGCCAAAGCGGCCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCaa                                                                                                                                                                                                                       <  1:368835/238‑1 (MQ=255)
                   tGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAGACCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCaa                                                                                                                                                                                                                       <  1:121473/238‑1 (MQ=255)
                                  ccGCGACAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGc                                                                                                                                                                                                       <  1:296259/239‑1 (MQ=255)
                                   cgcgCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAGcc                                                                                                                                                                                                      >  1:819908/1‑239 (MQ=255)
                                     cgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAg                                                                                                                                                                                                    >  1:993375/1‑239 (MQ=255)
                                                              cAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGTCGATCCGCGCAAAg                                                                                                                                                                           >  1:26167/1‑238 (MQ=255)
                                                              cAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACg                                                                                                                                                                                                                                        >  1:798518/1‑179 (MQ=255)
                                                              aaaCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACg                                                                                                                                                                                                                                        <  2:798518/178‑1 (MQ=255)
                                                                       aaTGCGGAACCAGAACAACGGGTCGATCCGCGAAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAAGACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGa                                                                                                                                                                 <  1:1148396/239‑1 (MQ=255)
                                                                              aaCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACCGGTCGATCCGCGCAAAGCCGCCGTCGCCGCGGCTATTGCCCGCGCTCAAGCCCGCAAACTGGAACAGCAACAGGCTGATTACGAGCCCGAAGAACCGGTCGAGCCGCGCAACGCCCCCGGCGCAGccaca                                                                                                                                                          >  2:874093/1‑236 (MQ=14)
                                                                              aaCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCTCCAAACCCCGCAAGCCGGACCAGCAACCGGCCAATGCCGAGCCCGACGACAAGGTCGATACGCGCAAAGCCGCCGTCGAAGCGGca                                                                                                                                                          >  2:368375/1‑238 (MQ=11)
                                                                              aaCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGACGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                                                                                                                                                                                                  >  1:570265/1‑200 (MQ=255)
                                                                               aCCCGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                                                                                                                         <  1:157573/239‑1 (MQ=255)
                                                                                  aGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGATCCAGAACACAGGCTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTagt                                                                                                                                                       >  2:1093172/1‑236 (MQ=255)
                                                                                     acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                                                                                                                    <  1:1122394/238‑1 (MQ=255)
                                                                                           ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAAAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGCGCCa                                                                                                                                             >  1:420682/1‑239 (MQ=255)
                                                                                           ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTGCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGCCGAAGCCGCTacacaccgacccg                                                                                                                                             >  1:162676/1‑227 (MQ=14)
                                                                                             acGATCCGCGAAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCAAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCa                                                                                                                                             <  1:532750/236‑1 (MQ=255)
                                                                                                 tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAATCCGCTATTGCCCGTGCAAAAGAACGCAAGCAGGAACAGCAACAGGCTAATGCCGAGCAAGAACAACAGGTCGATCCGCGCAAAg                                                                                                                                                                           >  1:188064/1‑203 (MQ=255)
                                                                                                  ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgcg                                                                                                                                      <  1:968480/239‑1 (MQ=255)
                                                                                                       cAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCGGCTATTGCCGACGCCAAAGCGCGCaaaa                                                                                                                                 >  1:142335/1‑238 (MQ=255)
                                                                                                             cgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                                                                           <  1:1093172/239‑1 (MQ=255)
                                                                                                                   cGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCAATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCaa                                                                                                                      <  1:230102/238‑1 (MQ=255)
                                                                                                                         cGCTATTGCCCATGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGAGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGc                                                                                                                 <  1:1098349/237‑1 (MQ=255)
                                                                                                                                  ccGGGCCAAAGCGCGCAAGCTGGAACAGAAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGc                                                                                                      <  1:320918/239‑1 (MQ=255)
                                                                                                                                     tGCCAAAGCGCGCAAGCTGGGACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGGAGCCGCGATTGCCCGGGCCAAAGCACCAAACAGCGAACAGCAACAGGCTAATGCCGAACCCCAAAAACCGGTCGATCCGCGCAAAGCCGCCGGTGCAGCAGCTATTGCACGAGCCAAAGGCAGCAAGCGGGGAGAGCAACCGCCtctatcggagcccg                                                                                                   >  2:182591/1‑226 (MQ=11)
                                                                                                                                     tGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAACAGGACCAGCAACCGGCTAAAGCTGAGCCat                                                                                                   >  2:1109679/1‑238 (MQ=255)
                                                                                                                                        cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGACCAGGTCGATCCGCGCAAAGCGGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGCCTAATGCGGAGCCagaa                                                                                                 <  1:903124/238‑1 (MQ=255)
                                                                                                                                        cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTCGACGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                                                                                                >  1:1125459/1‑239 (MQ=255)
                                                                                                                                         aaaGCGGGCAAGCTGGAACAGTAACAGGCTAATGCGGAACCAGAGACACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCTGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                                                                               <  1:761923/239‑1 (MQ=255)
                                                                                                                                                              cAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCCGAAGAACAGGTTGATCCGCGCaact                                                                           >  1:236732/1‑236 (MQ=255)
                                                                                                                                                                       aaTGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGa                                                                 <  1:368375/239‑1 (MQ=255)
                                                                                                                                                                        aTGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGAAAACGGGAACAGCAACAGGCTAATGCCGAGCCAGAATCACAGGTCGATCCGCGCAAAGCCGAAGACGAACCCGCTATTGCCGGCGCCACAGCGCCCAAAAGGGAACAGCAACCAGCCAGAGCCGAGCCAGGAGAACAGGTTGAACCGCGCCAATCTGCCggcca                                                                 >  2:1095070/1‑234 (MQ=14)
                                                                                                                                                                              aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGCGCCAACGCACACAAGCAGGCACACCAACAGGCCAATGCCGAGCAAGAACAGCAGGTCGATCCGGGCAAAGCCGCTGGTGAAGCAGCTCTGGCACGCGCCAAAGCGCGCACGACACAACAGCAACCGGCTCATGCTGCGACAGAACAAAACGCTCATCCGcacatcgcccgagacgcacacact                                                          >  2:33670/1‑219 (MQ=2)
                                                                                                                                                                              aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGCGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCGGCTATTGCCCGACCCAAAGCGCGCAACCCGGAACAGCAACCGGCAGATGCGGAGCCAGAAAAACAGGTCGATCCCCGCAACGAGACCCTCGAggcggcc                                                          >  2:895725/1‑238 (MQ=14)
                                                                                                                                                                              aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTACAGCCGAGGCCAAAAAACAGGTTGGTCGGCGCAACGCTGCCGttgcacccgcc                                                          >  2:900843/1‑229 (MQ=14)
                                                                                                                                                                                    aagaaCGGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCa                                                    <  1:1003544/239‑1 (MQ=255)
                                                                                                                                                                                     agaagaGGTCGAGCCACGGAAAACCGCCGTGGAAGCCGCTAGTACGCGTGCCACAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCa                                                    <  1:209239/238‑1 (MQ=17)
                                                                                                                                                                                     agaaCAGGACGACCCGCGCAAAACGGCCGTCGAAGCCGCTCTTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCa                                                    <  1:375740/238‑1 (MQ=255)
                                                                                                                                                                                            gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCaaa                                           >  1:65451/1‑239 (MQ=255)
                                                                                                                                                                                            gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCACGCAAACTGGAACAGCAACAGGCTAATGCCGAGCCAGAAAAACAGGTCGATCCGAGCAAAGCGGCCGTCGAGGCGGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACCGCAAACGGCTAATACGGAGCCAGAAGAACAGGTTGATCCGAGCAAACCAGCAGTCACAACGGCAATTGCCcgcgacc                                             >  2:775894/1‑235 (MQ=14)
                                                                                                                                                                                                aTCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCCCGCAAGCCGGCACAGCAACAGGCTAATGCAGCGCCAGAAGAACAGGGCGATCCGCGCAAACCGGCCGTCGTACCGGcgtttcgccgggcccagac                                         >  2:467074/1‑220 (MQ=14)
                                                                                                                                                                                                 tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCa                                        <  1:1109679/238‑1 (MQ=255)
                                                                                                                                                                                                 tCCGCCAAACGCCGCCGTCGAAGCCGCTATTCCCCGTGCCAGAGCACGCAAGCTGCAACAGCAACAGGCTAATCCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                                                                                                        <  1:357181/174‑1 (MQ=14)
                                                                                                                                                                                                    gcgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGTCGATCCGCGCAAAGCGGCCGTCGACGCGGCTATTGCCCGCGCCCAAGCGCGCAAACGCGCACACCAACAGGCTACTGCGGATCCAGAAGAACCAGCGGATCCCCGCAAAGCTGCCGGCGAAGCGGCTATTGCCCCCGCCAAAGCCCgaa                                    >  2:980046/1‑237 (MQ=14)
                                                                                                                                                                                                    gcgaAAAGCCGCCGTGGAAGCAGGTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAGCAGGCGAATGGTGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACg                                      <  1:308384/237‑1 (MQ=17)
                                                                                                                                                                                                     cgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACg                                                                                                                                                                                                                                        <  2:310443/44‑1 (MQ=38)
                                                                                                                                                                                                      gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCCCGCGCCAAAGCACGCAAg                                  >  1:477571/1‑239 (MQ=255)
                                                                                                                                                                                                      gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCACAGCCAGCaa                                   >  1:958329/1‑238 (MQ=255)
                                                                                                                                                                                                          aGCCGCCGTCGAAGCCGCTATGGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                              <  1:1054818/239‑1 (MQ=255)
                                                                                                                                                                                                             cgcagTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTGATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAAc                           <  1:42409/239‑1 (MQ=255)
                                                                                                                                                                                                               ccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACa                          >  1:580248/1‑238 (MQ=255)
                                                                                                                                                                                                                 gTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCa                       <  1:22875/239‑1 (MQ=255)
                                                                                                                                                                                                                         cGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCt                <  1:463422/238‑1 (MQ=255)
                                                                                                                                                                                                                          gCTATTGCCAGTCAAAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTaa              <  1:329285/239‑1 (MQ=255)
                                                                                                                                                                                                                            tATTGACCGTGCCAAAGCCAGCAAACTGGAACAGCCACAGGCAAATACCGAAACAGAACACCAGGACGATCCGCGCAAATCCGCCGTTGAAGCCGCTATTGACCGAGCCAAAGCGCGCAAACGGGAACAGCAACTGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATg            <  1:963742/239‑1 (MQ=14)
                                                                                                                                                                                                                                  ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTAc      <  1:1149463/239‑1 (MQ=255)
                                                                                                                                                                                                                                   cGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCCAAGCACGCAACCTGGCACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATACGCGCAAAGCCGCCGACGACCCCGACCTGCCCCGCGCCTACGCCCGCAAGCTGGagccgccacgcgcagacacgcaag      >  2:1137714/1‑215 (MQ=11)
                                                                                                                                                                                                                                     tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGCGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGCCGAAGCGGCTATTGCGCGCGCTCACACACACAACCCTCAACAGCAACCAGCTaaatctgtcacac   >  2:971187/1‑228 (MQ=14)
                                                                                                                                                                                                                                     tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCATCGGCTATTGCACGCGCCAAAGCCTGCAAGCTGGCACCGCAATAGGcactagaggtctcat   >  2:604478/1‑225 (MQ=17)
                                                                                                                                                                                                                                     cgCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGCACAGCAACCGGCTAACGCGGAGCCAGAAGAACAGGGTGATCCGCGCAAAGCAGCAGTCGCCGCGGCTATTCCCCGCGACAAGGCCCGAAAACTGGCCCAGACACAGTATAATGAGGCACCa  >  2:111777/2‑238 (MQ=14)
                                                                                                                                                                                                                                      gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCaga  >  1:75369/1‑239 (MQ=255)
                                                                                                                                                                                                                                           aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCa                                                                                                                                    <  1:954066/104‑1 (MQ=255)
                                                                                                                                                                                                                                           aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCa                                                                                                                                    >  2:954066/1‑104 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                        
ATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGA  >  NC_000913/1708576‑1709044

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: