Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,495,173 T→A 100% A40A (GCA→GCT yfdY ← DUF2545 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,495,1730TA97.2% 265.4 / ‑1.6 85A40A (GCA→GCTyfdYDUF2545 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (48/34);  total (49/36)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGCG  >  NC_000913/2495076‑2495269
                                                                                                 |                                                                                                
aGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTAACGCCAGTTAACCACATGATGAAGTAATAAATGAGAGc                                                                                                >  1:950998/1‑100 (MQ=255)
 gCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCa                                                                                               >  2:702970/1‑100 (MQ=255)
 gCGAATCATCAAGACACAAGCCGGGGCAAGAAACATAATTATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCa                                                                                               >  2:1869938/1‑100 (MQ=255)
   gAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGagcagc                                                                                             >  2:2427554/1‑100 (MQ=255)
    aaTCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCg                                                                                            <  2:2374875/100‑1 (MQ=255)
      tcatcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAg                                                                                          <  2:74829/100‑1 (MQ=255)
       catcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAgg                                                                                         >  2:2965944/1‑100 (MQ=255)
         tcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGa                                                                                        >  2:558365/1‑99 (MQ=255)
         tcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGAt                                                                                       <  1:900158/100‑1 (MQ=255)
          caAGTCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATc                                                                                      >  2:2913469/1‑100 (MQ=255)
          caAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATc                                                                                      >  1:482344/1‑100 (MQ=255)
           aaGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCa                                                                                     <  2:2779127/100‑1 (MQ=255)
            aGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAc                                                                                    >  1:1955396/1‑100 (MQ=255)
            aGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAc                                                                                    <  1:2917346/100‑1 (MQ=255)
            aGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAc                                                                                    >  2:1796091/1‑100 (MQ=255)
            aGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACAGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAc                                                                                    >  2:272995/1‑100 (MQ=255)
             gCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAcc                                                                                   >  2:879543/1‑100 (MQ=255)
               cacaAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCAt                                                                                 >  1:1466403/1‑100 (MQ=255)
                  aaGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGcc                                                                              >  2:2462143/1‑100 (MQ=255)
                  aaGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGcc                                                                              <  2:99293/100‑1 (MQ=255)
                  aaGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGcc                                                                              <  1:2678249/100‑1 (MQ=255)
                  aaGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGcc                                                                              <  2:1196912/100‑1 (MQ=255)
                  aaGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAAGATcggaagagc                                                                              >  1:2366904/1‑92 (MQ=255)
                   aGCCGGGGCAAGAAACATAAATATTCCCGTCACTAGTTCATTAACGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCg                                                                             <  1:410150/100‑1 (MQ=255)
                   aGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCg                                                                             >  1:1813368/1‑100 (MQ=255)
                    gCCGGGGCAAGACACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGCGAAGTAATAAATGAGAGCAGCAAGGATCACCATGCCga                                                                            >  2:1770755/1‑100 (MQ=255)
                    gCCGGGGCAAGAAACATATATATCCCCGTCACTAGTTTATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCga                                                                            <  1:87181/100‑1 (MQ=255)
                    gCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCga                                                                            >  2:1746070/1‑100 (MQ=255)
                          gCAAGAAAGATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCAGATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGaagaaag                                                                      <  2:862541/100‑1 (MQ=255)
                          gCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGaagaaag                                                                      >  1:714740/1‑100 (MQ=255)
                           cAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGaagaaaga                                                                     >  1:962891/1‑100 (MQ=255)
                           cAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGaagaaaga                                                                     >  1:1623363/1‑100 (MQ=255)
                              gAAACATAAATATCACCATCACTAATTCATTAACACCAGTTAAACACATGGAAAAGTAATAAATGAGAGCAGCGAAGATCACCATACCGAAGAAagaaga                                                                  >  2:5815/1‑100 (MQ=255)
                               aaaCATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAagag                                                                 <  2:2327894/100‑1 (MQ=255)
                               aaaCATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAagag                                                                 <  2:296168/100‑1 (MQ=255)
                                  cATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACAACATGGAGAAGTAATAAATGAGAGAAGCAAGGATCACCATGCCGAATAAAGAAGAGAcc                                                              >  2:1100696/1‑100 (MQ=255)
                                   aTAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAccg                                                             >  2:259815/1‑100 (MQ=255)
                                   aTAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAccg                                                             >  1:856740/1‑100 (MQ=255)
                                    tAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAccgc                                                            >  1:1455590/1‑100 (MQ=255)
                                      cgatctcccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAccgccg                                                          <  2:2413903/94‑1 (MQ=255)
                                       atatCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGa                                                         >  1:266277/1‑100 (MQ=255)
                                        tatCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGAt                                                        >  2:1712177/1‑100 (MQ=255)
                                         atCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCATCGAGGATCACCATGCCGAAGAAAAAAGAGACCGCAGATa                                                       >  1:461837/1‑100 (MQ=255)
                                           ccccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATAcc                                                     >  2:713075/1‑100 (MQ=255)
                                            cccGTCACTAGTTCATTACCGTCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCa                                                    <  1:520963/100‑1 (MQ=255)
                                            cccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTTAAGTAATAAATGAGAGCAGCGAGGATCACCATGCTGAAGTAAGAAGAGACCGCCGATa                                                       <  2:644045/97‑1 (MQ=255)
                                            cccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGagatcg                                                    >  1:2413932/1‑95 (MQ=255)
                                                 cACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTc                                               >  2:2625028/1‑100 (MQ=255)
                                                 cACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTc                                               <  1:1770965/100‑1 (MQ=255)
                                                   cTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAg                                             >  1:10708/1‑100 (MQ=255)
                                                     aGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTa                                           >  2:2556798/1‑100 (MQ=255)
                                                     aGTTCATTACCGCCAGTTAACCACATGATTAAGTAATAAATAAGAGCAGCGAGTATCACCATGCCGAATAAAGAAGATACCGCCGATACCACTTTCAGt                                            <  2:1412184/99‑1 (MQ=255)
                                                      gTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAATAAAGAAGATACCGCCTATACCACGTTCAGTAc                                          <  2:2137310/100‑1 (MQ=255)
                                                       ttcacgtccGCCAGTGAACCACATGGTGAAGTAAGAAGGGAGTGCCGCGTTGATCACCATGCCGACGAAAGAAGAGTCCGCCGATACCGCGttga                                              <  2:1739854/88‑3 (MQ=255)
                                                       tgcATTACCGCTAGTTAACCACATGGTTAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAcc                                         <  1:1234646/98‑1 (MQ=255)
                                                         cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTg                                       <  1:2914939/100‑1 (MQ=255)
                                                         cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTg                                       <  2:1980507/100‑1 (MQ=255)
                                                           ttACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATTCCGAAGAAAGAAGAGACCTCCGATACCACGTTCAGTACCTGAc                                     <  2:2454483/100‑1 (MQ=255)
                                                           ttACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGAc                                     >  2:2980981/1‑100 (MQ=255)
                                                             aCCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCg                                   <  1:2663016/100‑1 (MQ=255)
                                                               cGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAAACCGCCGATACCACGTTCAGTACCTGACCGAt                                 >  1:195368/1‑100 (MQ=255)
                                                                gctaGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATg                                <  2:2034364/97‑1 (MQ=255)
                                                                gCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATg                                <  2:206370/100‑1 (MQ=255)
                                                                 ccAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGt                               >  1:2603470/1‑100 (MQ=255)
                                                                  cAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTa                              <  2:1955373/100‑1 (MQ=255)
                                                                  cAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTa                              >  1:2386873/1‑100 (MQ=255)
                                                                      tAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCg                          >  1:1075105/1‑100 (MQ=255)
                                                                       aaCCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCg                          >  2:471709/1‑99 (MQ=255)
                                                                        aCCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCt                        <  2:2738152/100‑1 (MQ=255)
                                                                         ccACATGGTGAAGTAATAAATGAGAGCAGCGATGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCt                        <  1:490297/99‑1 (MQ=255)
                                                                         ccACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCtt                       >  1:1200756/1‑100 (MQ=255)
                                                                          cacaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTa                      >  1:3009459/1‑100 (MQ=255)
                                                                          cacaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCAGCTTa                      >  2:1736691/1‑100 (MQ=255)
                                                                               ggTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAg                 <  2:1517294/100‑1 (MQ=255)
                                                                                gTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGa                <  1:2965977/100‑1 (MQ=255)
                                                                                 tGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGCCCGATGTAGCCGATTACGCa                  <  2:2578934/97‑1 (MQ=255)
                                                                                 tGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGAc               <  1:308549/100‑1 (MQ=255)
                                                                                 tGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGAc               <  1:305463/100‑1 (MQ=255)
                                                                                    aGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAAt            >  1:1354952/1‑100 (MQ=255)
                                                                                      taataaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATac          <  1:800288/100‑1 (MQ=255)
                                                                                      taataaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATac          >  2:879498/1‑100 (MQ=255)
                                                                                       aataaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTAACGCAGACAATaca         >  1:212738/1‑100 (MQ=255)
                                                                                          aaaTGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAg      >  2:338763/1‑100 (MQ=255)
                                                                                           aaTGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAgg     >  1:576175/1‑100 (MQ=255)
                                                                                             tGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGc   >  2:2576328/1‑100 (MQ=255)
                                                                                             tGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGc   >  1:1666795/1‑100 (MQ=255)
                                                                                              gagaGCAGCGAGGATCACCATGTCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGCg  >  1:2516204/1‑100 (MQ=255)
                                                                                              gagaGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTTACCGATGTAGCCGCTTACGCAGACAATACACAGGGCg  <  2:1619397/100‑1 (MQ=255)
                                                                                                 |                                                                                                
AGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGCG  >  NC_000913/2495076‑2495269

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