Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,917,055 T→C Q144R (CAG→CGG)  atpA ← membrane‑bound ATP synthase, F1 sector, alpha‑subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,917,0550TC99.6% 703.7 / ‑4.7 222Q144R (CAG→CGG) atpAmembrane‑bound ATP synthase, F1 sector, alpha‑subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (103/118);  total (103/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

Only 100 of 232 total aligned reads displayed.

CCTGTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCTGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTACCGTCGATT  >  NC_000913/3916965‑3917144
                                                                                          |                                                                                         
ccTGTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGa                                                                                         <  2:3277829/93‑1 (MQ=255)
ccTGTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGa                                                                                         <  2:2046775/93‑1 (MQ=255)
ccTGTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGa                                                                                         >  1:6268703/1‑93 (MQ=255)
ccTGTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGa                                                                                         <  1:4701649/93‑1 (MQ=255)
ccTGTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGa                                                                                         >  1:4326331/1‑93 (MQ=255)
 cTGTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGAt                                                                                        >  1:1523234/1‑93 (MQ=255)
  tGTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATc                                                                                       >  2:3816963/1‑93 (MQ=255)
  tGTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATc                                                                                       <  2:154143/93‑1 (MQ=255)
   gTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCt                                                                                      <  2:3239955/93‑1 (MQ=255)
   gTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCt                                                                                      <  2:3024778/93‑1 (MQ=255)
    tCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTa                                                                                     <  2:5138327/93‑1 (MQ=255)
    tCTGACGGGCACCGATGATTAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCAGCCGATCTa                                                                                     <  2:4217223/93‑1 (MQ=255)
     cTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTAc                                                                                    >  2:3356214/1‑93 (MQ=255)
     cTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTAc                                                                                    >  2:1130554/1‑93 (MQ=255)
      tGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACg                                                                                   <  2:825720/93‑1 (MQ=255)
      tGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACg                                                                                   >  1:745483/1‑93 (MQ=255)
        aCGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGa                                                                                 <  1:511775/93‑1 (MQ=255)
             cACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGAc                                                                            <  1:2150753/93‑1 (MQ=255)
                cGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGtt                                                                         >  2:4951530/1‑93 (MQ=255)
                cGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGtt                                                                         >  1:3778872/1‑93 (MQ=255)
                 gatgatCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTc                                                                        >  1:3075741/1‑93 (MQ=255)
                 gatgatCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTc                                                                        >  1:1039206/1‑93 (MQ=255)
                   tgatCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGa                                                                      >  1:3008127/1‑93 (MQ=255)
                   tgatCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGa                                                                      >  1:5414537/1‑93 (MQ=255)
                      tCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATaa                                                                   >  1:5176813/1‑93 (MQ=255)
                      tCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATaa                                                                   <  1:2509789/93‑1 (MQ=255)
                      tCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATaa                                                                   <  2:4501261/93‑1 (MQ=255)
                       cAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAAc                                                                  <  2:5732046/93‑1 (MQ=255)
                         aTTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGc                                                                >  1:3783168/1‑93 (MQ=255)
                            cACGCTCACCACGACCGATTGGGAGCATGGAGTCAACGGCTTTATACCCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCg                                                             <  1:6016972/93‑1 (MQ=255)
                             ccGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCgg                                                            <  1:2547340/92‑1 (MQ=255)
                                 tGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCg                                                        >  2:4389841/1‑93 (MQ=255)
                                  gACCACGCCCGATTGGGATCATGGAGTCAGCGCCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGa                                                       <  1:3670192/93‑1 (MQ=255)
                                  gACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGa                                                       >  1:1203911/1‑93 (MQ=255)
                                  gACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGGTCGATAACGCCCGGAGCGa                                                       >  2:3787122/1‑93 (MQ=255)
                                   aCCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGAt                                                      <  2:3161689/93‑1 (MQ=255)
                                      aCGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGc                                                   <  2:363120/93‑1 (MQ=255)
                                       cGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCt                                                  <  1:1701029/93‑1 (MQ=255)
                                        gACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCtt                                                 >  2:4051879/1‑93 (MQ=255)
                                        gACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCtt                                                 >  1:572823/1‑93 (MQ=255)
                                         aCCGATTGGGATCATGGAGTCACCGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCTATTGCTTc                                                <  2:6154381/93‑1 (MQ=255)
                                         aCCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTc                                                <  1:2102118/93‑1 (MQ=255)
                                         aCCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTc                                                >  1:1898693/1‑93 (MQ=255)
                                         aCCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTc                                                <  1:4620568/93‑1 (MQ=255)
                                         aCCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTc                                                >  2:4592296/1‑93 (MQ=255)
                                         aCCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTc                                                <  2:4820540/93‑1 (MQ=255)
                                         aCCGATGGGGAGCATGTAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTc                                                <  2:1279958/93‑1 (MQ=255)
                                          ccGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCt                                               >  2:1266973/1‑93 (MQ=255)
                                          ccGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCt                                               >  2:5606613/1‑93 (MQ=255)
                                            gATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTAc                                             <  1:3180495/93‑1 (MQ=255)
                                            gATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTAc                                             >  1:86542/1‑93 (MQ=255)
                                             aTTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTAca                                            <  2:2769607/93‑1 (MQ=255)
                                              ttGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTAcag                                           >  2:4074120/1‑93 (MQ=255)
                                              ttGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTAcag                                           >  2:2474218/1‑93 (MQ=255)
                                                gggATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTAcagca                                         <  2:4642975/93‑1 (MQ=255)
                                                gggATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTAcagca                                         <  2:2086696/93‑1 (MQ=255)
                                                gggATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTAcagca                                         <  2:2339067/93‑1 (MQ=255)
                                                  gATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTAcagccgg                                       >  2:1802123/1‑90 (MQ=255)
                                                  gATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCaga                                       <  1:4374616/93‑1 (MQ=255)
                                                  gATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCaga                                       <  1:4929376/93‑1 (MQ=255)
                                                    tCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAgaga                                     <  1:5148171/93‑1 (MQ=255)
                                                      aTGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAg                                   >  2:448530/1‑93 (MQ=255)
                                                       tGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGc                                  <  1:397692/93‑1 (MQ=255)
                                                       tGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGc                                  <  1:3906635/93‑1 (MQ=255)
                                                       tGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGc                                  <  1:1867453/93‑1 (MQ=255)
                                                       tGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGc                                  <  1:448301/93‑1 (MQ=255)
                                                        ggATTCACCGGCTTTATACCCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGcc                                 <  2:5408022/93‑1 (MQ=255)
                                                        ggAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGcc                                 <  2:2147608/93‑1 (MQ=255)
                                                        ggAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACCGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGcc                                 >  1:4432162/1‑93 (MQ=255)
                                                         gAGTCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCcg                                >  2:744748/1‑93 (MQ=255)
                                                            tCAACGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCcgtcg                             >  2:3745326/1‑93 (MQ=255)
                                                             cAACGGCTTTATAACCGGTCTGTACTGGACGATCTACGGACTGACGTTCGATAACGACCGGAGCGAGTGCTTCTACAGCAGCGAGGCcgtcgg                            >  2:2801499/1‑92 (MQ=255)
                                                              aaCGGCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTg                           >  1:1592732/1‑93 (MQ=255)
                                                                 ggCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATc                        <  1:2953868/93‑1 (MQ=255)
                                                                  gCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATcc                       >  1:1489410/1‑93 (MQ=255)
                                                                  gCTTTATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATcc                       >  1:940505/1‑93 (MQ=255)
                                                                    tttATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAg                     >  1:269041/1‑93 (MQ=255)
                                                                     ttATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGc                    >  2:5530828/1‑93 (MQ=255)
                                                                     ttATAACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGc                    <  1:410131/93‑1 (MQ=255)
                                                                     taataCCCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGc                    <  1:6187692/91‑1 (MQ=255)
                                                                       ataACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCgg                  >  2:1445747/1‑93 (MQ=255)
                                                                        taACCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGa                 >  2:6257692/1‑93 (MQ=255)
                                                                         aaCCGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGAc                >  2:5985657/1‑93 (MQ=255)
                                                                           ccGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCt              <  2:1055952/93‑1 (MQ=255)
                                                                            cGGTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCtt             <  2:6038764/93‑1 (MQ=255)
                                                                             ggTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCttt            >  1:3350652/1‑93 (MQ=255)
                                                                              gTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACtttcc           >  1:1072148/1‑91 (MQ=255)
                                                                              gTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTa           >  2:1954280/1‑93 (MQ=255)
                                                                              gTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTa           >  1:3558419/1‑93 (MQ=255)
                                                                              gTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTa           >  1:5456248/1‑93 (MQ=255)
                                                                              gTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTa           <  2:2571408/93‑1 (MQ=255)
                                                                              gTCTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTa           <  1:1016139/93‑1 (MQ=255)
                                                                                cTGTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTAcc         <  2:4724839/93‑1 (MQ=255)
                                                                                  gTACCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTACCGt       <  1:3874053/93‑1 (MQ=255)
                                                                                    cccGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTACCGTCg     <  2:5863564/92‑1 (MQ=255)
                                                                                    aCCGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTACCGTCg     >  2:4074985/1‑93 (MQ=255)
                                                                                     ccGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTACCGTCGa    <  1:5975104/93‑1 (MQ=255)
                                                                                      cGGCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTACCGTCGAt   <  1:2299317/93‑1 (MQ=255)
                                                                                       ggCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTACCGTCGAtt  <  2:569043/93‑1 (MQ=255)
                                                                                       ggCCGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTACCGTCGAtt  <  2:3503428/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCTGTCTGACGGTCACCGATGATCAATTCACGCTGACCACGACCGATTGGGATCATGGAGTCAACGGCTTTATAACCGGTCTGTACCGGCTGATCTACGGACTGACGTTCGATAACGCCCGGAGCGATTGCTTCTACAGCAGAGAAGCCGTCGTGATCCAGCGGACCTTTACCGTCGATT  >  NC_000913/3916965‑3917144

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