Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,132 (C)8→9 pseudogene (3955/5422 nt) yneO ← AIDA‑I family autotransporter YneO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,1241.C99.5% 782.1 / ‑3.4 194pseudogene (3963/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (82/111);  total (82/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 226 total aligned reads displayed.

AACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATT  >  NC_000913/1594024‑1594219
                                                                                                     |                                                                                               
 aCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑c                                                                                                >  1:2083618/1‑101 (MQ=255)
 aCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑c                                                                                                <  1:2083263/101‑1 (MQ=255)
    gTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccc                                                                                             <  2:2083324/101‑1 (MQ=255)
       aGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccc                                                                                          >  2:2083620/1‑101 (MQ=255)
          aCATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccca                                                                                       >  1:2083622/1‑99 (MQ=255)
          aCATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccca                                                                                       >  2:2083624/1‑99 (MQ=255)
           cATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccccag                                                                                      <  2:2083332/101‑4 (MQ=255)
           cATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccccag                                                                                      <  2:2083286/101‑4 (MQ=255)
           cATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccccag                                                                                      <  1:2083381/101‑4 (MQ=255)
            aTTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccccagg                                                                                     >  1:2083628/1‑97 (MQ=255)
                 ccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATaa                                                                                 >  2:2083630/1‑101 (MQ=255)
                   ggTTTCCCCGCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACg                                                                               >  1:2083632/1‑101 (MQ=255)
                   ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACg                                                                               <  1:2083168/101‑1 (MQ=255)
                     tttCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACg                                                                               >  2:2083633/1‑99 (MQ=255)
                         cccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccc                                                                                        >  1:2083637/1‑84 (MQ=255)
                           cccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATc                                                                       <  2:2083470/101‑1 (MQ=255)
                            ccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATc                                                                       >  2:2083640/1‑100 (MQ=255)
                                cGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTg                                                                  >  2:2083642/1‑101 (MQ=255)
                                  cTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCa                                                                >  1:2083645/1‑101 (MQ=255)
                                   tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                               >  1:2083648/1‑101 (MQ=255)
                                   tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                               <  1:2083424/101‑1 (MQ=255)
                                    aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgt                                                              >  1:2083650/1‑101 (MQ=255)
                                     ttCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtg                                                             <  2:2083190/101‑1 (MQ=255)
                                      tCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtgg                                                            >  1:2083653/1‑101 (MQ=255)
                                       cATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtggt                                                           >  2:2083654/1‑101 (MQ=255)
                                        aTTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAAc                                                                                >  1:2083656/1‑79 (MQ=255)
                                           tGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACa                                                       >  1:2083658/1‑101 (MQ=255)
                                            gAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAt                                                      >  2:2083660/1‑101 (MQ=255)
                                            gAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccc                                                                                           <  2:2083240/63‑1 (MQ=255)
                                             aaCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgt                                                              <  2:2083257/92‑1 (MQ=255)
                                             aaCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAtt                                                     <  2:2083369/101‑1 (MQ=255)
                                              aCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAttt                                                    <  1:2083418/101‑1 (MQ=255)
                                               ccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTc                                                   <  2:2083452/101‑1 (MQ=255)
                                                cGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                                  <  1:2083366/101‑1 (MQ=255)
                                                cGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccc                                                                                        >  2:2083663/1‑61 (MQ=255)
                                                 gTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTccc                                                 >  1:2083664/1‑101 (MQ=255)
                                                 gTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccccc                                                                                       >  1:2083668/1‑60 (MQ=255)
                                                   gccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGa                                               >  2:2083671/1‑101 (MQ=255)
                                                   gccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGa                                               >  1:2083670/1‑101 (MQ=255)
                                                   gccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGa                                               <  1:2083245/101‑1 (MQ=255)
                                                     cgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACa                                             <  2:2083389/101‑1 (MQ=255)
                                                     cgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACa                                             <  1:2083204/101‑1 (MQ=255)
                                                       ccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                           >  2:2083675/1‑101 (MQ=255)
                                                       ccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                           >  1:2083676/1‑101 (MQ=255)
                                                       ccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                           <  1:2083284/101‑1 (MQ=255)
                                                         aTTCTACAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcac                                         >  1:2083679/1‑101 (MQ=255)
                                                          ttCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATca                                          <  2:2083360/99‑1 (MQ=255)
                                                           tctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                       <  2:2083335/101‑1 (MQ=255)
                                                            ctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtgg                                                            <  1:2083091/79‑1 (MQ=255)
                                                            ctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccc                                                                                              <  1:2083396/44‑1 (MQ=255)
                                                             tcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                     <  2:2083169/101‑1 (MQ=255)
                                                              ccAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGc                                    >  1:2083682/1‑101 (MQ=255)
                                                              ccAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGc                                    <  1:2083373/101‑1 (MQ=255)
                                                               cAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCt                                   <  1:2083358/101‑1 (MQ=255)
                                                                aGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTg                                  <  2:2083310/101‑1 (MQ=255)
                                                                 ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGt                                 <  1:2083322/101‑1 (MQ=255)
                                                                 ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGt                                 <  2:2083399/101‑1 (MQ=255)
                                                                  gTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                >  1:2083684/1‑101 (MQ=255)
                                                                   ttGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTa                               >  2:2083687/1‑101 (MQ=255)
                                                                   ttGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTa                               >  1:2083686/1‑101 (MQ=255)
                                                                     gccgccGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                             <  2:2083423/101‑1 (MQ=255)
                                                                      ccgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATgg                            <  1:2083347/101‑1 (MQ=255)
                                                                       cgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGa                           >  1:2083690/1‑101 (MQ=255)
                                                                        gccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGat                          >  1:2083693/1‑101 (MQ=255)
                                                                        gccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGat                          >  2:2083694/1‑101 (MQ=255)
                                                                         ccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGata                         >  2:2083697/1‑101 (MQ=255)
                                                                          cgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATTTCACCTGCAgtggtg                                                          >  1:2083700/1‑67 (MQ=255)
                                                                          cgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGatat                        >  1:2083698/1‑101 (MQ=255)
                                                                           gATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATAtt                       <  2:2083412/101‑1 (MQ=255)
                                                                           gATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATAtt                       <  1:2083454/101‑1 (MQ=255)
                                                                              aGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAt                    >  2:2083705/1‑101 (MQ=255)
                                                                               gTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtg                   <  1:2083499/101‑1 (MQ=255)
                                                                                  gTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATgtgt                >  1:2083706/1‑101 (MQ=255)
                                                                                  gTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATgtgt                <  1:2083431/101‑1 (MQ=255)
                                                                                   ttttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGtt               <  2:2083348/101‑1 (MQ=255)
                                                                                    tttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                                  <  2:2083264/65‑1 (MQ=255)
                                                                                    tttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg              <  2:2083302/101‑1 (MQ=255)
                                                                                       aGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACt           <  2:2083532/101‑1 (MQ=255)
                                                                                       aGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACt           <  1:2083483/101‑1 (MQ=255)
                                                                                       aGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACt           <  2:2083380/101‑1 (MQ=255)
                                                                                        gCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACtt          <  1:2083311/101‑1 (MQ=255)
                                                                                        gCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACtt          <  2:2083143/101‑1 (MQ=255)
                                                                                        gCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACtt          >  1:2083709/1‑101 (MQ=255)
                                                                                         cGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCGCCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTg         <  2:2083573/101‑1 (MQ=255)
                                                                                         cGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTg         >  2:2083710/1‑101 (MQ=255)
                                                                                         cGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTg         <  1:2083464/101‑1 (MQ=255)
                                                                                          gACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGc        >  1:2083712/1‑101 (MQ=255)
                                                                                          gACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGc        <  1:2083392/101‑1 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCAcc                                                                    <  2:2083243/39‑1 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtgt                  <  1:2083503/89‑1 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa      >  1:2083716/1‑101 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa      <  2:2083273/101‑1 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa      >  2:2083715/1‑101 (MQ=255)
                                                                                              cGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTg                                                                  <  2:2083434/39‑1 (MQ=255)
                                                                                              cGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtg                                                             <  1:2083511/44‑1 (MQ=255)
                                                                                              cGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTc                                                   <  1:2083294/54‑1 (MQ=255)
                                                                                              cGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc       <  1:2083356/98‑1 (MQ=255)
                                                                                               gCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTccc                                                 <  2:2083341/55‑1 (MQ=255)
                                                                                               gCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCatat   >  1:2083721/1‑101 (MQ=255)
                                                                                                cAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATAtt  <  1:2083403/101‑1 (MQ=255)
                                                                                                     |                                                                                               
AACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATT  >  NC_000913/1594024‑1594219

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