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.1% 963.9 / ‑0.4 233pseudogene (3963/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (128/103);  total (128/105)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

AACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATT  >  NC_000913/1594024‑1594219
                                                                                                     |                                                                                               
aaCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTg                                                                                                  >  2:2679852/1‑101 (MQ=255)
   ggTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccc                                                                                              >  2:2679856/1‑101 (MQ=255)
      cAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccc                                                                                           >  1:2679859/1‑101 (MQ=255)
       aGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccc                                                                                          >  1:2679862/1‑101 (MQ=255)
       aGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccc                                                                                          <  2:2679467/101‑1 (MQ=255)
         tACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccc                                                                                        >  1:2679866/1‑100 (MQ=255)
          aCATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccca                                                                                       >  2:2679869/1‑99 (MQ=255)
           cATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccccag                                                                                      >  1:2679873/1‑98 (MQ=255)
            aTTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccccagg                                                                                     >  2:2679874/1‑97 (MQ=255)
            aTTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccccagg                                                                                     >  1:2679876/1‑97 (MQ=255)
              tCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCa                                                                                    >  2:2679881/1‑101 (MQ=255)
               ccccGGTTCCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑ccccccccca                                                                                       <  2:2679568/96‑3 (MQ=255)
                cccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATa                                                                                  >  1:2679884/1‑101 (MQ=255)
                     tttCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGtt                                                                             <  2:2679415/101‑1 (MQ=255)
                     tttCCCCCCTCCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGtt                                                                             >  1:2679891/1‑101 (MQ=255)
                        ccccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCat                                                                          >  1:2679894/1‑101 (MQ=255)
                         cccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCata                                                                         <  1:2679497/101‑1 (MQ=255)
                              ttCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCAcc                                                                    >  2:2679897/1‑101 (MQ=255)
                              ttCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCAcc                                                                    >  1:2679899/1‑101 (MQ=255)
                               tCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGTGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCt                                                                   <  1:2679504/101‑1 (MQ=255)
                                  cTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCa                                                                >  1:2679903/1‑101 (MQ=255)
                                   tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCGGCATAACGTTCATATCACCTGCa                                                                >  2:2679908/1‑100 (MQ=255)
                                   tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                               >  1:2679907/1‑101 (MQ=255)
                                      tCATTTGAACCGGGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtgg                                                            <  2:2679556/101‑1 (MQ=255)
                                       cATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtggt                                                           >  1:2679911/1‑101 (MQ=255)
                                        aTTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtggtg                                                          >  1:2679915/1‑101 (MQ=255)
                                         tttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                               <  2:2679393/95‑1 (MQ=255)
                                         tttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGa                                                         <  1:2679334/101‑1 (MQ=255)
                                          ttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGAc                                                        >  1:2679919/1‑101 (MQ=255)
                                          ttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGAc                                                        >  2:2679918/1‑101 (MQ=255)
                                             aaCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAtt                                                     >  1:2679922/1‑101 (MQ=255)
                                             aaCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAtt                                                     <  1:2679399/101‑1 (MQ=255)
                                               ccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTc                                                   >  2:2679925/1‑101 (MQ=255)
                                                cGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                                  >  1:2679926/1‑101 (MQ=255)
                                                 gTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTccc                                                 >  2:2679931/1‑101 (MQ=255)
                                                  tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCg                                                >  1:2679933/1‑101 (MQ=255)
                                                  tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCg                                                <  1:2679668/101‑1 (MQ=255)
                                                   gccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTg                                                                                                  >  2:2679937/1‑50 (MQ=255)
                                                   gccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGa                                               >  1:2679936/1‑101 (MQ=255)
                                                     cgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCa                                                                <  1:2679578/82‑1 (MQ=255)
                                                     cgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACa                                             <  2:2679610/101‑1 (MQ=255)
                                                      gccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACAt                                            <  1:2679539/101‑1 (MQ=255)
                                                        cATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtgg                                                            >  2:2679944/1‑83 (MQ=255)
                                                        cATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATca                                          >  1:2679941/1‑101 (MQ=255)
                                                          ttCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATa                                                                                  <  2:2679491/59‑1 (MQ=255)
                                                          ttCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcaca                                        >  1:2679947/1‑101 (MQ=255)
                                                           tctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                       <  1:2679599/101‑1 (MQ=255)
                                                           tctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                       <  2:2679524/101‑1 (MQ=255)
                                                            ctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCAc                                                                     <  1:2679184/70‑1 (MQ=255)
                                                            ctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACAcc                                      <  1:2679625/101‑1 (MQ=255)
                                                            ctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACAcc                                      >  2:2679949/1‑101 (MQ=255)
                                                             tcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAt                                                      >  2:2679953/1‑84 (MQ=255)
                                                               cAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCt                                   >  1:2679956/1‑101 (MQ=255)
                                                               cAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCt                                   <  2:2679600/101‑1 (MQ=255)
                                                               cAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCt                                   <  1:2679742/101‑1 (MQ=255)
                                                                aGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTCCAGTGGTGACATTTccc                                                 >  1:2679959/1‑86 (MQ=255)
                                                                 ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                                  <  2:2679199/84‑1 (MQ=255)
                                                                 ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGt                                 <  2:2679403/101‑1 (MQ=255)
                                                                  gTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                >  1:2679962/1‑101 (MQ=255)
                                                                   ttGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTa                               <  1:2679528/101‑1 (MQ=255)
                                                                   ttGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTa                               <  1:2679463/101‑1 (MQ=255)
                                                                    tGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtggtg                                                          >  2:2679965/1‑73 (MQ=255)
                                                                    tGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTAt                              <  1:2679570/101‑1 (MQ=255)
                                                                     gccgccGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                             >  1:2679967/1‑101 (MQ=255)
                                                                      ccgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCAt                                                                                   <  2:2679587/46‑1 (MQ=255)
                                                                      ccgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCatat                                                                        <  1:2679450/57‑1 (MQ=255)
                                                                      ccgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATgg                            >  2:2679969/1‑101 (MQ=255)
                                                                       cgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGa                           >  2:2679972/1‑101 (MQ=255)
                                                                        gccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGat                          >  1:2679977/1‑101 (MQ=255)
                                                                         ccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGata                         <  2:2679483/101‑1 (MQ=255)
                                                                          cgATAGTTGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGatat                        <  2:2679623/101‑1 (MQ=255)
                                                                          cgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGatat                        >  2:2679979/1‑101 (MQ=255)
                                                                           gATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAgtgg                                                            >  2:2679985/1‑64 (MQ=255)
                                                                           gATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATAtt                       >  1:2679982/1‑101 (MQ=255)
                                                                            aTAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTa                      <  1:2679595/101‑1 (MQ=255)
                                                                            aTAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTa                      <  2:2679699/101‑1 (MQ=255)
                                                                             tAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTaa                     >  1:2679988/1‑101 (MQ=255)
                                                                               gTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtg                   >  2:2679990/1‑101 (MQ=255)
                                                                               gTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtg                   >  2:2679993/1‑101 (MQ=255)
                                                                                tGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtgt                  <  1:2679550/101‑1 (MQ=255)
                                                                                   ttttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGtt               <  1:2679305/101‑1 (MQ=255)
                                                                                   ttttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGtt               <  1:2679459/101‑1 (MQ=255)
                                                                                   ttttAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGtt               >  1:2679995/1‑101 (MQ=255)
                                                                                      tAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGAc            <  2:2679754/101‑1 (MQ=255)
                                                                                      tAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGAc            >  2:2679999/1‑101 (MQ=255)
                                                                                       aGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACt           >  1:2680000/1‑101 (MQ=255)
                                                                                       aGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACt           <  1:2679694/101‑1 (MQ=255)
                                                                                         cGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCg                                                >  2:2680003/1‑62 (MQ=255)
                                                                                           acacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc       >  1:2680008/1‑101 (MQ=255)
                                                                                           acacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc       <  1:2679446/101‑1 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTg                                                                  <  2:2679645/41‑1 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa      >  1:2680012/1‑101 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa      <  2:2679486/101‑1 (MQ=255)
                                                                                            cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa      >  2:2680010/1‑101 (MQ=255)
                                                                                             acGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCat     >  1:2680016/1‑101 (MQ=255)
                                                                                             acGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCat     <  2:2679540/101‑1 (MQ=255)
                                                                                             acGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa      <  1:2679720/100‑1 (MQ=255)
                                                                                              cGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtg                   <  2:2679679/85‑1 (MQ=255)
                                                                                               gCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCatat   <  1:2679638/101‑1 (MQ=255)
                                                                                               gCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCatat   >  2:2680019/1‑101 (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 < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑