Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,844 C→G 7.8% A642A (GCC→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,8440CG 7.8% 656.9 / 23.6 206A642A (GCC→GCG) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (106/80);  new base G (5/12);  total (112/94)
Fisher's exact test for biased strand distribution p-value = 4.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

GCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTG  >  NC_000913/1708698‑1708992
                                                                                                                                                   |                                                                                                                                                    
 cTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGc                                                                                                                                                  <  1:1342109/150‑1 (MQ=255)
 cTATGGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGc                                                                                                                                                  <  1:1539651/150‑1 (MQ=255)
   atcgCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCa                                                                                                                                                <  2:1702678/147‑1 (MQ=21)
   aTTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCa                                                                                                                                                <  1:2437521/150‑1 (MQ=255)
      gCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCTGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGaa                                                                                                                                             >  2:3909921/1‑150 (MQ=255)
       cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGaac                                                                                                                                            >  2:3737274/1‑150 (MQ=255)
       cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGaac                                                                                                                                            >  2:324970/1‑150 (MQ=255)
         cgcgCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGaacaa                                                                                                                                          <  2:722816/147‑1 (MQ=255)
          gTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGaacaac                                                                                                                                         <  1:736862/150‑1 (MQ=255)
          gTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGaacaac                                                                                                                                         <  1:737083/150‑1 (MQ=255)
          gTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGaacaac                                                                                                                                         <  1:965025/150‑1 (MQ=255)
          gTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGaacaac                                                                                                                                         <  2:2893812/150‑1 (MQ=255)
          gTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTCGAACAGCAACAG‑GCTAAAGCGGAGCCAGaactcc                                                                                                                                         >  2:2260935/1‑147 (MQ=11)
           tGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAacaaca                                                                                                                                        <  2:3002183/150‑1 (MQ=255)
            gCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAg                                                                                                                                       <  1:2726408/150‑1 (MQ=255)
               aaaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTc                                                                                                                                    <  1:81000/150‑1 (MQ=37)
                    gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAAC‑GGAAGAGCGATCATGCCGGGCAGGAAGAAAAGGTCGATca                                                                                                                               >  2:3795843/1‑149 (MQ=11)
                      gcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCcgc                                                                                                                             <  1:4022199/150‑1 (MQ=255)
                            tGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAg                                                                                                                       <  1:2649325/150‑1 (MQ=37)
                              gAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc                                                                                                                     >  2:1842777/1‑150 (MQ=255)
                              gAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc                                                                                                                     <  1:3826403/150‑1 (MQ=255)
                              gAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc                                                                                                                     >  2:3967485/1‑150 (MQ=255)
                                aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgc                                                                                                                   <  1:471961/150‑1 (MQ=255)
                                aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgc                                                                                                                   >  2:249875/1‑150 (MQ=255)
                                 cAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgcc                                                                                                                  <  1:1347722/150‑1 (MQ=255)
                                 cAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgcc                                                                                                                  >  1:3989919/1‑150 (MQ=255)
                                  aGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg                                                                                                                 >  1:1489412/1‑150 (MQ=255)
                                   gCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGTTCGGCGCAAAGCCGCCGTCGATGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACTG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGt                                                                                                                <  1:3260129/150‑1 (MQ=255)
                                    cAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                                               >  1:3394837/1‑150 (MQ=255)
                                       cAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGaa                                                                                                            <  2:233318/150‑1 (MQ=255)
                                        aGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                                                                                                           >  2:2583481/1‑150 (MQ=255)
                                          gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                                                         >  2:2862526/1‑150 (MQ=255)
                                            tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                                                       >  2:3833300/1‑150 (MQ=255)
                                            tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                                                       <  1:1218223/150‑1 (MQ=255)
                                                tcGGAACCAGAATCACATGTCGATCCGCGCAAAGCCGCCGTTGAAGCAGCTATTGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAG‑GCTAATGCCGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAtt                                                                                                   <  1:522846/149‑1 (MQ=11)
                                                  ggAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTTGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                                                                                 >  2:2907999/1‑150 (MQ=255)
                                                  ggAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                                                                                 <  1:2908942/150‑1 (MQ=255)
                                                    aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGcac                                                                                               >  1:3934618/1‑148 (MQ=255)
                                                     aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAATCTGGAACAGCAACAG‑GATAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGCCGCCGCTATTGCCCg                                                                                              >  2:2342016/1‑150 (MQ=255)
                                                     aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCg                                                                                              >  1:1995450/1‑150 (MQ=255)
                                                        agaagaACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGc                                                                                           >  2:1432084/1‑150 (MQ=255)
                                                        agaagaACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGc                                                                                           >  1:1284545/1‑150 (MQ=255)
                                                        agaagaACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGc                                                                                           >  2:116434/1‑150 (MQ=255)
                                                         gaggaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGcc                                                                                          >  2:3082394/4‑150 (MQ=255)
                                                              aCAGGTCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCTCGCAATCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                                                     <  1:3462245/150‑1 (MQ=255)
                                                               cAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcg                                                                                    >  1:3202638/1‑150 (MQ=255)
                                                                tggTAGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAG‑GCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgc                                                                                   <  1:2805053/149‑1 (MQ=11)
                                                                aGGTCGATCCGCGCAAAGTCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCAAAAGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgc                                                                                   >  2:15254/1‑150 (MQ=255)
                                                                      aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAAc                                                                             >  2:1273394/1‑150 (MQ=255)
                                                                      aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCAGTAGACGCCGCTATTGCCCGAGCAAAGGCCCTCAAAc                                                                             >  2:2092838/1‑150 (MQ=14)
                                                                        ccGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTTCCAAAGCACGCAAGCTGGAACAGCTACAG‑GCTAATGCCGAGCCAGCACAACAGGTCGATCCGCGCAAAGCCTCCGTCTAAGCCTCTATTGCCCGAGCCAAAGCGCGCAAACgg                                                                           <  1:1386601/150‑1 (MQ=14)
                                                                        ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACgg                                                                           <  2:2941371/150‑1 (MQ=255)
                                                                         cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCACACggg                                                                          >  1:2364865/1‑150 (MQ=255)
                                                                           cgcAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGaa                                                                        <  1:642599/150‑1 (MQ=255)
                                                                           cgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGaa                                                                        <  1:1006924/150‑1 (MQ=255)
                                                                             cAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                                                      <  2:489912/150‑1 (MQ=255)
                                                                              aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAg                                                                     <  1:4134725/150‑1 (MQ=255)
                                                                                   cgccgTCGTAGCAGCTATTGCCCGTGCCAACGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAAcc                                                                <  2:426283/150‑1 (MQ=255)
                                                                                   cgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑TCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAAcc                                                                >  1:3950612/1‑150 (MQ=255)
                                                                                   cgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAAcc                                                                >  2:1586546/1‑150 (MQ=255)
                                                                                    gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGCACAGCAACCg                                                               >  2:3588387/1‑150 (MQ=255)
                                                                                     ccgTCGAAGCCGCTATTGCCCGTGCCATAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCACGAGCCAAAGCGCGCAAGCGGGAACAGCAACCgg                                                              >  2:1921025/1‑150 (MQ=255)
                                                                                     ccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCACCCgg                                                              >  2:2410381/1‑150 (MQ=255)
                                                                                       gTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCt                                                            <  1:4236939/150‑1 (MQ=255)
                                                                                           aaGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATg                                                        >  1:4019269/1‑150 (MQ=255)
                                                                                            aGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGc                                                       <  2:1219053/150‑1 (MQ=255)
                                                                                             gCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCg                                                      >  1:3663336/1‑150 (MQ=255)
                                                                                               cGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                                                    <  1:1586817/150‑1 (MQ=255)
                                                                                               cGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                                                    >  1:1252460/1‑150 (MQ=255)
                                                                                                  tATTGCCCGTGCCAACGCACGCTAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGcc                                                 <  1:2151865/150‑1 (MQ=255)
                                                                                                  tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGcc                                                 <  2:3989919/150‑1 (MQ=255)
                                                                                                    ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCag                                               >  2:1371962/1‑150 (MQ=255)
                                                                                                       cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                                            >  2:1666293/1‑150 (MQ=255)
                                                                                                       cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                                            >  1:910782/1‑150 (MQ=255)
                                                                                                       cccGGGCCAACGCTCGCAAGCTGGAACAGCATCAG‑GCTACTGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATGGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                                            <  2:4238249/150‑1 (MQ=14)
                                                                                                        ccGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                           >  2:2917466/1‑150 (MQ=255)
                                                                                                        ccGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGAGCAGCGCAAAGCCGCCGTCGAAGCCGCTCTTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                           >  2:3159751/1‑150 (MQ=255)
                                                                                                          gTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAAc                                         >  1:3571074/1‑150 (MQ=255)
                                                                                                            gCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg                                       <  2:2873245/150‑1 (MQ=255)
                                                                                                            gCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg                                       >  1:3120342/1‑150 (MQ=255)
                                                                                                             ccAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAgg                                      >  1:4138619/1‑150 (MQ=255)
                                                                                                               aaaGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGtt                                    >  2:51738/1‑150 (MQ=255)
                                                                                                                aaGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTg                                   >  1:1784600/1‑150 (MQ=255)
                                                                                                                aaGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTg                                   >  1:1119164/1‑150 (MQ=255)
                                                                                                                aaGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAGGCGCGAAAACGGGAACAGCAGACGGCGCTTGCGGGGCCAGAAGAAAAGGTTg                                   >  2:326739/1‑150 (MQ=17)
                                                                                                                  gCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGAt                                 >  1:1085967/1‑150 (MQ=255)
                                                                                                                  gCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGAt                                 >  2:1200124/1‑150 (MQ=255)
                                                                                                                    aCGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc                               >  1:218618/1‑150 (MQ=255)
                                                                                                                      gCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc                             <  1:2583481/150‑1 (MQ=255)
                                                                                                                            tGGAACAGCAACAG‑GCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAg                       <  2:2149379/150‑1 (MQ=255)
                                                                                                                             ggAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGc                      >  2:1286452/1‑150 (MQ=255)
                                                                                                                              gAACAGCAACAG‑GCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCt                     <  2:3436189/150‑1 (MQ=255)
                                                                                                                                 cAGCAACAG‑GCTAATGCCGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAACCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGcc                  <  2:795876/150‑1 (MQ=255)
                                                                                                                                    cAACAG‑GCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCg                                                      <  2:2936255/111‑1 (MQ=255)
                                                                                                                                    cAACAG‑GCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCg                                                      >  1:2936255/1‑111 (MQ=255)
                                                                                                                                       cAG‑GCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGaa            <  1:1855385/150‑1 (MQ=255)
                                                                                                                                        cacgCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGaa            <  2:3768150/147‑1 (MQ=255)
                                                                                                                                           tatAATGCTGAGCCAGAACTCCCGGTCGATGAGCGCAAAGCCGCCGTCGAAGCCGCTTTTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGATAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCg         <  2:3560228/148‑1 (MQ=17)
                                                                                                                                             tAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGACCAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGc       >  1:2665827/1‑150 (MQ=255)
                                                                                                                                                 gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAtt   <  1:3837859/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAG‑GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTG  >  NC_000913/1708698‑1708992

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑