Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,917,345 A→G C47C (TGT→TGC 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,3450AG84.7% 716.9 / 119.7 308C47C (TGT→TGCatpAmembrane‑bound ATP synthase, F1 sector, alpha‑subunit
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base G (119/141);  total (144/164)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

Only 100 of 314 total aligned reads displayed.

GACCCATAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAG  >  NC_000913/3917254‑3917436
                                                                                           |                                                                                           
 aCCCATAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCtgcatgca                                                                                           <  1:3415827/93‑1 (MQ=255)
  cccATAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCaa                                                                                          <  2:3129779/93‑1 (MQ=255)
  cccATAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCaa                                                                                          >  2:6904230/1‑93 (MQ=255)
  cccATAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCaa                                                                                          <  1:3629437/93‑1 (MQ=255)
   ccATAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAAt                                                                                         <  2:1995681/93‑1 (MQ=255)
    cATAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATc                                                                                        <  1:7071704/93‑1 (MQ=255)
    cATAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATc                                                                                        <  2:6746804/93‑1 (MQ=255)
    cATAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATc                                                                                        <  2:6164868/93‑1 (MQ=255)
      tAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCgg                                                                                      <  2:2603948/93‑1 (MQ=255)
       aacaacCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGc                                                                                     <  1:3181924/93‑1 (MQ=255)
         caacCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCa                                                                                   <  2:1385351/93‑1 (MQ=255)
         caacCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCa                                                                                   >  1:5065034/1‑93 (MQ=255)
         caacCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCa                                                                                   >  1:4841866/1‑93 (MQ=255)
         caacCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCa                                                                                   <  2:4365833/93‑1 (MQ=255)
            ccGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGc                                                                                <  2:7515866/93‑1 (MQ=255)
              gCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCg                                                                              >  2:5359618/1‑93 (MQ=255)
               cACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGt                                                                             >  1:1577592/1‑93 (MQ=255)
               cACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGt                                                                             <  2:3314008/93‑1 (MQ=255)
                 ccTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGa                                                                           <  2:4009937/93‑1 (MQ=255)
                  cTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGaa                                                                          <  1:3115808/93‑1 (MQ=255)
                  cTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGaa                                                                          >  1:836258/1‑93 (MQ=255)
                  cTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGaa                                                                          <  1:7200891/93‑1 (MQ=255)
                   tACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAAt                                                                         <  1:5754855/93‑1 (MQ=255)
                    aCAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATg                                                                        <  1:1812329/93‑1 (MQ=255)
                     cAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGc                                                                       <  2:1516479/93‑1 (MQ=255)
                      agagTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCg                                                                      <  2:3925513/93‑1 (MQ=255)
                       gagTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCgg                                                                     <  1:3233796/93‑1 (MQ=255)
                        agTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCTTGCATACAATCGGCCAGGCCGTGAATGCGGa                                                                    >  1:6944192/1‑93 (MQ=255)
                        agTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGa                                                                    >  2:7342308/1‑93 (MQ=255)
                             cgcTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAaca                                                               <  1:1287809/93‑1 (MQ=255)
                             cgcTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCCGCATGCAATCGGCCAGGCCGTGAATGCGGATAaca                                                               <  2:4525124/93‑1 (MQ=255)
                              gcTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAacac                                                              <  1:2158979/93‑1 (MQ=255)
                              gcTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAacac                                                              >  1:5351062/1‑93 (MQ=255)
                              gcTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAacac                                                              <  1:5426385/93‑1 (MQ=255)
                               cTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACAcc                                                             <  1:3811080/93‑1 (MQ=255)
                                tCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCg                                                            >  1:2042451/1‑93 (MQ=255)
                                 cGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGt                                                           <  1:2918095/93‑1 (MQ=255)
                                  gAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTc                                                          <  1:102270/93‑1 (MQ=255)
                                   aGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCa                                                         >  1:529613/1‑93 (MQ=255)
                                    ggTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCAc                                                        >  1:4449158/1‑93 (MQ=255)
                                    ggTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCAc                                                        <  1:6898668/93‑1 (MQ=255)
                                     gTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACt                                                       <  2:2708191/93‑1 (MQ=255)
                                     gTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACt                                                       <  1:1208425/93‑1 (MQ=255)
                                      ttCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACtt                                                      >  1:6602599/1‑93 (MQ=255)
                                      ttCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACtt                                                      >  2:2569860/1‑93 (MQ=255)
                                       tCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTa                                                     <  1:6313604/93‑1 (MQ=255)
                                       tCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTa                                                     >  2:7302280/1‑93 (MQ=255)
                                       tCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTa                                                     >  2:7106813/1‑93 (MQ=255)
                                          gTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAg                                                  >  1:6416518/1‑93 (MQ=255)
                                          gTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAg                                                  <  2:1207359/93‑1 (MQ=255)
                                             cGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGaaa                                               <  2:3218023/93‑1 (MQ=255)
                                             cGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGaaa                                               <  2:3442358/93‑1 (MQ=255)
                                              gATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAAc                                              >  2:5284598/1‑93 (MQ=255)
                                               aTAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACa                                             <  1:2794198/93‑1 (MQ=255)
                                                 aGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAAt                                           <  2:2861714/93‑1 (MQ=255)
                                                 aGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAAt                                           <  2:483270/93‑1 (MQ=255)
                                                  gCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATa                                          <  2:4182095/93‑1 (MQ=255)
                                                  gCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATa                                          >  2:5946671/1‑93 (MQ=255)
                                                   cGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAg                                         <  2:348347/93‑1 (MQ=255)
                                                    gTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGt                                        <  1:4140190/93‑1 (MQ=255)
                                                     tAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTa                                       >  2:1697522/1‑93 (MQ=255)
                                                     tAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTa                                       >  2:5054235/1‑93 (MQ=255)
                                                     tAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTa                                       >  1:3519153/1‑93 (MQ=255)
                                                      aaCGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTAc                                      <  2:6480614/93‑1 (MQ=255)
                                                       aCGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTAcc                                     <  2:5177048/93‑1 (MQ=255)
                                                       aCGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTAcc                                     <  1:6497962/93‑1 (MQ=255)
                                                       aCGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTAcc                                     <  1:625174/93‑1 (MQ=255)
                                                       aCGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTAcc                                     >  2:6523789/1‑93 (MQ=255)
                                                        cGGTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCt                                    <  1:1642075/93‑1 (MQ=255)
                                                          gTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTc                                  <  1:2579464/93‑1 (MQ=255)
                                                          gTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTc                                  <  1:1904118/93‑1 (MQ=255)
                                                          gTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTc                                  >  2:2483623/1‑93 (MQ=255)
                                                          gTTACCCGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTc                                  <  1:5538620/93‑1 (MQ=255)
                                                            tACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGt                                >  1:2739805/1‑93 (MQ=255)
                                                              cccGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTtg                              <  1:6688160/93‑1 (MQ=255)
                                                              cccGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTtg                              <  1:6071576/93‑1 (MQ=255)
                                                               ccGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTtgt                             >  2:492646/1‑93 (MQ=255)
                                                               ccGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTtgt                             <  2:5869001/93‑1 (MQ=255)
                                                                cGGCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTtgtg                            >  2:4746713/1‑93 (MQ=255)
                                                                 ggCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGCGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGa                           <  1:7373771/93‑1 (MQ=255)
                                                                  gCAGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAg                          <  2:4090844/93‑1 (MQ=255)
                                                                    aGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCt                        <  1:3965538/93‑1 (MQ=255)
                                                                    aGGGAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCt                        >  1:4243014/1‑93 (MQ=255)
                                                                      ggAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTtc                      >  1:4578183/1‑93 (MQ=255)
                                                                      ggAGATCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTtc                      <  1:5606981/93‑1 (MQ=255)
                                                                       gagaTCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTtca                     <  2:5655288/93‑1 (MQ=255)
                                                                       gagaTCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTtca                     >  2:6994916/1‑93 (MQ=255)
                                                                          aTCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTtcactc                  <  2:5081703/93‑1 (MQ=255)
                                                                          aTCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTtcactc                  >  2:1089041/1‑93 (MQ=255)
                                                                          aTCATTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTtcactc                  >  2:2204251/1‑93 (MQ=255)
                                                                             aTTTCACCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTcaca               <  2:1783769/93‑1 (MQ=255)
                                                                                  aCCCTGCATGCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACAtt          >  1:946578/1‑93 (MQ=255)
                                                                                   cccTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTg         >  1:142607/1‑93 (MQ=255)
                                                                                   cccTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTg         >  1:2344914/1‑93 (MQ=255)
                                                                                      tgcatgcaATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAAc      <  2:2662527/93‑1 (MQ=255)
                                                                                      tgcatgcaATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAAc      >  2:5485840/1‑93 (MQ=255)
                                                                                      tgcatgcaATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAAc      >  2:5009771/1‑93 (MQ=255)
                                                                                       gCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACt     <  1:7497985/93‑1 (MQ=255)
                                                                                         aaaCAATCGGCCAGGCCGTGAATGGGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGa   <  1:4703794/91‑1 (MQ=255)
                                                                                          tACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAg  >  2:6638457/1‑93 (MQ=255)
                                                                                           |                                                                                           
GACCCATAACAACCGCACCTACAGAGTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAG  >  NC_000913/3917254‑3917436

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