Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 669,863 T→C A322A (GCA→GCG holA ← DNA polymerase III delta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913669,8630TC95.6% 463.1 / 7.3 159A322A (GCA→GCGholADNA polymerase III delta subunit
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base C (92/60);  total (97/62)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-01

Only 100 of 163 total aligned reads displayed.

AACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAA  >  NC_000913/669773‑669952
                                                                                          |                                                                                          
aaCAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGc                                                                                          <  2:1678389/93‑1 (MQ=255)
aaCAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGACAGGGGGTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGc                                                                                          >  1:278698/1‑93 (MQ=255)
   agagCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCca                                                                                       <  1:3096340/93‑1 (MQ=255)
   agagCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCca                                                                                       <  2:1176131/93‑1 (MQ=255)
   agagCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCca                                                                                       >  1:672266/1‑93 (MQ=255)
   agagCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCca                                                                                       <  1:1594411/93‑1 (MQ=255)
    gagCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCcac                                                                                      >  2:4096419/1‑93 (MQ=255)
     agCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCcaca                                                                                     <  2:2564127/93‑1 (MQ=255)
     agCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCcaca                                                                                     <  1:2858778/93‑1 (MQ=255)
     agCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCcaca                                                                                     <  1:4498592/93‑1 (MQ=255)
     agCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCcaca                                                                                     >  2:4277439/1‑93 (MQ=255)
      gCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCacac                                                                                    <  2:514907/93‑1 (MQ=255)
      gCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCacac                                                                                    <  2:1914642/93‑1 (MQ=255)
       ccTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACact                                                                                   >  1:1720345/1‑93 (MQ=255)
        cTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACactg                                                                                  <  1:1424377/93‑1 (MQ=255)
        cTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACactg                                                                                  <  2:2661623/93‑1 (MQ=255)
        cTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGCTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACactg                                                                                  >  1:894280/1‑93 (MQ=255)
         tGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACactga                                                                                 >  1:2914987/1‑93 (MQ=255)
         tGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACactga                                                                                 <  1:2435673/93‑1 (MQ=255)
           tAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACactgact                                                                               >  1:163672/1‑93 (MQ=255)
            aaaGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACactgactg                                                                              >  2:2073804/1‑93 (MQ=255)
              aGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACtgactgac                                                                            <  2:4365275/93‑1 (MQ=255)
                aTTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCg                                                                          <  1:3625738/93‑1 (MQ=255)
                   tCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTaa                                                                       >  1:1681337/1‑93 (MQ=255)
                   tCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGCaa                                                                       >  1:3146937/1‑93 (MQ=255)
                     atatCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATc                                                                     >  2:1248843/1‑93 (MQ=255)
                      tatCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCt                                                                   <  1:2359085/93‑1 (MQ=255)
                       atCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCtt                                                                  >  1:1972604/1‑93 (MQ=255)
                       atCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCtt                                                                  <  1:544430/93‑1 (MQ=255)
                       atCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCtt                                                                  <  1:1070233/93‑1 (MQ=255)
                             cGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGa                                                            >  2:1482715/1‑93 (MQ=255)
                             cGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGa                                                            <  2:2448778/93‑1 (MQ=255)
                               tCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAgg                                                          >  2:1928010/1‑93 (MQ=255)
                               tCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAgg                                                          >  2:867689/1‑93 (MQ=255)
                                cGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAggg                                                         <  2:3435830/93‑1 (MQ=255)
                                cGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAggg                                                         <  1:2871964/93‑1 (MQ=255)
                                cGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAggg                                                         <  2:1038192/93‑1 (MQ=255)
                                cGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAggg                                                         <  2:2866398/93‑1 (MQ=255)
                                cGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTACCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAggg                                                         <  1:2510295/93‑1 (MQ=255)
                                cGATAAATACGTCCGACAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCTCACTGACTGACCGTAATCTTGTTTGAggg                                                         <  2:4687554/93‑1 (MQ=255)
                                 gATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGt                                                        >  2:694768/1‑93 (MQ=255)
                                  aTAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATACCCCTCCAACTCCGCCCCACCTGACTGACCCGAATCTTTGTTTAAggtgg                                                       >  1:2207725/1‑90 (MQ=255)
                                   tAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGa                                                      >  1:3503463/1‑93 (MQ=255)
                                   tAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGa                                                      <  2:673210/93‑1 (MQ=255)
                                     aaTACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGt                                                    >  1:740088/1‑93 (MQ=255)
                                     aaTACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGt                                                    >  2:2762/1‑93 (MQ=255)
                                     aaTACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGt                                                    >  2:2495742/1‑93 (MQ=255)
                                     aaTACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGt                                                    >  2:1405125/1‑93 (MQ=255)
                                      aTACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGtt                                                   <  2:3174829/93‑1 (MQ=255)
                                      aTACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGtt                                                   >  2:4236633/1‑93 (MQ=255)
                                       tACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTc                                                  >  1:1584229/1‑93 (MQ=255)
                                       tACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTc                                                  >  2:1712314/1‑93 (MQ=255)
                                        aCGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTcc                                                 >  2:3146524/1‑93 (MQ=255)
                                        aCGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTcc                                                 <  1:1005500/93‑1 (MQ=255)
                                         cggccGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCg                                                <  2:1325381/90‑1 (MQ=255)
                                         cGTCCGCCAGGGGTTTATGGCACAACCGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCg                                                >  2:2737603/1‑93 (MQ=255)
                                          gTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGt                                               >  2:865269/1‑93 (MQ=255)
                                           tCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGtt                                              >  1:4559871/1‑93 (MQ=255)
                                           tCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGtt                                              <  1:1050156/93‑1 (MQ=255)
                                                cAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCgtgt                                         <  1:4372919/93‑1 (MQ=255)
                                                 aGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTc                                        >  1:1384068/1‑93 (MQ=255)
                                                 aGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTc                                        <  2:1223664/93‑1 (MQ=255)
                                                  ggggTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCa                                       >  2:4142771/1‑93 (MQ=255)
                                                  ggggTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCa                                       >  2:3796269/1‑93 (MQ=255)
                                                   gggTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAg                                      >  1:2022273/1‑93 (MQ=255)
                                                   gggTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAg                                      >  2:4525369/1‑93 (MQ=255)
                                                      tttATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAg                                   >  1:1081492/1‑93 (MQ=255)
                                                      tttATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAg                                   >  1:4404469/1‑93 (MQ=255)
                                                        tATGGCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGtt                                 >  2:2464269/1‑93 (MQ=255)
                                                           ggCACAACAGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCa                              <  1:4026083/93‑1 (MQ=255)
                                                            gCACAACAGAAGAGATACCCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCAc                             <  2:3205701/93‑1 (MQ=255)
                                                               caacaGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGATTTGCACGGc                          >  2:1830819/1‑93 (MQ=255)
                                                               caacaGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGc                          >  2:3583224/1‑93 (MQ=255)
                                                               caacaGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGc                          <  2:771514/93‑1 (MQ=255)
                                                                aacaGAAGAGATATCCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTCCCGTTCGTGTCAGGAGTTGCACGGcc                         <  2:3833352/93‑1 (MQ=255)
                                                                aacaGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGcc                         >  1:2471698/1‑93 (MQ=255)
                                                                 acaGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCt                        <  2:3320819/93‑1 (MQ=255)
                                                                 acaGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCt                        <  2:562172/93‑1 (MQ=255)
                                                                  caGAAGAGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTg                       >  1:1388171/1‑93 (MQ=255)
                                                                   agaagaGATAACCCGTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGa                      >  2:2863423/1‑93 (MQ=255)
                                                                    gaagaGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGAc                     >  1:2956847/1‑93 (MQ=255)
                                                                    gaagaGATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGAc                     >  2:2134329/1‑93 (MQ=255)
                                                                      agagATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGt                   <  1:1199985/93‑1 (MQ=255)
                                                                      agagATAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGt                   <  2:3942533/93‑1 (MQ=255)
                                                                         gaTAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAAc                >  1:3353305/1‑93 (MQ=255)
                                                                         gaTAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAAc                <  2:1808031/93‑1 (MQ=255)
                                                                          aTAACCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACt               <  1:674252/93‑1 (MQ=255)
                                                                           tAACCCTTCCAGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGAGTGAGATCCGTTCGTGTCAGGAGTTGCACGGGCTGACGTAACTg              >  1:268961/1‑93 (MQ=255)
                                                                             aCCCTTCCAGCTCCGCGCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGACGTGCACGGCCTGACGTAACTGCg            >  2:3523951/1‑93 (MQ=255)
                                                                             aCCCTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCg            >  2:2948658/1‑93 (MQ=255)
                                                                              cccTTCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGt           >  2:3344700/1‑93 (MQ=255)
                                                                                 ttCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTg        >  2:2224942/1‑93 (MQ=255)
                                                                                 ttCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTg        >  1:2299254/1‑93 (MQ=255)
                                                                                  tCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGa       >  1:3022812/1‑93 (MQ=255)
                                                                                  tCCAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGa       <  2:219561/93‑1 (MQ=255)
                                                                                   ccAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGAc      <  2:2548768/93‑1 (MQ=255)
                                                                                   ccAGCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGAc      <  1:3918338/93‑1 (MQ=255)
                                                                                      gCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTa   >  1:4548061/1‑93 (MQ=255)
                                                                                      gCTCCGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGGTCCGTTCGTGTCAGGAGGTGCACGGCCTGACGTACCTGCGTCTGACTTa   >  1:353529/1‑93 (MQ=255)
                                                                                       ctccGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTaa  <  1:2655042/93‑1 (MQ=255)
                                                                                          |                                                                                          
AACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAA  >  NC_000913/669773‑669952

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