Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,709,038 T→A 8.3% V707E (GTA→GAA)  rsxC → SoxR [2Fe‑2S] reducing system protein RsxC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,709,0380TA 8.3% 221.8 / 3.3 86V707E (GTA→GAA) rsxCSoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base T (41/37);  new base A (3/5);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAGCG  >  NC_000913/1708889‑1709185
                                                                                                                                                     |                                                                                                                                                   
cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGt                                                                                                                                                     >  2:926815/1‑150 (MQ=255)
     ttGCCCGATCCAAAGCGCGCAAACTGTAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCag                                                                                                                                                <  1:1767809/150‑1 (MQ=255)
         ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGAACGCAACAAAGAACTCAAGCTGGAACAGCAAAAGGCTAATGCGGAACCagaaga                                                                                                                                            >  2:1789992/1‑150 (MQ=17)
                    cgcgcAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGACGCGGCTATTGCACGCGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCGAAATCGGTACCAGAAGCACAGCTAGATc                                                                                                                                 >  1:807169/1‑150 (MQ=17)
                      cgcAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCcg                                                                                                                               >  2:923497/1‑150 (MQ=255)
                        cAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCCCGCGCACAAGCACGCAAGCTGGAACAGCAACAGGATAATGACGTACAATAAGAACAGGTAGATCagcg                                                                                                                             >  2:2485022/1‑150 (MQ=17)
                        cAAACGGCAACAGCAACCCCCTAATGCCGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCAAAAGCAACAAAGCACGAACAGCAACAGGATAAAGCCGTACAAGAACAACATGTTAAACcgcc                                                                                                                             >  2:181216/1‑149 (MQ=19)
                              ggAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGc                                                                                                                       >  2:28539/1‑150 (MQ=255)
                                aaCAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTCCCAGAAGAACAGGTTGCTCAGCGCAACGCgg                                                                                                                     >  1:810668/1‑150 (MQ=18)
                                aaCAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCAAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCgg                                                                                                                     >  2:1551465/1‑150 (MQ=255)
                                    gCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCAAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGt                                                                                                                 >  2:360220/1‑150 (MQ=255)
                                     cAACCGGCTAATGCGGAGCCAGAATAACAGGTTGATCCGCGCAAAGCAGCAGTAGAAGCGGCTAATGCACGCGCAAACGCAAGAAACCACGCACAGAAAAAGGATAAAGAGGTACCAAAAAACCAGGTACATCACCGAACAGCGGCatat                                                                                                                >  2:2190624/1‑147 (MQ=9)
                                      aaCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTg                                                                                                               >  1:1156364/1‑150 (MQ=255)
                                         cGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGccc                                                                                                            >  1:944652/1‑149 (MQ=255)
                                          ggCTAACGCCGAGCCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcgc                                                                                                           <  1:928119/150‑1 (MQ=255)
                                             tAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGc                                                                                                        >  2:785883/1‑150 (MQ=255)
                                               aTGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCCAAGCCGCTATCGCACGCGCCACAGCACCCAAGCGGGAACAGCAACAGGCTAATGCGGACCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTa                                                                                                      <  1:1127582/150‑1 (MQ=11)
                                                tGCGGAGCCAGAACAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTAt                                                                                                     <  1:1789992/150‑1 (MQ=255)
                                                  cGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGCACGCGCAAAAGCAAAAAAACTCGAACAGCAAAAGGATAATACGGAACAAGAAAAACAATTAGATAACCGAAAAGCGGCAGTTACCGAGTATATTg                                                                                                   >  2:2371594/1‑150 (MQ=12)
                                                        cAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcgcg                                                                                             >  2:738093/1‑150 (MQ=255)
                                                                  ggTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaa                                                                                   >  2:163012/1‑150 (MQ=255)
                                                                  ggTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaa                                                                                   >  1:1953147/1‑150 (MQ=255)
                                                                      gATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaaaaaa                                                                               >  2:1078052/1‑150 (MQ=255)
                                                                       aTCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAg                                                                              >  1:1561582/1‑150 (MQ=255)
                                                                       aTCCGCGCAAAGCCGCCGTCGAAGCTGCTATTGCACGAGCCAAAGCACGCAAGCTGGAACAGCAACCGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAg                                                                              <  1:731992/150‑1 (MQ=255)
                                                                         ccGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgcc                                                                            <  2:2477703/150‑1 (MQ=255)
                                                                          ctcccacaGCCGCCGTCGAAGCCGCTATAGCCCGCTCCAAAGCACTCAACCTGCAACAGCAACTCGCAAATGCGGACCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgccg                                                                           <  2:2373825/143‑1 (MQ=17)
                                                                             gcAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGccc                                                                        <  1:1942519/150‑1 (MQ=255)
                                                                                aaGCTGCCGTCTAAGCCGCTATTGCACGCGCCAAAGCACGCAAGCGGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCcagc                                                                     <  1:1696469/150‑1 (MQ=255)
                                                                                   ccgCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCgg                                                                                                                     <  1:202646/97‑1 (MQ=255)
                                                                                    tGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGaa                                                                 >  1:2037470/1‑150 (MQ=255)
                                                                                      ccGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAgg                                                               <  1:320820/150‑1 (MQ=255)
                                                                                       cGTCGAAGCGGCTATTGCCCGCGCCCAAGCCCGCAACCGGGACCAGCAACAGGCTAATGGGGTCCCAGAAGAACATGTCGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGt                                                              <  1:1354080/150‑1 (MQ=255)
                                                                                       cGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGt                                                              >  2:174439/1‑150 (MQ=255)
                                                                                       cGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGt                                                              >  1:1696783/1‑150 (MQ=255)
                                                                                        gTCGAAGCGGCTATTGCACGCTCCAAAGCACTCAATCTGGAACAGCAAAAGGCTAATGCGGTTCCAGAAGAACAGGTTGATCCACGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGtt                                                             <  2:1989277/150‑1 (MQ=37)
                                                                                        gTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGtt                                                             >  2:2149530/1‑150 (MQ=255)
                                                                                             aGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTaaa                                                        >  2:869119/1‑150 (MQ=255)
                                                                                               cgtcctttGCGCGCGCCTACGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACg                                                      <  1:2050041/144‑1 (MQ=25)
                                                                                               cGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACg                                                      <  1:28539/150‑1 (MQ=255)
                                                                                                ggCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGa                                                     >  1:583808/1‑150 (MQ=255)
                                                                                                ggCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCCCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGa                                                     >  2:187217/1‑150 (MQ=255)
                                                                                                 gCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAg                                                    <  1:1865590/150‑1 (MQ=255)
                                                                                                  cTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAgg                                                   >  2:344266/1‑150 (MQ=255)
                                                                                                   tATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTACCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGa                                                  <  1:35442/150‑1 (MQ=255)
                                                                                                     ttGCACGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACt                                                <  2:1401561/150‑1 (MQ=255)
                                                                                                       gCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCACAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTaa                                              >  1:1875172/1‑150 (MQ=255)
                                                                                                          cgtgcCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATg                                           <  1:1080083/147‑1 (MQ=255)
                                                                                                          cgcgCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATg                                           >  1:2031107/1‑150 (MQ=255)
                                                                                                          cgcgCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCACAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATg                                           >  1:1008983/1‑150 (MQ=255)
                                                                                                           gcgcCAAAGCACGCAAGCGGGAACAGCAACAGGCTAATGCGGTTCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATgg                                          <  2:2456770/150‑1 (MQ=255)
                                                                                                              cgaacgCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGTCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTAt                                       <  2:536261/145‑1 (MQ=255)
                                                                                                                aaaGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTc                                     <  1:1843229/150‑1 (MQ=255)
                                                                                                                 aaGCACGCTATCTGGAACATCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCa                                    <  1:2371594/150‑1 (MQ=255)
                                                                                                                    cACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGaa                                 >  1:2019745/1‑150 (MQ=255)
                                                                                                                    cACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGaa                                 >  2:1376271/1‑150 (MQ=255)
                                                                                                                     aCGCAATCTGGAACAGCAACAGGCTAATGCGGATCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAt                                <  2:1697181/150‑1 (MQ=255)
                                                                                                                     aCGCAAGCTGTAACAGCAAGAGGCTAATGCGGTTCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAt                                <  2:299918/150‑1 (MQ=255)
                                                                                                                     aCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAt                                <  2:1953147/150‑1 (MQ=255)
                                                                                                                      ctcAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAta                               <  1:181216/148‑1 (MQ=255)
                                                                                                                      ctcAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAta                               <  1:189980/148‑1 (MQ=255)
                                                                                                                        cAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAtagc                             >  2:1128842/1‑150 (MQ=255)
                                                                                                                            cTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAtagctagc                         <  1:919432/150‑1 (MQ=255)
                                                                                                                            cTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAtagctagc                         <  1:2485022/150‑1 (MQ=255)
                                                                                                                                aaCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTccc                     >  1:563176/1‑150 (MQ=255)
                                                                                                                                 aCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTcccc                    <  2:328675/150‑1 (MQ=255)
                                                                                                                                 aCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTcccc                    <  1:367006/150‑1 (MQ=255)
                                                                                                                                  cAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCt                   >  2:1019086/1‑150 (MQ=255)
                                                                                                                                    gCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTta                 >  2:2241414/1‑150 (MQ=255)
                                                                                                                                    gCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAg                                                                              <  1:659423/89‑1 (MQ=255)
                                                                                                                                    gCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAg                                                                              >  2:659423/1‑89 (MQ=255)
                                                                                                                                      aaCAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTtata               <  1:1804723/150‑1 (MQ=255)
                                                                                                                                        cAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATAcc             <  2:880776/150‑1 (MQ=255)
                                                                                                                                        aaGGCTAATGCGGATCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATAcc             <  2:1981059/149‑1 (MQ=255)
                                                                                                                                          ggCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCa           >  1:1139053/1‑150 (MQ=255)
                                                                                                                                          ggCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCa           <  2:1369323/150‑1 (MQ=255)
                                                                                                                                          ggCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCa           >  1:2037580/1‑150 (MQ=255)
                                                                                                                                          ggCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCa           <  2:810668/150‑1 (MQ=255)
                                                                                                                                          ggCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCACAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCa           >  2:1200963/1‑150 (MQ=255)
                                                                                                                                           gCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCACAATAGCTAGCTCCCCTTATACCCaa          >  2:1390102/1‑149 (MQ=255)
                                                                                                                                              aaTGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAAc       <  1:941520/150‑1 (MQ=255)
                                                                                                                                               aTGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAAcc      <  1:1053626/150‑1 (MQ=255)
                                                                                                                                               aTGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAAcc      <  2:557230/150‑1 (MQ=255)
                                                                                                                                                tGCGGTACCAGAAAAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCa     >  1:891988/1‑150 (MQ=255)
                                                                                                                                                 gCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAg    <  2:2239114/150‑1 (MQ=255)
                                                                                                                                                  cGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAgc   >  2:1409311/1‑150 (MQ=255)
                                                                                                                                                   ggTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAgcg  >  2:286476/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAGCG  >  NC_000913/1708889‑1709185

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