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,1730TA96.6% 278.4 / 0.4 88A40A (GCA→GCTyfdYDUF2545 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (48/37);  total (49/39)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

AGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGC  >  NC_000913/2495076‑2495268
                                                                                                 |                                                                                               
aGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGc                                                                                               <  2:874940/100‑1 (MQ=255)
  cGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAg                                                                                             >  1:498004/1‑100 (MQ=255)
       catcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCg                                                                                           <  2:2478272/97‑1 (MQ=255)
       catcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAgg                                                                                        <  1:1819059/100‑1 (MQ=255)
       catcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAgg                                                                                        <  2:397084/100‑1 (MQ=255)
       catcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAgg                                                                                        >  2:710611/1‑100 (MQ=255)
        atcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGa                                                                                       >  1:988513/1‑100 (MQ=255)
        atcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGa                                                                                       >  1:1451861/1‑100 (MQ=255)
         tcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCGCGTCACTAGTTCATTACCGCCAGTTAACAACATGGTGAAGTAATAAATGAGAGCAGCGAGGa                                                                                       >  2:1704166/1‑99 (MQ=255)
          caAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATc                                                                                     >  1:1270468/1‑100 (MQ=255)
          caAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATc                                                                                     >  2:1246066/1‑100 (MQ=255)
            aGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAc                                                                                   >  1:2303098/1‑100 (MQ=255)
            aGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAc                                                                                   >  2:1811013/1‑100 (MQ=255)
              ccACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCa                                                                                 <  1:1508910/100‑1 (MQ=255)
              ccACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCa                                                                                 <  1:364360/100‑1 (MQ=255)
              ccACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCa                                                                                 >  1:1096987/1‑100 (MQ=255)
               cacaAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGCAGTAATAAATGAGAGCAGCGAGGATCACCAt                                                                                >  1:1941473/1‑100 (MQ=255)
               cacaAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCAt                                                                                >  1:954489/1‑100 (MQ=255)
               cacaAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGCAATAAATGAGAGCAGCGAGGATCACCAt                                                                                >  2:2013998/1‑100 (MQ=255)
                   aGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCg                                                                            <  1:1336498/100‑1 (MQ=255)
                   aGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCg                                                                            >  1:1064884/1‑100 (MQ=255)
                     ccGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCgaa                                                                          >  2:668464/1‑100 (MQ=255)
                      cGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGCGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCgaag                                                                         >  1:919328/1‑100 (MQ=255)
                       tgggCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGCTAACCACATGGTGAAGTAATAAACGAGAGCAGCGAGGATCACCATGCCgaag                                                                         >  2:2192730/2‑99 (MQ=255)
                         ggCAAGAAATATAAATATCCCCGTCACTAGTTTATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGaagaaa                                                                      <  1:1792488/100‑1 (MQ=255)
                           cAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGaagaaaga                                                                    >  1:2144726/1‑100 (MQ=255)
                           cAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGaagaaaga                                                                    >  2:823536/1‑100 (MQ=255)
                           cAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGaagaaaga                                                                    >  1:127909/1‑100 (MQ=255)
                             aGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAaagaag                                                                  <  1:2010939/100‑1 (MQ=255)
                             aGAAACATAAATATCCCCGTCACTAGGTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAagaaagaa                                                                   >  2:887047/1‑99 (MQ=255)
                                aaCATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAgaga                                                               >  1:697352/1‑100 (MQ=255)
                                aaCATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAgaga                                                               >  1:618278/1‑100 (MQ=255)
                                  cATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAcc                                                             >  2:1933700/1‑100 (MQ=255)
                                  cATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAcc                                                             <  1:1927700/100‑1 (MQ=255)
                                    tAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAccgc                                                           >  1:710208/1‑100 (MQ=255)
                                     aaaTATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAccgcc                                                          <  2:33771/100‑1 (MQ=255)
                                      aaTATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAccgccg                                                         <  1:1466931/100‑1 (MQ=255)
                                       atatCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGTTCACCATGCCGAAGAAAGAAGAGACCGCCGa                                                        <  2:2152796/100‑1 (MQ=255)
                                         atCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATa                                                      >  1:2560706/1‑100 (MQ=255)
                                           ccccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATAcc                                                    >  2:1449328/1‑100 (MQ=255)
                                            cccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCa                                                   <  2:437989/100‑1 (MQ=255)
                                             ccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCAc                                                  >  1:833727/1‑100 (MQ=255)
                                             ccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCAc                                                  <  2:2293675/100‑1 (MQ=255)
                                              cGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACg                                                 <  2:992942/100‑1 (MQ=255)
                                              cGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACg                                                 >  1:1779703/1‑100 (MQ=255)
                                                tCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGtt                                               >  2:1562246/1‑100 (MQ=255)
                                                 cACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTc                                              <  2:2041743/100‑1 (MQ=255)
                                                  aCTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCa                                             <  1:2245068/100‑1 (MQ=255)
                                                  aCTAGTTCATTACCGCCAGTTAACCACACGGTGAAGTAATAAATGAGAGCAGCGAGGATCCTCATGCCGAAAAAAGAAGAGAGCGCTGATACCa                                                   <  2:1343231/94‑1 (MQ=255)
                                                     aGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTa                                          >  1:2348700/1‑100 (MQ=255)
                                                     aGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTa                                          <  2:168040/100‑1 (MQ=255)
                                                      gTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAc                                         >  1:2231360/1‑100 (MQ=255)
                                                      gTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAc                                         >  1:1194394/1‑100 (MQ=255)
                                                      gTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAc                                         <  2:1403744/100‑1 (MQ=255)
                                                       ttCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAcn                                        <  2:618271/100‑2 (MQ=255)
                                                       ttCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAcc                                        >  2:817269/1‑100 (MQ=255)
                                                         cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTg                                      <  2:1194386/100‑1 (MQ=255)
                                                             aCCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCg                                  <  1:1325960/100‑1 (MQ=255)
                                                              ccGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGa                                 >  1:2002457/1‑100 (MQ=255)
                                                              ccGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGa                                 >  1:462398/1‑100 (MQ=255)
                                                              ccGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGa                                 <  1:1045074/100‑1 (MQ=255)
                                                              ccGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGACCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGa                                 >  2:1513092/1‑100 (MQ=255)
                                                                gCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATg                               <  2:127907/100‑1 (MQ=255)
                                                                    gTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGc                           <  1:1792883/100‑1 (MQ=255)
                                                                      tAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATAACACGTTCAGTACCTGACCGATGTAGCCg                         >  1:2137893/1‑100 (MQ=255)
                                                                       aaCCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGccac                        >  2:1762260/1‑98 (MQ=255)
                                                                         ccACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCtt                      >  1:1570171/1‑100 (MQ=255)
                                                                          cacaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTa                     <  1:2255678/100‑1 (MQ=255)
                                                                          cacaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTa                     >  1:1872737/1‑100 (MQ=255)
                                                                          cacaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTa                     >  1:191133/1‑100 (MQ=255)
                                                                           acaTGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTAc                    <  2:833720/100‑1 (MQ=255)
                                                                           acaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGCTGTAGCCGCtt                      >  2:189144/1‑98 (MQ=255)
                                                                           acaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTAc                    >  1:855629/1‑100 (MQ=255)
                                                                           acaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTAc                    <  2:1735442/100‑1 (MQ=255)
                                                                              tGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCa                 <  2:2033511/100‑1 (MQ=255)
                                                                              tGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCa                 >  1:2457888/1‑100 (MQ=255)
                                                                              tGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCa                 <  1:1292583/100‑1 (MQ=255)
                                                                               ggTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAg                <  1:1074928/100‑1 (MQ=255)
                                                                                gTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACAACGTTCCGTACCTGACCGATGTAGCCGCTTACGCAGa               >  2:1993994/1‑100 (MQ=255)
                                                                                 tGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGAc              >  2:1651074/1‑100 (MQ=255)
                                                                                  gAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGAc              <  1:459829/99‑1 (MQ=255)
                                                                                  gAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACa             >  1:950325/1‑100 (MQ=255)
                                                                                  gAAGTAATAAAGGAGAGCAGCGAGGATAACCATGCCGAAGAAAGAAGATACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACa             <  2:248369/100‑1 (MQ=255)
                                                                                   aaGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACaa            <  2:2002236/100‑1 (MQ=255)
                                                                                     gTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATa          <  1:2379651/100‑1 (MQ=255)
                                                                                       aataaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATaca        >  2:950728/1‑100 (MQ=255)
                                                                                          aaaTGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAg     <  2:1225513/100‑1 (MQ=255)
                                                                                             tGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGc  <  2:1096979/100‑1 (MQ=255)
                                                                                                 |                                                                                               
AGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGC  >  NC_000913/2495076‑2495268

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