Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,707,029 A→G V37V (GTA→GTG)  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,707,0290AG99.3% 475.5 / ‑4.8 141V37V (GTA→GTG) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (84/55);  total (84/57)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 141 total aligned reads displayed.

CTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAG  >  NC_000913/1706933‑1707125
                                                                                                |                                                                                                
ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGccc                                                                                               >  1:1928974/1‑100 (MQ=255)
ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGccc                                                                                               >  1:1909694/1‑100 (MQ=255)
 tctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTACAACGGTACACCCCTGCGCCAGGTGcccc                                                                                              >  1:2098525/1‑100 (MQ=255)
  ctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGcccc                                                                                              <  1:3382985/99‑1 (MQ=255)
   tGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGcccc                                                                                              <  1:505079/98‑1 (MQ=255)
      aTTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGGCAGGTGCCCCTGgcg                                                                                         >  1:1084979/1‑100 (MQ=255)
      aTTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGgcg                                                                                         <  2:1356411/100‑1 (MQ=255)
       ttCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGCACCTGgcgc                                                                                        <  2:1118615/100‑1 (MQ=255)
        tcgaaaaaaaTAAAAGCTGGGATTTCCACGGCGCCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGCACCTGgcgc                                                                                        <  1:2977295/95‑1 (MQ=255)
        tCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGgcgc                                                                                        >  1:1495628/1‑99 (MQ=255)
        tCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCa                                                                                       >  2:2032261/1‑100 (MQ=255)
        tCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCa                                                                                       <  1:405427/100‑1 (MQ=255)
         cAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCa                                                                                       <  1:2867499/99‑1 (MQ=255)
         cAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAg                                                                                      <  2:1576861/100‑1 (MQ=255)
          aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCGAACGGTACACCCCTGCGCCAGGGGCACCTGGCGCAgg                                                                                     <  1:915695/100‑2 (MQ=255)
          aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGCACCTGGCGCAgg                                                                                     <  2:1905368/100‑2 (MQ=255)
          aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGCACCTGGCGCAgg                                                                                     <  2:1354785/100‑2 (MQ=255)
           gAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCg                                                                                    >  2:2659541/1‑100 (MQ=255)
           gAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCg                                                                                    >  1:2921918/1‑100 (MQ=255)
           gAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCg                                                                                    >  2:608837/1‑100 (MQ=255)
           gAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCg                                                                                    >  2:1180251/1‑100 (MQ=255)
            aaaaaaTAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGt                                                                                   <  2:1944578/100‑1 (MQ=255)
            aaaaaaTAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGt                                                                                   >  2:2640644/1‑100 (MQ=255)
                aataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTg                                                                               <  1:3423517/100‑1 (MQ=255)
                aataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTg                                                                               <  2:1231686/100‑1 (MQ=255)
                aataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTg                                                                               <  2:2713132/100‑1 (MQ=255)
                aataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTg                                                                               <  1:2954599/100‑1 (MQ=255)
                aataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTg                                                                               <  2:86852/100‑1 (MQ=255)
                  taaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGtt                                                                             <  2:3198434/100‑1 (MQ=255)
                       tCTGGGATTTCAATGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTcc                                                                        <  1:1130985/100‑1 (MQ=255)
                       cttccgATCTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTcc                                                                        <  2:2311579/95‑1 (MQ=255)
                          gggATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACt                                                                     >  1:2672789/1‑100 (MQ=255)
                           ggATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTg                                                                    <  1:229841/100‑1 (MQ=255)
                            gATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGa                                                                   >  2:3047341/1‑100 (MQ=255)
                             aTTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGaa                                                                  <  2:1975497/100‑1 (MQ=255)
                             aTTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGaa                                                                  <  1:3306465/100‑1 (MQ=255)
                             aTTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGa                                                                   >  1:2430872/1‑99 (MQ=255)
                               ttCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAAc                                                                <  2:2168538/100‑1 (MQ=255)
                                tCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACa                                                               >  2:756107/1‑100 (MQ=255)
                                 cAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCCGCGTTTTGTTATTCCagatcggaag                                                              >  1:2311547/1‑91 (MQ=255)
                                 cAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAg                                                              >  2:2129629/1‑100 (MQ=255)
                                 cAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAg                                                              >  1:3242314/1‑100 (MQ=255)
                                        ggcATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGGCAGGTGCCCCTGGCGCAGCGTTTCGTTATTCCACTGCACCAGCCTAGTg                                                       >  1:373523/1‑100 (MQ=255)
                                          catccatcCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGgc                                                     >  2:1653946/1‑100 (MQ=255)
                                             ccatccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGATTTGTTATTCCACTGAAACAGCATATTGGCGCt                                                  >  1:1948335/1‑100 (MQ=255)
                                              catccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCGGGTGACCCGGGCGCAGCGCCGTGTTATCCCACTGAAACAGCGTATCGgcg                                                    >  1:1370600/1‑97 (MQ=255)
                                               atccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGa                                                >  2:647708/1‑100 (MQ=255)
                                                 ccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAg                                              <  1:2560114/100‑1 (MQ=255)
                                                  caCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAgg                                             >  2:2798263/1‑100 (MQ=255)
                                                  caCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAgg                                             <  2:3315063/100‑1 (MQ=255)
                                                  caCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAgg                                             <  2:390372/100‑1 (MQ=255)
                                                    ccGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTg                                           >  2:697912/1‑100 (MQ=255)
                                                    ccGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTg                                           <  2:54484/100‑1 (MQ=255)
                                                      ggAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAg                                         <  2:3089104/100‑1 (MQ=255)
                                                      ggAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAg                                         >  1:10000/1‑100 (MQ=255)
                                                        agaTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGtt                                       <  1:1072224/100‑1 (MQ=255)
                                                        agaTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGtt                                       >  1:3247828/1‑100 (MQ=255)
                                                         gaTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtg                                      >  2:3003591/1‑100 (MQ=255)
                                                          aTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgt                                     >  1:2084422/1‑100 (MQ=255)
                                                           tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                                    >  1:1329169/1‑100 (MQ=255)
                                                           tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                                    <  1:674270/100‑1 (MQ=255)
                                                           tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                                    >  2:2523411/1‑100 (MQ=255)
                                                            gAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGc                                   >  2:349345/1‑100 (MQ=255)
                                                            gAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGc                                   >  1:2882307/1‑100 (MQ=255)
                                                             aaaaCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCg                                  <  2:2875394/100‑1 (MQ=255)
                                                              aaaCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGt                                 >  1:1103566/1‑100 (MQ=255)
                                                                 cccAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAg                              <  2:67258/100‑1 (MQ=255)
                                                                 cccAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAg                              <  2:2142150/100‑1 (MQ=255)
                                                                 cccAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAg                              >  2:3410189/1‑100 (MQ=255)
                                                                 cccAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAg                              >  2:1249162/1‑100 (MQ=255)
                                                                  ccAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCACATTGGCGCTGAAGGTGAGTTGTGCGTTAGc                             >  1:2319640/1‑100 (MQ=255)
                                                                    aGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTTTTATTCCACTTAAACAGCATATTTGCTCTGAAGGTGAGTTGGGCGTTAGCGt                           >  2:158196/1‑100 (MQ=255)
                                                                       ccAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCgg                        >  2:395290/1‑100 (MQ=255)
                                                                       ccAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCgg                        >  2:2629632/1‑100 (MQ=255)
                                                                        cAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGc                       >  2:3348032/1‑100 (MQ=255)
                                                                         aaCGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGc                       <  2:1313408/99‑1 (MQ=255)
                                                                          aCGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGa                     >  1:2384759/1‑100 (MQ=255)
                                                                             acaaaacgcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATaa                  >  2:1588975/9‑100 (MQ=38)
                                                                              tACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATaaa                 <  2:2168321/100‑1 (MQ=255)
                                                                              tACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATaaa                 >  2:2189564/1‑100 (MQ=255)
                                                                               aaaacgcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAg                >  1:3419444/7‑100 (MQ=38)
                                                                                  ccccTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGTCGATCAAGTa              >  1:406171/1‑99 (MQ=255)
                                                                                  ccccTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAt             >  2:941547/1‑100 (MQ=255)
                                                                                  ccccTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAt             >  2:229006/1‑100 (MQ=255)
                                                                                  ccccTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATACTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGt               <  2:2921369/98‑1 (MQ=255)
                                                                                   ngcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAtt            >  1:1249237/3‑100 (MQ=255)
                                                                                   cgcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGAGAAAGTAtt            >  1:656646/3‑100 (MQ=255)
                                                                                    gcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTg           >  2:2693101/2‑100 (MQ=255)
                                                                                    gcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTg           >  2:2835527/2‑100 (MQ=255)
                                                                                    ccTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTg           >  2:2712311/1‑100 (MQ=255)
                                                                                    ccTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTg           >  2:2712107/1‑100 (MQ=255)
                                                                                     cTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTgc          >  2:2423263/1‑100 (MQ=255)
                                                                                     cTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTg           >  2:2704315/1‑99 (MQ=255)
                                                                                         gcaaGGTGCCCCTGGCGCGGCGTTTTGTGATTCCACTGAAACAGCATATGGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCgg      <  1:1632571/97‑1 (MQ=255)
                                                                                         gcCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCgg      >  1:3234094/1‑100 (MQ=255)
                                                                                          ccAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGc     >  2:576543/1‑100 (MQ=255)
                                                                                          ccAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGc     >  2:3201031/1‑100 (MQ=255)
                                                                                           cAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGcc    <  1:978387/100‑1 (MQ=255)
                                                                                            aGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCa   >  2:2871454/1‑100 (MQ=255)
                                                                                             ggTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAg  <  2:313254/100‑1 (MQ=255)
                                                                                                |                                                                                                
CTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAG  >  NC_000913/1706933‑1707125

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