Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,706,969 C→T 6.7% N17N (AAC→AAT)  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,706,9690CT 6.7% 166.2 / 2.4 60N17N (AAC→AAT) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base T (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

ACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAA  >  NC_000913/1706826‑1707110
                                                                                                                                               |                                                                                                                                             
aCTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCCGAAAAAATAAAATCTGGGATTTCAAcggcgg                                                                                                                                          >  2:1327128/1‑149 (MQ=255)
 cTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTTCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAAAcgcggc                                                                                                                                         >  2:493617/1‑149 (MQ=255)
   gtaTCTCGTTGCAACCGGTCGCTGAAACACCTGGCTCCTGGAGATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAGGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGcat                                                                                                                                       <  1:383131/147‑1 (MQ=255)
         cGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCatccatcc                                                                                                                                 >  2:93019/1‑149 (MQ=255)
             gCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCCTCCATCCCCCg                                                                                                                             >  1:428601/1‑149 (MQ=255)
               aaCCGGTCGCAGAAACACCTGACTCCTGGATATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCCGAAAAAATAAAAACTGGGATTTAACCGGCGGCCTCCATCCACCGga                                                                                                                           >  1:758753/1‑149 (MQ=255)
                           aaaCACCTGACTCCTGGAAATGGGATCTGAACCCCATTCCCGTGCGTATCATTCCCGTGGAACACCAGGCTTAAGTTATTCTCTGCCTTCAGAAAAAATAAAATCCGGAAATTCCACGGCGGCAGCCACCCACCGGAGATCGAAACCCa                                                                                                               >  2:888829/1‑149 (MQ=255)
                           aaaCACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCACCGGCGGCATCCATCCACCGGAGATGAAAACCCa                                                                                                               >  1:357213/1‑149 (MQ=255)
                           aaaCACCCTACTCCTCGAAAGCGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCa                                                                                                               <  2:655834/149‑1 (MQ=255)
                                    aCTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGAGGAAAACCCAGTCCAACgg                                                                                                      >  2:898131/1‑149 (MQ=255)
                                        cTGGAAATGGGATCTGAACACCATTCCCGCGCGTATCCTTCCCGTGGAACACCATGCTTAAGTTACTCTCTGCATTCAGAAAAAATAAAATCTGGGATGTCAACGGCGGCATCCATCCACCGGAGACGAAAAGCCAGTCGCACGGTaga                                                                                                  >  2:621781/1‑147 (MQ=255)
                                        cTGCAAAGGGGATCTGAACCCCATTCCGGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTaca                                                                                                  <  2:1221625/149‑1 (MQ=255)
                                               tGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCGATGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTgc                                                                                           >  2:734343/1‑149 (MQ=255)
                                               tGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTAGTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAATGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTgc                                                                                           >  2:734241/1‑149 (MQ=255)
                                                   aTCTGAACACCATTCCCGTTCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCa                                                                                       <  1:528675/149‑1 (MQ=255)
                                                      tGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGt                                                                                    >  1:1232914/1‑149 (MQ=255)
                                                      tGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGt                                                                                    >  1:1114279/1‑149 (MQ=255)
                                                        aaCACCATTCCCGTGCGTATCATTCCCGCGGCACACCATGCCTAAGCTACTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGCTGAAAACCCAGCCCCCCCCTCCCCCCCTCCGCCAGGCAc                                                                                  >  2:165341/1‑149 (MQ=255)
                                                              aTTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTgg                                                                            <  1:93019/149‑1 (MQ=255)
                                                              aTTCCCGTGCGTATCATTCCCGTGCAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTgg                                                                            >  1:1081520/1‑149 (MQ=255)
                                                               ttCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGCAAAACTAAAATCTGGGATTTCAACGGAGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTCACCCTGgc                                                                           >  1:300723/1‑149 (MQ=255)
                                                                  ccGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCa                                                                        >  1:1193078/1‑149 (MQ=255)
                                                                    gTGCGTATCATCCCCGTGGAACACCATGCTTCATTTTATACCTTCATTTAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCGGTCCAACGGTACACAACTGCGCCAGGGGTACCTGGCGCAgg                                                                      <  2:305241/149‑2 (MQ=255)
                                                                     ggggTAGCATCCCCGTGGAACACCATGCTTAAGTTACTCCCTGCCTTCCGAAAAAATAAAATCTGGGAGTGCAACGGCGGCAGCCTGGCTCCGGTTTTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCg                                                                     <  1:345578/146‑1 (MQ=255)
                                                                               tCCCGTGGAACTCCATGCTTAAGTTATGCTCTGCATTCAGAAAAAATAAAATCGGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTAtt                                                           <  2:822266/149‑1 (MQ=255)
                                                                               tCCCGTGGAACATCATGCTTAAGTGATTCTCTGCATTCTGAAAAAATAAAAGCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTAtt                                                           <  2:624997/149‑1 (MQ=255)
                                                                                 ccGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACGCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTCGCGCCCCGTTTTGTTGTTcc                                                         >  2:424518/1‑149 (MQ=255)
                                                                                 ccGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCCGGCGCAGCGGTTTGTTATTcc                                                         >  2:1186385/1‑149 (MQ=255)
                                                                                  cGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCCCCGTTTTGTTATTccc                                                        >  2:267030/1‑148 (MQ=255)
                                                                                  cGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCa                                                        <  1:1327128/149‑1 (MQ=255)
                                                                                   gTGGAACACCATGCTTAAGTTATGCCCTGCATTCAGAAAAAAGAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCAc                                                       <  1:206896/149‑1 (MQ=255)
                                                                                       cacacCATGCTTAAGTTCTTCTCTGCATTCAGAAAAAATAAAAGCTGGGCTTTCAACGGCGGCATCCATCCACCTGAGATGACAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGaa                                                   <  2:464622/148‑1 (MQ=255)
                                                                                       aaCACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGaa                                                   >  2:181850/1‑149 (MQ=255)
                                                                                        acacCCTGCTTAAGTCATTCTCTGCATTCCGAAAAAAAAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGATATGAAAACCCAGTCCAACGGTACTCCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGaaa                                                  <  1:506545/149‑1 (MQ=255)
                                                                                            cATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGc                                              <  1:1028876/149‑1 (MQ=255)
                                                                                               gCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCTGCGTTTTGTTATCCCACTGAAACAGCata                                           <  1:916648/149‑1 (MQ=255)
                                                                                                  tAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTg                                        <  1:181850/149‑1 (MQ=255)
                                                                                                    aGTTATCCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGgc                                      <  1:583481/149‑1 (MQ=255)
                                                                                                     gTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGAAGCGTGTTGTTATTCCGCGCACACGGCATATTGgcg                                     >  2:845862/1‑149 (MQ=255)
                                                                                                      gtATTCTCTGCATTCGGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATCGAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGgcgc                                    <  2:1086725/148‑1 (MQ=255)
                                                                                                           ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAATGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCTTTGAAACAGCATATTGGCGCTGAAg                               <  1:734241/149‑1 (MQ=255)
                                                                                                           ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAATGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCCCTGAAACAGCATATTGGCGCTGAAg                               <  1:734343/149‑1 (MQ=255)
                                                                                                           ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAg                               <  2:1180027/149‑1 (MQ=255)
                                                                                                           ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCACCCCCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAg                               <  1:757399/149‑1 (MQ=255)
                                                                                                             ctGCATTCAGAAAAAATATCCGCAGGAAGTCCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGt                             <  2:300723/149‑1 (MQ=255)
                                                                                                               gCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAACCAGGACAGGGGCGGTGAAggggc                           >  1:1110430/1‑146 (MQ=255)
                                                                                                               gCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGGAGGTGa                           >  1:12260/1‑149 (MQ=255)
                                                                                                                  ttCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGtt                        >  1:392211/1‑149 (MQ=255)
                                                                                                                     aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                     >  1:327897/1‑149 (MQ=255)
                                                                                                                     aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                     >  1:735672/1‑149 (MQ=255)
                                                                                                                     aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                     >  1:541394/1‑149 (MQ=255)
                                                                                                                       aaaaaaTAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGGGAATGAAACGCAGCCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCCCTGAAACAGCATATGGGCGCTGGAGGTGAGTtgtggg                   >  1:706222/1‑147 (MQ=255)
                                                                                                                       aaaaaaTAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGGGCAGCGTCTTGTTCTTCCACTGACCCGGCAGGCTGGCGCTGAAGGGGGGGtgtggg                   >  2:972791/1‑147 (MQ=255)
                                                                                                                         aaaaaaaaatCTGGGATTTCAACGGCGGCATCCATCCCCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGtt                 <  1:493617/149‑1 (MQ=255)
                                                                                                                          aaataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTa                <  2:1081520/149‑1 (MQ=255)
                                                                                                                          aaataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTa                >  1:584601/1‑149 (MQ=255)
                                                                                                                            ataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGc              <  1:1226255/149‑1 (MQ=255)
                                                                                                                                aatCTGGGTGTTCGACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTCCCCCTGGCGCAGCGTTTTGTTAGTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCg          <  2:210350/149‑1 (MQ=255)
                                                                                                                                    tGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAGCAGCAAATTGGCGCGGAAGGTGGGTTTCGGGTTGCCGTCggagc      >  2:991513/1‑146 (MQ=255)
                                                                                                                                    tGGGATTCCAACGGCGGCATCCATCTACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGa      <  1:898131/149‑1 (MQ=255)
                                                                                                                                        aTTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATaaa  <  2:357213/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
ACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAA  >  NC_000913/1706826‑1707110

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑