Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,662 C→G 9.9% Q582E (CAA→GAA)  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,6620CG 9.9% 135.8 / 5.5 58Q582E (CAA→GAA) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (16/34);  new base G (4/2);  total (20/38)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCC  >  NC_000913/1708514‑1708808
                                                                                                                                                    |                                                                                                                                                   
aGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGaaca                                                                                                                                                    <  1:1572498/150‑1 (MQ=255)
 gcgcAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAAACGGCTAATCCGCAACCAGaacac                                                                                                                                                   >  1:1731049/1‑149 (MQ=255)
   gcAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAacaaca                                                                                                                                                 <  1:1337924/150‑1 (MQ=255)
    cAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAg                                                                                                                                                >  1:1752443/1‑150 (MQ=255)
     aaGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAgg                                                                                                                                               <  2:1032017/150‑1 (MQ=255)
       tccAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACAGCCGTAGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGATCGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑c                                                                                                                                            >  1:946833/2‑150 (MQ=255)
        ccAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGCACAACAGGT‑cc                                                                                                                                           >  1:1703678/1‑149 (MQ=255)
         cAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAGTGCGGCACCAGATCAACAGGA‑CGa                                                                                                                                          >  2:559957/1‑150 (MQ=255)
           gagCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCTCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATc                                                                                                                                        <  1:1461572/150‑1 (MQ=255)
           gagCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATc                                                                                                                                        <  1:922829/150‑1 (MQ=255)
               gtatcAGGCAGACCTGCAGCAAACTAACGACGAAGCAACCGTTGCTGATCCACGTAAAACGGCCGTTGAAGCAGCTATCGCCCGCGCCATAGCGCGCAAGCTGGAACAGCATCAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCcgcg                                                                                                                                    <  2:1238868/146‑1 (MQ=25)
               gAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCcgcg                                                                                                                                    >  1:1698208/1‑150 (MQ=255)
               gAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCcgcg                                                                                                                                    <  1:1594703/150‑1 (MQ=255)
                aaaCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACCGCAACAGGCTACTGCGGAACCAGACCAACAGGT‑CGATCCgcgc                                                                                                                                   >  1:571632/1‑150 (MQ=255)
                aaaCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGG‑CGATCCgcgc                                                                                                                                   >  1:1741679/1‑150 (MQ=255)
                 aaCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCa                                                                                                                                  <  1:1772748/150‑1 (MQ=255)
                   cAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCaaa                                                                                                                                >  2:1765851/1‑150 (MQ=255)
                     ggCAGACCTTCAGCCAACTCCCGACGCAGAAACCGTGGCTGATCCTCGCAAACCTGCCGTGGAAGCAGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGT‑CGATCCGCGCAAAgc                                                                                                                              <  2:90515/150‑1 (MQ=14)
                         gAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAgccgcc                                                                                                                          <  2:1009128/150‑1 (MQ=255)
                              gcagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACAAGAACAACAGAT‑CGAAACACGAAAAgacgcaggagc                                                                                                                     >  1:1337473/1‑144 (MQ=25)
                                 gcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGACCCAGAACAACAGGA‑CGATCCGCGCAAAGCCGCAGTCGAAGc                                                                                                                  >  1:1032501/1‑150 (MQ=255)
                                 gcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGAGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCAGCCGTCGAAGc                                                                                                                  >  1:1639310/1‑150 (MQ=255)
                                    aaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                                                               <  1:875747/150‑1 (MQ=255)
                                     aCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                                                                              <  1:1652897/150‑1 (MQ=255)
                                               cagcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg                                                                                                    <  2:778749/150‑1 (MQ=255)
                                               cagcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg                                                                                                    <  1:1037096/150‑1 (MQ=255)
                                                   aaCCGTTGCTGAGCCCCGTAAAACTGCCGTTGTAGCAGCTATCGCCCGCGCCAACGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGCACCACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaa                                                                                                <  2:307632/150‑1 (MQ=17)
                                                     aacatGCCGATCCCCGCAAACCTGCAGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCTCAAGCTTGAACAGCATCAGGCTAATGCGGTACCAGAAGAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAg                                                                                              <  2:1650275/146‑1 (MQ=11)
                                                       gTTGCTGATCCACGTAAAACTGCCTGTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcg                                                                                            <  1:436907/150‑1 (MQ=255)
                                                       gTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcg                                                                                            <  2:1752443/150‑1 (MQ=255)
                                                         tGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGGCAAAgagcg                                                                                          >  1:1650275/1‑150 (MQ=255)
                                                         tGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                                          <  1:683630/150‑1 (MQ=255)
                                                             taTCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                                                      <  1:127404/149‑1 (MQ=255)
                                                                ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCAGGTAGGAAAAGG‑CGCTCCGAGCGAAACAGCCGAAGAAGCAGGTTTTGCCCGTGCAACAGCGCGCaatata                                                                                   >  2:1492457/1‑146 (MQ=12)
                                                                ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                   <  1:705294/150‑1 (MQ=255)
                                                                 cACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTgg                                                                                  <  2:1794308/150‑1 (MQ=255)
                                                                 cACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTgg                                                                                  <  2:1354900/150‑1 (MQ=255)
                                                                  ccGTAAAACTGCCTTTGAAGCAGCTAGCTCCCGCGCCAAAGCTCGCAAGCTGGAACATCATCAGGCTAATTCGGAACCAGAACAACATGG‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                                 <  1:1311565/149‑1 (MQ=17)
                                                                    gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                               <  1:1765851/150‑1 (MQ=255)
                                                                         aCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGT‑AGATCCGCGCAAAGCCGCAGTAGCAGCGGCTATTGCCCGCGCAAACGCCCGCAACCACGCACAGCac                                                                          >  1:839106/1‑149 (MQ=14)
                                                                           tGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAA‑AGGTAGGATCAGCGCACGGCTGACGTCGAGTCCGGAGTTGAGCGAGGCAAAGCGCGGAAGCTGGAAAAGAAACa                                                                        >  1:813606/1‑150 (MQ=9)
                                                                             ccGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCACGCTGGCACAGCAACAgg                                                                      >  1:1405249/1‑150 (MQ=255)
                                                                                 tGATGCAGCTCTCTCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                                                                  <  1:1207031/150‑1 (MQ=255)
                                                                                  ctagaagcTATCTCCCGTCCCAAAGCTCGCAAGCTGGAACATCATCAGGCTAATCAGGACCCAGCACAACAGGT‑CGATCCCCGAAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAAt                                                                 <  1:326706/145‑1 (MQ=11)
                                                                                     gcagcaATCGCCCGCCCCTAAGCGCGCAAGCTGGAACATCAACAGGCTAATGCGGAACCAGCACAACAGGG‑CGATCCCCGCAAAGCCGCCGTCGAAGCCGCTATGGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCg                                                              <  1:337334/150‑1 (MQ=14)
                                                                                              gCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGccc                                                                                                       >  2:649251/1‑100 (MQ=37)
                                                                                               cccGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCTGAAGCAGAACAACAGGG‑CGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaag                                                    <  2:1440803/150‑1 (MQ=14)
                                                                                                    gcCTAAGCTCGGAAGCTGGAACATCAACAGGCTACTCCGGAAGCAGAACAACTTGT‑AGATCCGCGCAAATCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg                                               <  2:259689/150‑1 (MQ=14)
                                                                                                     ccAAACCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                              <  2:1566170/150‑1 (MQ=255)
                                                                                                      cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                                                                                         >  1:1869890/1‑90 (MQ=33)
                                                                                                       cacccgcgcAAGCTGGAACAGCTACAGGATAATGCGGAACCAGAAACCGAGGG‑CGATCCCCGGAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                            <  2:62887/146‑1 (MQ=14)
                                                                                                       aaatgtcgcAATCTGGAACATCAACAGGATAATTCGGAACCAGAAAAACATGC‑CGATCCTCGCAAATCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAATCAGAAGAACAGGTc                                            <  1:1253784/144‑1 (MQ=14)
                                                                                                              gcAAGCTGGAACAGCAACAGGCTAATGAGGAACCAGAACAACAGGT‑CGATCCTCGTAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACATCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc                                     <  1:1679040/150‑1 (MQ=14)
                                                                                                                 atcTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCa                                  <  1:1681520/148‑1 (MQ=255)
                                                                                                                    tGGAACAGCTCCAGGCAACTTCGGAACCAGAACAACAGGC‑CGATCCTCGCAAACCCGCCGGCGCAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAg                               <  2:1725055/150‑1 (MQ=14)
                                                                                                                     ggAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGT‑CGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGGAGGCTTGCGCAGCAGACGTATTAGGAGGAGCCATAAGAACAGGTAGATACGCGAAAAga                              >  2:1679664/1‑149 (MQ=2)
                                                                                                                           gCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                                >  2:480163/1‑78 (MQ=21)
                                                                                                                           gCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                                <  1:480163/78‑1 (MQ=21)
                                                                                                                             aaCAGGCTAATTCGGAACCAGAACAACAGGT‑CGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCg                      <  2:1048366/150‑1 (MQ=255)
                                                                                                                                aGGCTAATCCGGACCCAGAACAACAGGC‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                   <  2:577992/150‑1 (MQ=255)
                                                                                                                                                 gAACAACAGGT‑CGATCCTCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc  <  2:1396984/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGT‑CGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCC  >  NC_000913/1708514‑1708808

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