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,2560GA98.8% 872.3 / ‑3.0 245E65K (GAA→AAA) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (127/115);  total (130/115)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 272 total aligned reads displayed.

AGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTC  >  NC_000913/2621166‑2621345
                                                                                          |                                                                                         
aGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaa                                                                                         >  2:7664167/1‑93 (MQ=255)
  gATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaaga                                                                                       <  1:6739410/93‑1 (MQ=255)
  gATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaaga                                                                                       >  1:5317819/1‑93 (MQ=255)
  gATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaaga                                                                                       >  2:4264075/1‑93 (MQ=255)
  gATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaaga                                                                                       <  1:1896886/93‑1 (MQ=255)
   aTGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaagaa                                                                                      <  1:7082072/93‑1 (MQ=255)
    tGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAAc                                                                                     <  2:8751160/93‑1 (MQ=255)
    tGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAAc                                                                                     <  2:3683546/93‑1 (MQ=255)
     gCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACa                                                                                    <  1:4992413/93‑1 (MQ=255)
     gCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACa                                                                                    <  1:4622958/93‑1 (MQ=255)
      cGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACaa                                                                                   <  1:1470016/93‑1 (MQ=255)
       gTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAg                                                                                  <  1:8588505/93‑1 (MQ=255)
       gTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAg                                                                                  <  2:5724019/93‑1 (MQ=255)
       gTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAg                                                                                  <  1:4368270/93‑1 (MQ=255)
       gTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAg                                                                                  >  2:196701/1‑93 (MQ=255)
       gTTTCCTGGGAATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAg                                                                                  >  1:5275684/1‑93 (MQ=255)
          tCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCt                                                                               <  2:6271360/93‑1 (MQ=255)
              ggggATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCaa                                                                           >  1:7264280/1‑93 (MQ=255)
               gggATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAAc                                                                          >  2:1555716/1‑93 (MQ=255)
                 gATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaa                                                                                         <  1:3649561/76‑1 (MQ=255)
                 gATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAActc                                                                        <  1:8689022/93‑1 (MQ=255)
                   tCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACtctc                                                                      <  1:5722843/93‑1 (MQ=255)
                   tCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACtctc                                                                      >  1:6096852/1‑93 (MQ=255)
                     tATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCAt                                                                    >  1:6733432/1‑93 (MQ=255)
                     tATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCAt                                                                    <  2:2371226/93‑1 (MQ=255)
                     tATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCAt                                                                    >  2:8387673/1‑93 (MQ=255)
                     tATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCAt                                                                    >  2:7398079/1‑93 (MQ=255)
                     tATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCAt                                                                    <  2:4968975/93‑1 (MQ=255)
                      aTTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCAtt                                                                   <  2:3377631/93‑1 (MQ=255)
                       ttCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTc                                                                  >  2:6656261/1‑93 (MQ=255)
                        tCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTcc                                                                 <  1:8329779/93‑1 (MQ=255)
                         ccAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTccc                                                                <  2:7973761/93‑1 (MQ=255)
                          cAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCg                                                               <  2:4526034/93‑1 (MQ=255)
                          cAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCg                                                               <  2:4722278/93‑1 (MQ=255)
                          cAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCg                                                               >  1:2980720/1‑93 (MQ=255)
                           aaTAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGTTCCAACTCTCATTCCCGc                                                              <  2:8625605/93‑1 (MQ=255)
                           aaTAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGc                                                              >  2:241413/1‑93 (MQ=255)
                              aaCCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATGAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAt                                                           <  2:7181568/93‑1 (MQ=255)
                               aCCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAtt                                                          >  1:600434/1‑93 (MQ=255)
                                  ttGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTAc                                                       >  2:5416422/1‑93 (MQ=255)
                                  ttGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTAc                                                       >  1:1080672/1‑93 (MQ=255)
                                  ttGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTAc                                                       >  1:4496476/1‑93 (MQ=255)
                                  ttGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTAATTCCCGCCATTTAc                                                       <  1:736600/93‑1 (MQ=255)
                                   tgatgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACt                                                      >  1:482725/1‑93 (MQ=255)
                                    gatgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTg                                                     <  1:6166500/93‑1 (MQ=255)
                                    gatgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTg                                                     >  2:2027708/1‑93 (MQ=255)
                                    gatgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTg                                                     >  2:61370/1‑93 (MQ=255)
                                     attaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCACCTCTCATTCCCGCCATTTACTgg                                                    <  1:5474369/90‑1 (MQ=255)
                                     atgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTgg                                                    <  2:1602207/93‑1 (MQ=255)
                                     atgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTgg                                                    >  1:4112527/1‑93 (MQ=255)
                                        aGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCa                                                 <  2:2848369/93‑1 (MQ=255)
                                        aGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCa                                                 <  2:5052203/93‑1 (MQ=255)
                                        aGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAAATCTCATTCCCGCCATTTACTGGGCa                                                 >  2:5989850/1‑93 (MQ=255)
                                         gTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCaa                                                <  2:217808/93‑1 (MQ=255)
                                          ttCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCaaa                                               >  2:3834376/1‑93 (MQ=255)
                                           tCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCaaaa                                              >  1:3394854/1‑93 (MQ=255)
                                             tataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAAtt                                            >  1:966575/1‑93 (MQ=255)
                                             tataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAAtt                                            >  2:3129599/1‑93 (MQ=255)
                                             tataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAAtt                                            >  1:801420/1‑93 (MQ=255)
                                             tataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAAtt                                            <  1:7730795/93‑1 (MQ=255)
                                             tataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAAtt                                            >  1:3252486/1‑93 (MQ=255)
                                               taAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCa                                          >  2:1088420/1‑93 (MQ=255)
                                                 aaGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGt                                        <  1:1754063/93‑1 (MQ=255)
                                                  aGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTc                                       >  1:6690623/1‑93 (MQ=255)
                                                     ccGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCg                                    >  1:7886233/1‑93 (MQ=255)
                                                     ccGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCg                                    >  2:7328082/1‑93 (MQ=255)
                                                      cGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCgg                                   >  1:8163455/1‑93 (MQ=255)
                                                       gCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCggg                                  >  2:2541289/1‑93 (MQ=255)
                                                        cTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGt                                 >  1:6864505/1‑93 (MQ=255)
                                                         ttCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAACATTCAGTCCCGGGTg                                >  2:2323902/1‑93 (MQ=255)
                                                         ttCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTg                                >  1:1004062/1‑93 (MQ=255)
                                                          tCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGc                               <  1:932880/93‑1 (MQ=255)
                                                             ctgaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGa                            >  1:6266826/1‑93 (MQ=255)
                                                              tgaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGaa                           >  2:8169514/1‑93 (MQ=255)
                                                               gaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGaaa                          <  1:4036729/93‑1 (MQ=255)
                                                                aactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAg                         >  1:7533928/1‑93 (MQ=255)
                                                                 actgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGc                        <  1:4710688/93‑1 (MQ=255)
                                                                 actgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGc                        >  1:2041567/1‑93 (MQ=255)
                                                                 actgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGc                        <  1:2326499/93‑1 (MQ=255)
                                                                  ctgaaGCGACGCATCATTATTAGCGAAGAACAAGGGTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCACGGGTGCTGAAAGcc                       >  1:5557062/1‑93 (MQ=255)
                                                                     aaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGAt                    >  2:8524321/1‑93 (MQ=255)
                                                                     aaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGAt                    >  1:3161605/1‑93 (MQ=255)
                                                                        cGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAg                 <  2:5947731/93‑1 (MQ=255)
                                                                         gACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAgg                <  1:7666261/93‑1 (MQ=255)
                                                                          aCGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGa               <  1:1245216/93‑1 (MQ=255)
                                                                          aCGCATAATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGa               <  1:5146945/93‑1 (MQ=255)
                                                                            gCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAAt             >  2:5660727/1‑93 (MQ=255)
                                                                            gCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAAt             >  2:1197551/1‑93 (MQ=255)
                                                                            gCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGACATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAAt             >  1:1424932/1‑93 (MQ=255)
                                                                             catcatTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAAtt            >  2:7023069/1‑93 (MQ=255)
                                                                              atcatTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTc           >  1:4210852/1‑93 (MQ=255)
                                                                              atcatTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTc           <  1:688184/93‑1 (MQ=255)
                                                                               tcatTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCg          >  1:2805806/1‑93 (MQ=255)
                                                                               tcatTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCg          <  2:4020999/93‑1 (MQ=255)
                                                                               tcatTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCg          >  1:6118514/1‑93 (MQ=255)
                                                                                  ttattaGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACg       <  1:5491620/93‑1 (MQ=255)
                                                                                  ttattaGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACg       >  1:3770588/1‑93 (MQ=255)
                                                                                   tattaGGAAAGAAAAAGGGTCCAATGTTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACgg      <  1:5988702/93‑1 (MQ=255)
                                                                                   tattaGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACgg      >  1:6592831/1‑93 (MQ=255)
                                                                                      taGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCAGGGTGCTGAAAGCCGATCAGGAATTCGACGGCct   <  1:7012699/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTC  >  NC_000913/2621166‑2621345

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