Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,1241.C58.8% 595.1 / inf 250pseudogene (3963/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base . (51/52);  new base C (65/82);  total (116/134)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 291 total aligned reads displayed.

AACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGG  >  NC_000913/1594024‑1594224
                                                                                                     |                                                                                                    
aaCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTg                                                                                                       >  2:3065599/1‑101 (MQ=255)
   ggTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccc                                                                                                   <  1:3065261/101‑1 (MQ=255)
    gTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccc                                                                                                  <  1:3065087/101‑1 (MQ=255)
    gTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccc                                                                                                  <  2:3064833/101‑1 (MQ=255)
      cAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccc                                                                                                >  2:3065604/1‑101 (MQ=255)
      cAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccc                                                                                                <  1:3065229/101‑1 (MQ=255)
       aGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccc                                                                                               <  2:3065024/101‑1 (MQ=255)
        gTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccc                                                                                              <  1:3065001/101‑1 (MQ=255)
         tACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccc                                                                                             <  2:3065227/101‑2 (MQ=255)
         tACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccc                                                                                             >  1:3065607/1‑100 (MQ=255)
           cATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccc                                                                                              >  2:3065609/1‑98 (MQ=255)
           cATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACACAGTG‑cccccccccag                                                                                           <  1:3065155/101‑4 (MQ=255)
              tCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCa                                                                                         >  2:3065613/1‑100 (MQ=255)
               ccccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccca                                                                                            <  2:3065220/96‑3 (MQ=255)
                 ccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATaa                                                                                      <  2:3065212/101‑1 (MQ=255)
                 ccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccc                                                                                               >  1:3065617/1‑91 (MQ=255)
                  cGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAAc                                                                                     <  1:3065252/101‑1 (MQ=255)
                   ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccca                                                                                            <  2:3065245/92‑3 (MQ=255)
                       tCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCat                                                                               >  2:3065621/1‑101 (MQ=255)
                        ccccccTTCGCTAATCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCat                                                                               >  1:3065624/1‑101 (MQ=255)
                         cccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccc                                                                                              >  1:3065627/1‑84 (MQ=255)
                          ccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCatat                                                                             <  2:3065437/101‑1 (MQ=255)
                          ccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCatat                                                                             <  1:3065421/101‑1 (MQ=255)
                          ccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCatat                                                                             <  2:3065189/101‑1 (MQ=255)
                           cccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATc                                                                            <  1:3065275/101‑1 (MQ=255)
                            ccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCAc                                                                          >  2:3065629/1‑101 (MQ=255)
                               tCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCt                                                                        <  1:3065182/101‑1 (MQ=255)
                                ctcTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTg                                                                       >  1:3065633/3‑101 (MQ=255)
                                  cTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCa                                                                     >  2:3065638/1‑101 (MQ=255)
                                  cTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                                    >  1:3065637/1‑101 (MQ=255)
                                   tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccc                                                                                              >  1:3065641/1‑74 (MQ=255)
                                    aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccc                                                                                             >  1:3065647/1‑73 (MQ=255)
                                    aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAgtg                                                                  >  1:3065643/1‑101 (MQ=255)
                                      tCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtgg                                                                 >  1:3065650/1‑101 (MQ=255)
                                         tttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGa                                                              <  1:3065195/101‑1 (MQ=255)
                                         tttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGAc                                                             <  1:3065101/101‑1 (MQ=255)
                                         tttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGAc                                                             >  1:3065654/1‑101 (MQ=255)
                                            gacccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                                    <  2:3065203/89‑1 (MQ=255)
                                             aaCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAtt                                                          >  1:3065657/1‑101 (MQ=255)
                                             aaCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAttt                                                         <  2:3065268/101‑1 (MQ=255)
                                              aCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAttt                                                         <  2:3065328/101‑1 (MQ=255)
                                              aCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTc                                                        <  2:3065342/101‑1 (MQ=255)
                                               ccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCatat                                                                             <  2:3065102/80‑1 (MQ=255)
                                               ccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                                       >  2:3065659/1‑101 (MQ=255)
                                                cGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCata                                                                              <  2:3065281/77‑1 (MQ=255)
                                                  tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGa                                                    >  1:3065663/1‑101 (MQ=255)
                                                   gccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAc                                                   >  2:3065665/1‑101 (MQ=255)
                                                      gccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACAt                                                 >  2:3065669/1‑101 (MQ=255)
                                                       ccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                                >  2:3065671/1‑101 (MQ=255)
                                                        cATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcac                                              >  1:3065673/1‑101 (MQ=255)
                                                          ttCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcaca                                             >  2:3065676/1‑101 (MQ=255)
                                                            ctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACAcc                                           <  2:3065402/101‑1 (MQ=255)
                                                              ccAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGc                                         >  2:3065680/1‑101 (MQ=255)
                                                              ccAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCt                                        <  1:3065335/101‑1 (MQ=255)
                                                               cAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCt                                        >  1:3065683/1‑101 (MQ=255)
                                                                 ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGt                                      >  2:3065688/1‑101 (MQ=255)
                                                                 ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                     >  1:3065687/1‑101 (MQ=255)
                                                                 ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                     <  1:3065551/101‑1 (MQ=255)
                                                                 ggTTCCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGc                                                                      <  2:3065234/68‑1 (MQ=255)
                                                                  gTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                     <  2:3065411/101‑1 (MQ=255)
                                                                    tGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTAt                                   >  1:3065694/1‑101 (MQ=255)
                                                                    tGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                                  >  2:3065693/1‑101 (MQ=255)
                                                                     gccgccGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                                  >  2:3065699/1‑101 (MQ=255)
                                                                     gccgccGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATgg                                 <  1:3065395/101‑1 (MQ=255)
                                                                      ccgccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTccc                                                      <  1:3064887/79‑1 (MQ=255)
                                                                       cgccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGat                               <  1:3065349/101‑1 (MQ=255)
                                                                       cgccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGat                               >  1:3065701/1‑101 (MQ=255)
                                                                        gccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGata                              >  2:3065704/1‑101 (MQ=255)
                                                                        gccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGata                              >  1:3065705/1‑101 (MQ=255)
                                                                         ccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGatat                             <  2:3065455/101‑1 (MQ=255)
                                                                          cgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCa                                                                     >  1:3065709/1‑61 (MQ=255)
                                                                           gATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATAtt                            >  2:3065714/1‑101 (MQ=255)
                                                                           gATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATAtt                            >  1:3065713/1‑101 (MQ=255)
                                                                            aTAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTa                           <  1:3065071/101‑1 (MQ=255)
                                                                            aTAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTaa                          <  2:3065368/101‑1 (MQ=255)
                                                                             tAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAt                         <  1:3065143/101‑1 (MQ=255)
                                                                             tAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAt                         <  2:3065198/101‑1 (MQ=255)
                                                                              aGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAt                         <  1:3065286/101‑1 (MQ=255)
                                                                              aGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtg                        <  1:3065294/101‑1 (MQ=255)
                                                                               gTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtgt                       <  1:3065174/101‑1 (MQ=255)
                                                                                tGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtgt                       >  1:3065717/1‑101 (MQ=255)
                                                                                 ggTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtgtg                      <  1:3065307/101‑1 (MQ=255)
                                                                                  gTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGata                              <  2:3065168/91‑1 (MQ=255)
                                                                                    tttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg                   <  2:3065353/101‑1 (MQ=255)
                                                                                    tttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg                   <  2:3065311/101‑1 (MQ=255)
                                                                                    tttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg                   <  1:3065133/101‑1 (MQ=255)
                                                                                    tttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg                   <  2:3065073/101‑1 (MQ=255)
                                                                                     ttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGa                  >  2:3065721/1‑101 (MQ=255)
                                                                                       aGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACt                >  1:3065723/1‑101 (MQ=255)
                                                                                        gCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTg              <  1:3065378/101‑1 (MQ=255)
                                                                                           acacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc            >  2:3065727/1‑101 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtg                                                                  <  1:3065367/46‑1 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACAcc                                           <  2:3065389/69‑1 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtgtg                      <  2:3065297/90‑1 (MQ=255)
                                                                                             acGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCat          >  2:3065730/1‑101 (MQ=255)
                                                                                              cGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCata         >  1:3065732/1‑101 (MQ=255)
                                                                                               gCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCatat        >  1:3065736/1‑101 (MQ=255)
                                                                                                cAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATAtt       <  1:3065264/101‑1 (MQ=255)
                                                                                                   tG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATg   <  1:3065237/101‑1 (MQ=255)
                                                                                                    g‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATgg  >  1:3065740/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
AACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGG  >  NC_000913/1594024‑1594224

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