Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,621,256 G→A E65K (GAA→AAA)  ppk → polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,621,2560GA100.0% 836.1 / NA 211E65K (GAA→AAA) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (119/92);  total (119/92)

Only 100 of 221 total aligned reads displayed.

AGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTC  >  NC_000913/2621166‑2621345
                                                                                          |                                                                                         
aGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaa                                                                                         <  2:307693/93‑1 (MQ=255)
 ggATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaag                                                                                        >  1:4234861/1‑93 (MQ=255)
  gATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaaga                                                                                       >  2:5847270/1‑93 (MQ=255)
    tGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAAc                                                                                     <  1:4508852/93‑1 (MQ=255)
     gCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACa                                                                                    >  1:278508/1‑93 (MQ=255)
      cGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACaa                                                                                   >  1:5187445/1‑93 (MQ=255)
      cGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCAGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACaa                                                                                   >  2:2927834/1‑93 (MQ=255)
         ttCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGc                                                                                >  1:1246153/1‑93 (MQ=255)
         ttCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGc                                                                                >  1:3362270/1‑93 (MQ=255)
           ccTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTc                                                                              <  2:3870915/93‑1 (MQ=255)
           ccTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTc                                                                              <  2:3124611/93‑1 (MQ=255)
           ccTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTc                                                                              <  1:2092964/93‑1 (MQ=255)
            cTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTcc                                                                             >  1:2857513/1‑93 (MQ=255)
            cTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTcc                                                                             >  2:3614386/1‑93 (MQ=255)
            cTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTcc                                                                             <  1:5370274/93‑1 (MQ=255)
            cTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTcc                                                                             >  2:2140937/1‑93 (MQ=255)
              ggggATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCaa                                                                           <  2:1835623/93‑1 (MQ=255)
              ggggATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCaa                                                                           >  1:3607669/1‑93 (MQ=255)
                 gATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAActc                                                                        <  2:2394976/93‑1 (MQ=255)
                   tCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACtctc                                                                      <  2:5335525/93‑1 (MQ=255)
                   tCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACtctc                                                                      <  1:4340244/93‑1 (MQ=255)
                   tCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACtctc                                                                      <  1:35348/93‑1 (MQ=255)
                    cTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCa                                                                     >  2:4324586/1‑93 (MQ=255)
                     tATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCAt                                                                    <  2:4756025/93‑1 (MQ=255)
                      aTTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCAtt                                                                   <  1:5992775/93‑1 (MQ=255)
                       ttCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTc                                                                  <  2:4135939/93‑1 (MQ=255)
                       ttCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTc                                                                  >  1:613258/1‑93 (MQ=255)
                        tCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTcc                                                                 <  2:6278606/93‑1 (MQ=255)
                         ccAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTccc                                                                >  1:3041202/1‑93 (MQ=255)
                         ccAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTccc                                                                <  1:4380859/93‑1 (MQ=255)
                          cAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCg                                                               >  1:5258785/1‑93 (MQ=255)
                          cAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCg                                                               >  1:4553884/1‑93 (MQ=255)
                            aTAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGcc                                                             <  2:386085/93‑1 (MQ=255)
                              aaCCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAt                                                           >  1:2519532/1‑93 (MQ=255)
                              aaCCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAt                                                           <  1:4139147/93‑1 (MQ=255)
                               aCCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAtt                                                          <  2:1965955/93‑1 (MQ=255)
                               aCCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAtt                                                          >  2:2992408/1‑93 (MQ=255)
                               aCCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAtt                                                          >  1:972511/1‑93 (MQ=255)
                               aCCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAtt                                                          >  1:5127418/1‑93 (MQ=255)
                                ccTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAttt                                                         >  2:5828489/1‑93 (MQ=255)
                                 cTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTa                                                        <  1:4757518/93‑1 (MQ=255)
                                  ttGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTAc                                                       >  1:3377676/1‑93 (MQ=255)
                                   tgatgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACt                                                      <  2:1100305/93‑1 (MQ=255)
                                    gatgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTg                                                     <  2:5470175/93‑1 (MQ=255)
                                       gaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGc                                                  >  2:144114/1‑93 (MQ=255)
                                        aGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCa                                                 >  2:3320874/1‑93 (MQ=255)
                                          ttCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCaaa                                               <  1:2678713/93‑1 (MQ=255)
                                            cTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAAt                                             >  2:4427240/1‑93 (MQ=255)
                                             tataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAAtt                                            >  2:4817800/1‑93 (MQ=255)
                                             tataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAAtt                                            <  1:105239/93‑1 (MQ=255)
                                              ataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCGACCATTTACTAGGAAAAttta                                           >  1:3139822/1‑92 (MQ=255)
                                              ataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTc                                           >  1:3581180/1‑93 (MQ=255)
                                              ataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTc                                           >  1:2882859/1‑93 (MQ=255)
                                              ataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTc                                           <  2:3300740/93‑1 (MQ=255)
                                               taAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCa                                          <  2:863290/93‑1 (MQ=255)
                                                aaaGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAg                                         <  2:1307430/93‑1 (MQ=255)
                                                  aGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTc                                       >  2:4066024/1‑93 (MQ=255)
                                                   gTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTcc                                      <  2:2193947/93‑1 (MQ=255)
                                                      cGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCTCGCCATTTACTGGGCAAAATTCAGTCCCgg                                   >  1:564752/1‑93 (MQ=255)
                                                      cGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCgg                                   >  2:2082501/1‑93 (MQ=255)
                                                       gCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTACAACTCTCATTCCCGACATTTACTGGGCAAAATTCAGTCCCggg                                  >  2:2831206/1‑93 (MQ=255)
                                                         ttCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTg                                >  2:763882/1‑93 (MQ=255)
                                                           cGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCt                              <  2:3419429/93‑1 (MQ=255)
                                                           cGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCt                              <  2:1913224/93‑1 (MQ=255)
                                                           cGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCt                              >  1:2905908/1‑93 (MQ=255)
                                                           cGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCt                              <  1:4748749/93‑1 (MQ=255)
                                                            gCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTg                             >  2:4934754/1‑93 (MQ=255)
                                                             ctgaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGa                            <  1:4492971/93‑1 (MQ=255)
                                                              tgaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGaa                           <  1:1547711/93‑1 (MQ=255)
                                                               gaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGaaa                          >  2:5231815/1‑93 (MQ=255)
                                                               gaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGaaa                          >  2:1738012/1‑93 (MQ=255)
                                                                aactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAg                         >  1:1439867/1‑93 (MQ=255)
                                                                aactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAg                         >  2:1383888/1‑93 (MQ=255)
                                                                 actgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGc                        >  1:4455508/1‑93 (MQ=255)
                                                                 actgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGc                        >  2:6172845/1‑93 (MQ=255)
                                                                  ctgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGcc                       <  1:5595212/93‑1 (MQ=255)
                                                                    gaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGa                     >  2:5423045/1‑93 (MQ=255)
                                                                     aaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGAt                    <  1:4260386/93‑1 (MQ=255)
                                                                     aaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGAt                    >  1:214870/1‑93 (MQ=255)
                                                                     aaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGAt                    <  1:5798402/93‑1 (MQ=255)
                                                                      aGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATc                   <  2:2651048/93‑1 (MQ=255)
                                                                       gCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCa                  <  1:134828/93‑1 (MQ=255)
                                                                        cGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAg                 <  2:1588432/93‑1 (MQ=255)
                                                                        cGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAg                 <  1:5418705/93‑1 (MQ=255)
                                                                        cGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAg                 >  1:6074522/1‑93 (MQ=255)
                                                                         gACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAgg                >  2:4153395/1‑93 (MQ=255)
                                                                         gACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAgg                <  1:3853575/93‑1 (MQ=255)
                                                                           cGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGaa              <  1:363867/93‑1 (MQ=255)
                                                                           cGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGaa              >  1:3426121/1‑93 (MQ=255)
                                                                            gCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAAt             >  1:870346/1‑93 (MQ=255)
                                                                               tcatTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCg          >  1:2326256/1‑93 (MQ=255)
                                                                                catTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGa         >  2:2383435/1‑93 (MQ=255)
                                                                                 attattAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGAc        >  1:977340/1‑93 (MQ=255)
                                                                                 attattAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGAc        >  2:3848468/1‑93 (MQ=255)
                                                                                    attaGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGc     >  2:2968827/1‑93 (MQ=255)
                                                                                     ttaGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGcc    <  1:5888557/93‑1 (MQ=255)
                                                                                     ttaGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGcc    <  2:6433164/93‑1 (MQ=255)
                                                                                     ttaGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGcc    <  1:31441/93‑1 (MQ=255)
                                                                                       aGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCctc  >  1:725116/1‑93 (MQ=255)
                                                                                       aGCAAAGAACAAGGCTCAAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCctc  <  2:1997228/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTC  >  NC_000913/2621166‑2621345

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