Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5660TC28.2% 449.1 / 231.5 309E139E (GAA→GAG) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base T (119/102);  new base C (49/38);  total (169/140)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 312 total aligned reads displayed.

GCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAAT  >  NC_000913/986467‑986665
                                                                                                   |                                                                                                   
  cATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGTCGTATCACCACCAAATTc                                                                                                   >  2:4272764/1‑100 (MQ=255)
  cATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTc                                                                                                   >  2:4131562/1‑100 (MQ=255)
   aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCt                                                                                                  >  2:844602/1‑100 (MQ=255)
   aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCt                                                                                                  <  1:1541351/100‑1 (MQ=255)
     cAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAGCACCACCAAATTCTgg                                                                                                >  1:4761189/1‑100 (MQ=255)
      aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                               <  1:367203/100‑1 (MQ=255)
        ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgcag                                                                                             >  1:3896564/1‑100 (MQ=255)
        ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGACGCAAACCCGACCAACGAAGAAGTCATAGCTGTATGAAGTAACACCCCCAAATCCTAgcag                                                                                             >  1:2573308/1‑100 (MQ=255)
         cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                            <  1:194513/100‑1 (MQ=255)
           gACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCat                                                                                          >  1:4538096/1‑100 (MQ=255)
               aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                      >  2:3510052/1‑100 (MQ=255)
               aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                      >  2:4927616/1‑100 (MQ=255)
                aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGGCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTATGGCAGCATATCgg                                                                                     <  2:954520/100‑1 (MQ=255)
                aagaagTTGGAGTTACGATAGGGAGCAACGCCGCCAACACGACCGACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                     <  1:3086033/100‑1 (MQ=255)
                  gaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtg                                                                                   <  1:5026902/100‑1 (MQ=255)
                  gaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtg                                                                                   >  2:4233343/1‑100 (MQ=255)
                    agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTa                                                                                 <  2:3047526/100‑1 (MQ=255)
                    agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTa                                                                                 >  1:2805213/1‑100 (MQ=255)
                     gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTaa                                                                                >  1:3172081/1‑100 (MQ=255)
                        ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAccc                                                                             >  1:2955721/1‑100 (MQ=38)
                         gAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCa                                                                            >  1:2326320/1‑100 (MQ=255)
                           gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGt                                                                          >  1:1810530/1‑100 (MQ=38)
                             tACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGc                                                                        <  2:600434/100‑1 (MQ=255)
                              aCGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGca                                                                       >  2:3273869/1‑100 (MQ=255)
                               cGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcat                                                                      >  2:1517642/1‑100 (MQ=255)
                                gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatc                                                                     <  2:3854328/100‑1 (MQ=38)
                                gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                                     <  2:4326143/100‑1 (MQ=255)
                                 aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatca                                                                    <  1:3156801/100‑1 (MQ=255)
                                    ggTAGCAACGCCGCGAAGATGACCAACGTTGAAATCTTCGCTGTATGCAGTATCACTATCAAATTCTGCCAGCTTGTCGCCGTGCACCTTTGCATCATaa                                                                 >  2:2565089/1‑100 (MQ=255)
                                    ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                                 >  2:1875311/1‑100 (MQ=255)
                                     gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                                >  1:340571/1‑100 (MQ=38)
                                     gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                                <  2:2404935/100‑1 (MQ=38)
                                      tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaan                                                               <  1:4115681/100‑2 (MQ=255)
                                       aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                              >  1:833846/1‑100 (MQ=38)
                                       aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                              <  1:2030166/100‑1 (MQ=255)
                                       aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                              >  2:4524215/1‑100 (MQ=255)
                                         cAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGGCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCAGCATAAACcac                                                            <  2:3660210/100‑1 (MQ=255)
                                         cAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcac                                                            <  2:2885648/100‑1 (MQ=255)
                                          caCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcaca                                                           >  1:940589/2‑100 (MQ=255)
                                           ccgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                          >  1:2412240/2‑100 (MQ=38)
                                           aCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                          <  1:4598733/100‑1 (MQ=255)
                                            cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                         >  1:3758925/1‑100 (MQ=38)
                                             gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                        >  1:1299892/1‑100 (MQ=38)
                                              ccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGt                                                       <  2:2935316/100‑1 (MQ=255)
                                               cgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTa                                                      >  2:890508/1‑100 (MQ=255)
                                                gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                     <  2:820721/100‑1 (MQ=38)
                                                gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                     <  2:2076007/100‑1 (MQ=38)
                                                gccAACACGACCAACGAAGAACTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                     <  2:1800066/100‑1 (MQ=255)
                                                       cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                              <  2:712102/100‑1 (MQ=255)
                                                       cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                              <  1:3375535/100‑1 (MQ=255)
                                                       cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                              <  1:2196404/100‑1 (MQ=255)
                                                       cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                              <  1:1350777/100‑1 (MQ=255)
                                                       cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACAACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                              <  2:4613366/100‑1 (MQ=255)
                                                         aCCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                                            <  2:3388059/100‑1 (MQ=38)
                                                            aaCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTaa                                         >  2:1595300/1‑100 (MQ=255)
                                                            aaCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTaa                                         >  1:1246088/1‑100 (MQ=255)
                                                             aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAAt                                        <  2:925711/100‑1 (MQ=255)
                                                               gaagaaGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCg                                      <  1:989443/100‑1 (MQ=255)
                                                                 agaagTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                                    <  2:2394013/100‑1 (MQ=38)
                                                                 agaacTCATCGCTGTATGCAGTATCACCACCAAATTGTGGCAGCATATCAGTGTAACCCAGTGCATCATGAACCACACCGTAGTTACGGCCGTGATCGaa                                    >  2:4773298/1‑100 (MQ=255)
                                                                  gaagTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                                   <  2:3012226/100‑1 (MQ=38)
                                                                  gaagTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                                   <  2:4020638/100‑1 (MQ=38)
                                                                  gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                                   <  2:1168085/100‑1 (MQ=255)
                                                                   aagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCGTCATAAACCACACCGTAGTTACGGCCGTAATCGAAAg                                  >  1:4006201/1‑100 (MQ=255)
                                                                   aagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAg                                  <  2:4365734/100‑1 (MQ=255)
                                                                      tCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTGACGGCCGTAATCGAAAGAAc                               <  2:2302065/100‑1 (MQ=255)
                                                                       cATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                              >  1:3341425/1‑100 (MQ=38)
                                                                       cATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                              <  1:1042725/100‑1 (MQ=255)
                                                                        aTCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCa                             >  2:5016512/1‑100 (MQ=255)
                                                                          cGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATGGCGGTGTAACCCAGTGCATCATAAACCACCCCGCAGTTACGGCCGTAATCGAAAGAACCAAc                           <  2:4896134/100‑1 (MQ=255)
                                                                          agCTGTAGGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAAc                           <  2:1071710/99‑1 (MQ=38)
                                                                           gCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg                          >  1:2067952/1‑100 (MQ=38)
                                                                             tGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTc                        <  1:4838725/100‑1 (MQ=255)
                                                                                gtGCAGTATCACCACCAAATTCTGGCAGCATATCGTTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc                     >  1:1900629/2‑100 (MQ=255)
                                                                                 tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                    >  1:2335821/1‑100 (MQ=255)
                                                                                 tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                    <  2:1680362/100‑1 (MQ=255)
                                                                                 tGCAGTATCACCACCAAATCCTGGCAGCATATCGGTGTAACCGAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                    >  1:4392824/1‑100 (MQ=255)
                                                                                  gCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGGACCAACGTCAGCGt                   >  2:2125084/1‑100 (MQ=255)
                                                                                  gCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGt                   >  1:3512374/1‑100 (MQ=255)
                                                                                   cAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTa                  >  2:3160147/1‑100 (MQ=255)
                                                                                   cAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTa                  <  1:2134508/100‑1 (MQ=255)
                                                                                        tCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTaa             <  1:1757136/100‑1 (MQ=255)
                                                                                         caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTGCGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg            >  1:4927061/1‑100 (MQ=255)
                                                                                         caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg            >  1:3324572/1‑100 (MQ=255)
                                                                                           ccgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc          >  1:795725/4‑100 (MQ=255)
                                                                                           ccgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc          <  2:1203254/97‑1 (MQ=255)
                                                                                           ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc          <  1:4692840/100‑1 (MQ=255)
                                                                                           ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc          >  1:4661730/1‑100 (MQ=255)
                                                                                            cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc         >  1:1118294/3‑100 (MQ=255)
                                                                                            cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc         >  1:742852/3‑100 (MQ=255)
                                                                                            cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc         >  1:4344761/3‑100 (MQ=255)
                                                                                            cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc         >  1:4241355/3‑100 (MQ=255)
                                                                                            cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc         >  1:3024905/3‑100 (MQ=255)
                                                                                            caccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc         <  1:58475/100‑1 (MQ=255)
                                                                                             gccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGCAAGAACAAACGTCAGCGTATGTAAGAcca        >  1:5105524/2‑99 (MQ=255)
                                                                                             accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc        >  1:4153287/1‑100 (MQ=255)
                                                                                             accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc        >  1:3271020/1‑100 (MQ=255)
                                                                                                aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcgcg     <  1:3469256/100‑1 (MQ=255)
                                                                                                aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcgcg     <  1:2433230/100‑1 (MQ=255)
                                                                                                aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAAAATATTTAAGAcaagcg     <  2:1380654/100‑6 (MQ=39)
                                                                                                   |                                                                                                   
GCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAAT  >  NC_000913/986467‑986665

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