Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,0230GA41.1% 117.4 / 197.0 180G320G (GGC→GGTompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base G (60/46);  new base A (25/49);  total (85/95)
Fisher's exact test for biased strand distribution p-value = 3.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 183 total aligned reads displayed.

GGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCG  >  NC_000913/985924‑986122
                                                                                                   |                                                                                                   
ggTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTgcg                                                                                                     >  1:2589013/1‑100 (MQ=255)
   gtcgccTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCAccc                                                                                                  >  1:1483744/1‑100 (MQ=255)
     ngtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCAc                                                                                                >  1:393115/2‑100 (MQ=255)
     cgtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATCGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCAc                                                                                                >  1:653616/1‑100 (MQ=255)
     cgtcTGAACCTACGCACAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCAc                                                                                                >  2:1885642/1‑100 (MQ=255)
      gtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACt                                                                                               <  2:862470/100‑1 (MQ=255)
      gtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACt                                                                                               <  2:2378568/100‑1 (MQ=255)
      gtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACt                                                                                               <  2:372728/100‑1 (MQ=255)
      gtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACt                                                                                               >  1:2247319/1‑100 (MQ=255)
       tcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCAACtt                                                                                              >  1:787621/1‑100 (MQ=255)
       tatGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACtt                                                                                              >  2:2408751/3‑100 (MQ=255)
         tGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCa                                                                                            >  2:1744345/1‑100 (MQ=255)
            aCCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAag                                                                                         <  2:2137546/100‑1 (MQ=255)
              cTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAagta                                                                                       <  2:911713/100‑1 (MQ=255)
               tACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAagtag                                                                                      <  2:1430337/100‑1 (MQ=255)
               tACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAagtag                                                                                      <  2:687640/100‑1 (MQ=255)
               caCGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATCGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACGCACTTCAAagtag                                                                                      >  1:2378033/2‑100 (MQ=37)
               atctcccAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAagtag                                                                                      <  2:2444701/96‑1 (MQ=255)
                aCGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAagtagt                                                                                     <  2:1098380/100‑1 (MQ=255)
                aCGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAagtagt                                                                                     <  2:1885582/100‑1 (MQ=255)
                 cGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGtt                                                                                    >  1:2588388/1‑100 (MQ=255)
                  gCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTc                                                                                   <  2:75239/100‑1 (MQ=255)
                   cccAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCa                                                                                  >  1:769154/1‑100 (MQ=255)
                   cccAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCa                                                                                  >  1:911788/1‑100 (MQ=255)
                   cccAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCa                                                                                  >  1:2166450/1‑100 (MQ=255)
                    ccAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCAc                                                                                 >  1:1360649/1‑100 (MQ=255)
                     cAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTGTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCAcc                                                                                >  1:879998/1‑100 (MQ=255)
                          tgttgtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATc                                                                           <  1:1790659/100‑1 (MQ=255)
                           gttgtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCa                                                                          >  1:1709019/1‑100 (MQ=255)
                            ttgtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCan                                                                         <  1:458441/100‑2 (MQ=255)
                              gtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACa                                                                       <  1:1733122/100‑1 (MQ=255)
                              gtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACa                                                                       >  1:1593931/1‑100 (MQ=255)
                              gtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACa                                                                       >  1:2183837/1‑100 (MQ=255)
                              gtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACa                                                                       <  2:1463256/100‑1 (MQ=255)
                              gtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACa                                                                       >  1:555508/1‑100 (MQ=255)
                              gtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACa                                                                       <  2:844342/100‑1 (MQ=255)
                               tCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACAt                                                                      >  1:2040096/1‑100 (MQ=255)
                                cAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTTTTTTTTTAGATTGCGCCCACTTCAAAGTAGTTCACGAGTTCAACATc                                                                     >  2:1691715/1‑100 (MQ=25)
                                    aTCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCg                                                                 <  2:2048307/100‑1 (MQ=255)
                                    aTCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCg                                                                 >  1:1854069/1‑100 (MQ=255)
                                    aTCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCg                                                                 <  2:301562/100‑1 (MQ=255)
                                     tCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGa                                                                <  2:351433/100‑1 (MQ=255)
                                      cGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGAt                                                               <  2:1254517/100‑1 (MQ=255)
                                       gATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTTAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACGAGATCAACATCACCGATa                                                              >  2:1231840/1‑100 (MQ=255)
                                        aTCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATAc                                                             <  2:1547474/100‑1 (MQ=255)
                                         tCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATAcc                                                            <  1:1284990/100‑1 (MQ=255)
                                              ttGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTa                                                       <  1:106612/100‑1 (MQ=255)
                                               tgatgaTGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTAc                                                      <  2:2095144/100‑1 (MQ=255)
                                                 atgatgTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGt                                                    >  2:1594759/1‑100 (MQ=255)
                                                 atgatgTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGt                                                    <  1:280864/100‑1 (MQ=255)
                                                    atgTAGTCAACATAGGTGTACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCtt                                                 <  2:636213/100‑1 (MQ=255)
                                                     tgTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCttt                                                <  2:1572023/100‑1 (MQ=255)
                                                      gTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTc                                               >  1:1369350/1‑100 (MQ=255)
                                                       tAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTcc                                              <  2:1349935/100‑2 (MQ=255)
                                                       tAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCg                                              >  1:1990415/1‑100 (MQ=255)
                                                        aGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGc                                             >  1:532356/1‑100 (MQ=255)
                                                         gTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCt                                            >  1:119163/1‑100 (MQ=255)
                                                              cATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGa                                       >  1:363407/1‑100 (MQ=255)
                                                                tAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGAtt                                     >  1:1011894/1‑100 (MQ=255)
                                                                  ggTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTg                                   >  1:2257257/1‑100 (MQ=255)
                                                                  ggTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTg                                   <  2:2206179/100‑1 (MQ=255)
                                                                  ggTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTCTAGATTTg                                   >  1:2522087/1‑100 (MQ=255)
                                                                   gTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTgg                                  <  1:1945088/100‑1 (MQ=255)
                                                                   gTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTgg                                  <  2:1794211/100‑1 (MQ=255)
                                                                   gTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTGAGATTTgg                                  <  2:2489106/100‑1 (MQ=255)
                                                                   gTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTgg                                  <  2:1268038/100‑1 (MQ=255)
                                                                   gTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTgg                                  <  2:442709/100‑1 (MQ=255)
                                                                   gTGGACATGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTgg                                  <  2:1200551/100‑1 (MQ=255)
                                                                    tGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGgt                                 >  1:132723/1‑100 (MQ=255)
                                                                     ggACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGgtg                                <  2:1942202/100‑1 (MQ=255)
                                                                      gACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTGTGCGTCTTTCGCTTTCGATTTGgtgt                               >  2:2570488/1‑100 (MQ=255)
                                                                      gACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGgtgt                               <  1:2431513/100‑1 (MQ=255)
                                                                       aCATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTa                              <  2:2575853/100‑1 (MQ=255)
                                                                         aTGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAg                            <  1:1622363/100‑1 (MQ=255)
                                                                         aTGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAg                            <  2:1308811/100‑1 (MQ=255)
                                                                          tGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGc                           <  2:1645935/100‑1 (MQ=255)
                                                                          tGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGc                           <  2:438724/100‑1 (MQ=255)
                                                                          tGTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGc                           <  2:2541083/100‑1 (MQ=255)
                                                                           gTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCg                          <  2:1982256/100‑1 (MQ=255)
                                                                           gTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCg                          <  2:635241/100‑1 (MQ=255)
                                                                           gTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCg                          <  2:2149388/100‑1 (MQ=255)
                                                                           gTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCg                          <  2:668330/100‑1 (MQ=255)
                                                                           gTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCg                          >  1:1234926/1‑100 (MQ=255)
                                                                           gTTTTTGTTGAAGTAGTAGGTTGCACCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTGCTACGTCTTTCGCTTTAGATTTGGTGTAAGCg                          <  2:2332978/100‑1 (MQ=255)
                                                                             ttttttttAATTATTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTTGTGTAAGCGAt                        >  2:2235490/1‑100 (MQ=255)
                                                                             ttttGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGAt                        <  2:1829414/100‑1 (MQ=255)
                                                                             ttttGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGAt                        <  1:50131/100‑1 (MQ=255)
                                                                                tgttgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATgga                     >  2:558877/1‑100 (MQ=255)
                                                                                 gttgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATggac                    >  1:116674/1‑100 (MQ=255)
                                                                                 gttgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTCGGTGTAAGCGATggac                    >  1:2177003/1‑100 (MQ=255)
                                                                                 gttgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCGTCTACGTCTTTCACTTTAGAGTTGGTGTAAGCGATggac                    >  1:2473800/1‑100 (MQ=255)
                                                                                  ttgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATggacg                   <  2:438672/100‑1 (MQ=255)
                                                                                   tgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATggacgg                  >  2:2375701/1‑100 (MQ=255)
                                                                                   tgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATggacgg                  >  2:1401044/1‑100 (MQ=255)
                                                                                            agGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGAcc         <  2:1519503/100‑1 (MQ=255)
                                                                                             ggTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCg        <  2:2405779/100‑1 (MQ=255)
                                                                                             ggTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCg        <  1:299595/100‑1 (MQ=255)
                                                                                              gTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGa       >  1:694540/1‑100 (MQ=255)
                                                                                               ttGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGaa      >  1:1133639/1‑100 (MQ=255)
                                                                                                   gcCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCg  >  2:377763/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCG  >  NC_000913/985924‑986122

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