Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,558,583 A→G V124A (GTC→GCC)  b1486 ← putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,558,5830AG99.0% 300.8 / ‑4.7 97V124A (GTC→GCC) b1486putative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (50/46);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

ATGACGCCCAGCCCAAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATA  >  NC_000913/1558493‑1558672
                                                                                          |                                                                                         
aTGCCGCCCAGCCCAAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCt                                                                                         <  1:3266110/93‑1 (MQ=255)
 tGACCCCCAGCCCAAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTg                                                                                        <  2:2488996/93‑1 (MQ=255)
   aCGCCCAGCCCAAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCa                                                                                      <  2:980430/93‑1 (MQ=255)
   aCGCCCAGCCCAAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCa                                                                                      <  2:1102523/93‑1 (MQ=255)
    cGCCCAGCCCAAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCag                                                                                     <  2:3272365/93‑1 (MQ=255)
       ccagcccaAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCagaga                                                                                  <  2:4663029/93‑1 (MQ=255)
       ccagcccaAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCagaga                                                                                  >  2:3819917/1‑93 (MQ=255)
       ccagcccaAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCagaga                                                                                  >  2:1073352/1‑93 (MQ=255)
        cgacccaAGCCAGAACGCAGGTTTGGAGATTCGGGTAATGGCCATTATTCGCACCGGATGATCAATCCAGCGATTTCGCCAGGCTGCagagag                                                                                 <  2:1001628/90‑1 (MQ=255)
        cagcccaAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCagagag                                                                                 <  1:4467020/93‑1 (MQ=255)
         agcccaAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGt                                                                                <  2:4246699/93‑1 (MQ=255)
         agcccaAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGt                                                                                <  2:1961490/93‑1 (MQ=255)
          gcccaAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTa                                                                               <  1:623492/93‑1 (MQ=255)
          gcccaAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTa                                                                               >  2:606019/1‑93 (MQ=255)
          gcccaAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTa                                                                               <  1:2158739/93‑1 (MQ=255)
            ccaAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCTCCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATg                                                                             >  1:223584/1‑93 (MQ=255)
              aaGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGcc                                                                           >  1:449485/1‑93 (MQ=255)
              aaGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGcc                                                                           <  2:1521895/93‑1 (MQ=255)
               aGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGAGCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGccc                                                                          <  2:2574847/93‑1 (MQ=255)
                 ccAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAg                                                                        >  2:3883105/1‑93 (MQ=255)
                  cAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGc                                                                       <  1:1293070/93‑1 (MQ=255)
                   aGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCg                                                                      >  1:3086346/1‑93 (MQ=255)
                    gAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCgg                                                                     <  2:4203582/93‑1 (MQ=255)
                      aCGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGa                                                                   <  2:3862783/93‑1 (MQ=255)
                      aCGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGAGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGa                                                                   <  1:29193/93‑1 (MQ=255)
                       cGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGAt                                                                  >  2:2495151/1‑93 (MQ=255)
                       cGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGAt                                                                  <  1:3232247/93‑1 (MQ=255)
                       cGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGAt                                                                  <  1:4333194/93‑1 (MQ=255)
                       cGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGAt                                                                  <  1:3805227/93‑1 (MQ=255)
                        gCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATg                                                                 >  1:3734410/1‑93 (MQ=255)
                         cAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGc                                                                <  2:3346529/93‑1 (MQ=255)
                          aGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGcc                                                               <  1:4342613/93‑1 (MQ=255)
                           ggTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCa                                                              >  2:787953/1‑93 (MQ=255)
                             tgtgGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAAt                                                            >  2:1512220/1‑93 (MQ=255)
                               tgGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCa                                                          >  2:1232053/1‑93 (MQ=255)
                               tgGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCa                                                          >  1:666960/1‑93 (MQ=255)
                               tgGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCa                                                          <  1:2913889/93‑1 (MQ=255)
                                 gagaTTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAgg                                                        <  1:4404637/93‑1 (MQ=255)
                                  agaTTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGa                                                       <  2:1339156/93‑1 (MQ=255)
                                  agaTTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGa                                                       >  2:2129211/1‑93 (MQ=255)
                                   gaTTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCTCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAg                                                      <  1:2520601/93‑1 (MQ=255)
                                   gaTTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAg                                                      >  1:2991163/1‑93 (MQ=255)
                                   gaTTCCGGTAATGGCCATTATTAGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAg                                                      >  1:393656/1‑93 (MQ=255)
                                    aTTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGt                                                     <  2:2511424/93‑1 (MQ=255)
                                     ttCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTg                                                    >  1:3111606/1‑93 (MQ=255)
                                       ccGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGcc                                                  >  2:2863043/1‑93 (MQ=255)
                                          gTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCcagc                                               <  1:4498039/93‑1 (MQ=255)
                                           tAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCcagca                                              >  1:3120285/1‑93 (MQ=255)
                                            aaTGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCTCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCcagcag                                             <  1:1538394/93‑1 (MQ=255)
                                             aTGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCcagcagc                                            >  1:3784079/1‑93 (MQ=255)
                                              tGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCcagcagca                                           >  1:3297737/1‑93 (MQ=255)
                                               ggCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCcagcagcag                                          >  1:4422771/1‑93 (MQ=255)
                                                gCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAgg                                         <  2:2234082/93‑1 (MQ=255)
                                                gCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAgg                                         >  2:4240435/1‑93 (MQ=255)
                                                  cgttattCGCGCTGGATGATCTAGCCAGCGACTGCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATACCAATCAGGACTGCCAGGAGCAGGGc                                       <  2:3734646/91‑1 (MQ=255)
                                                  cATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGc                                       >  1:3927489/1‑93 (MQ=255)
                                                  cATTATTCACACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGc                                       >  2:4572767/1‑93 (MQ=255)
                                                   attattCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCa                                      <  1:2214959/93‑1 (MQ=255)
                                                   attattCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCa                                      >  1:2911560/1‑93 (MQ=255)
                                                   attattCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCa                                      <  1:1930789/93‑1 (MQ=255)
                                                    ttattCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCAc                                     >  1:220257/1‑93 (MQ=255)
                                                    ttattCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCAc                                     >  1:1955338/1‑93 (MQ=255)
                                                        tCGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCaaa                                 <  1:1404619/93‑1 (MQ=255)
                                                         cGCACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCaaaa                                <  2:3866290/93‑1 (MQ=255)
                                                           cACCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGc                              >  1:2555200/1‑93 (MQ=255)
                                                            aCCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGcc                             <  2:2334520/93‑1 (MQ=255)
                                                            aCCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGcc                             >  1:1454937/1‑93 (MQ=255)
                                                            aCCAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGcc                             >  1:2022148/1‑93 (MQ=255)
                                                             ccAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCa                            >  1:1495643/1‑93 (MQ=255)
                                                              cAGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAg                           >  2:2713826/1‑93 (MQ=255)
                                                               aGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGc                          <  1:1853418/93‑1 (MQ=255)
                                                               aGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGc                          >  2:1265151/1‑93 (MQ=255)
                                                               aGATGATCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGc                          >  1:7572/1‑93 (MQ=255)
                                                                gatgatCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCt                         >  1:1413294/1‑93 (MQ=255)
                                                                 atgatCCAGCCAGCTATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTc                        <  2:2006100/93‑1 (MQ=255)
                                                                 atgatCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTc                        >  1:389901/1‑93 (MQ=255)
                                                                    atCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTcccg                     >  1:2833469/1‑91 (MQ=255)
                                                                    atCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAg                     <  1:4106735/93‑1 (MQ=255)
                                                                     tCCAGCCAGCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGc                    >  1:1444655/1‑93 (MQ=255)
                                                                      ccagccagCGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCg                   <  2:119282/93‑1 (MQ=255)
                                                                         gccagcGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCg                >  2:3094517/1‑93 (MQ=255)
                                                                         gccagcGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCg                >  2:3416997/1‑93 (MQ=255)
                                                                          gcagcGATTTAGCCAGGCTGCAGAGAGTATGCCGAGCGGGATGCCAATCAGGAGCGGCAGCAGCAGGGTACCAAAAGCCAGATTCAGGGTCGc               >  2:1541640/2‑93 (MQ=255)
                                                                          ccagcGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGc               <  1:2328512/93‑1 (MQ=255)
                                                                           cagcGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGcc              <  2:3238620/93‑1 (MQ=255)
                                                                            agcGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCg             >  2:333282/1‑93 (MQ=255)
                                                                            agcGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCg             >  2:330410/1‑93 (MQ=255)
                                                                             gcGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCgg            <  1:107733/93‑1 (MQ=255)
                                                                             gcGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCgg            <  1:3449752/93‑1 (MQ=255)
                                                                             gcGATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCgg            >  1:2763359/1‑93 (MQ=255)
                                                                               gATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGa          >  2:850680/1‑93 (MQ=255)
                                                                               gATTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGa          >  2:3104185/1‑93 (MQ=255)
                                                                                aTTTCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGaa         >  2:283960/1‑93 (MQ=255)
                                                                                  ttCGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGaaaa       >  2:3717783/1‑93 (MQ=255)
                                                                                   tCGCCAGGCTGCAGAGAGTAGGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGaaaaa      <  1:1975932/93‑1 (MQ=255)
                                                                                    cGCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAat     >  1:4048960/1‑93 (MQ=255)
                                                                                     gCCAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAata    <  1:2575352/93‑1 (MQ=255)
                                                                                       cAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAAtata  >  2:2577824/1‑93 (MQ=255)
                                                                                       cAGGCTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAAtata  <  2:3931009/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATGACGCCCAGCCCAAGCCAGAACGCAGGTGTGGAGATTCCGGTAATGGCCATTATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATA  >  NC_000913/1558493‑1558672

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