Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,431,9630AG51.0% 72.8 / 131.6 96S11P (TCC→CCC) msrApeptide‑methionine (S)‑S‑oxide reductase MsrA
Reads supporting (aligned to +/- strand):  ref base A (7/40);  new base G (49/0);  total (56/40)
Fisher's exact test for biased strand distribution p-value = 3.68e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-15
Rejected as consensus: Frequency below/above cutoff threshold.

AGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTT  >  NZ_CP009273/4431874‑4432051
                                                                                         |                                                                                         
aGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGa                                                                                           <  1:1137756/90‑1 (MQ=255)
 ggTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGaa                                                                                          <  1:1717688/90‑1 (MQ=255)
      aTTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAg                                                                                     >  1:1173793/1‑90 (MQ=255)
        tGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGAt                                                                                   <  1:1048492/90‑1 (MQ=255)
        tGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGAt                                                                                   <  1:1459961/90‑1 (MQ=255)
          gTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAAACAGATGc                                                                                 >  2:1295417/1‑90 (MQ=255)
           tCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATggt                                                                                >  1:113525/1‑88 (MQ=255)
           tCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCt                                                                               >  2:1425069/1‑90 (MQ=255)
            cATTGAGTGACCGTTGACCGCATGCAGCGTGTCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAAAAGCtt                                                                              >  1:2150191/1‑90 (MQ=255)
            cATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCtt                                                                              >  2:1056791/1‑90 (MQ=255)
            cATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGAAGCtt                                                                              >  2:507702/1‑90 (MQ=255)
            cATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGAAGCtt                                                                              >  1:1175549/1‑90 (MQ=255)
              ttGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGAAGGtttt                                                                            >  1:1005320/1‑90 (MQ=255)
                gAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGACGGGGGAAACCAGATGCTTTTTa                                                                          <  1:1925283/90‑1 (MQ=255)
                 aGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTAt                                                                         <  2:166715/90‑1 (MQ=255)
                  gTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGTGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATc                                                                        <  1:1444563/90‑1 (MQ=255)
                  gTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATc                                                                        <  1:599288/90‑1 (MQ=255)
                   tGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGGAGCAGATGCtttttttta                                                                       >  1:468345/1‑86 (MQ=255)
                   tGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAGCCGCGGCtttttttta                                                                       >  1:1812591/1‑86 (MQ=255)
                   tGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCCAATGCTTTTTTTCa                                                                       >  2:131882/1‑90 (MQ=255)
                   tGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGGTGCCTTTTATCa                                                                       >  2:411652/1‑90 (MQ=255)
                   tGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCtttttttta                                                                       >  2:669571/1‑86 (MQ=255)
                   tGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATTCTTTTTATCa                                                                       >  2:2069403/1‑90 (MQ=255)
                      ccGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGTCGGGGGAAACCAGATGCTTTTTATCAAAt                                                                    <  2:1836041/90‑1 (MQ=255)
                      ccGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGTCGGGGGAAACCAGATGCTTTTTATCAAAt                                                                    <  2:1282705/90‑1 (MQ=255)
                      ccGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGAATCGGCGGGGGAAACCAGATGCTTTTTATCAAAt                                                                    <  1:1204878/90‑1 (MQ=255)
                       cGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGTTGGTTTTTTTTaaaaa                                                                   >  1:140912/1‑88 (MQ=255)
                       cGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGAAGGTTTTTTTTaaaaa                                                                   >  1:2021420/1‑88 (MQ=255)
                        gTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGGTTTTTTTCAAATaa                                                                  >  2:687861/1‑90 (MQ=255)
                        gTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCTTTTTTTTAAATaa                                                                  >  2:848242/1‑90 (MQ=255)
                          tGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTTAaaaaaaa                                                                >  2:533802/1‑89 (MQ=255)
                            aCCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAAAAGCTTTTTTTTAAATAAAATc                                                              >  2:1076567/1‑90 (MQ=255)
                               gCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCTTTTTTTCAAAAAAAATCATg                                                           >  2:211626/1‑90 (MQ=255)
                               gCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTCAAAAAAAATCATg                                                           >  2:544897/1‑90 (MQ=255)
                               gCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGAAGGTTTTTTTTTAATAAAATCATg                                                           >  1:1907763/1‑90 (MQ=255)
                               gCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAAAAGCTTTTTTTTAAAAAAAATCATg                                                           >  1:994577/1‑90 (MQ=255)
                                 aTGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGgt                                                         <  2:1993970/90‑1 (MQ=255)
                                  tGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGAt                                                                                   >  1:1810338/1‑64 (MQ=255)
                                  tGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGAt                                                                                   >  1:1282633/1‑64 (MQ=255)
                                  tGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGAAGGTTTTTTTTAAATAAAATCATgggg                                                        >  1:1785802/1‑88 (MQ=255)
                                  tGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGAt                                                                                   <  2:1282633/64‑1 (MQ=255)
                                  tGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGAt                                                                                   <  2:1810338/64‑1 (MQ=255)
                                     aGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGGTTTTTTTTAAATAAAATCATGGGGTCg                                                     >  1:1888348/1‑90 (MQ=255)
                                      gCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTATCAAATAAAATCCTGGGGTCGc                                                    >  2:1602037/1‑90 (MQ=255)
                                      gCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGc                                                    >  2:1501291/1‑90 (MQ=255)
                                      gCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGAGGCTTTTTTTCAAAAAAAATCATGGGGTCGc                                                    >  2:865577/1‑90 (MQ=255)
                                        gTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGctc                                                  <  1:1153154/90‑1 (MQ=255)
                                        gTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGctc                                                  <  1:1112764/90‑1 (MQ=255)
                                         tGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTATCAAATAAACTCATGGTGTTGcgct                                                 >  2:1913069/1‑90 (MQ=255)
                                          ggCTACGGGCATCGAGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTATCAAAAAAACTCATGGGGTTGCtctc                                                >  2:1835528/1‑90 (MQ=255)
                                           gCTACGGGCATAGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTcc                                               <  2:291879/90‑1 (MQ=255)
                                              aCGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGGTTTTTTTTAAAAAAAATCATGGGGTCGCTCTcccca                                            >  1:1164957/1‑88 (MQ=255)
                                              aCGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGa                                            >  2:2010296/1‑90 (MQ=255)
                                              aCGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTCAAAAAAAATCATGGGGGCGCTCTCCCGa                                            >  2:164373/1‑90 (MQ=255)
                                                gggCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCTTTTTTTTAAATAAACTCATGGGGTCGCTCTCCCGaaa                                          >  2:1688003/1‑90 (MQ=255)
                                                    aTCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGGTTTTTTTTAAAAAAATTCCTTGTGTCGCTCTCCCCGAAATct                                      >  2:1328693/1‑89 (MQ=255)
                                                    aTCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTTAAATAAAATCATGGGGTCGCTCTCCCCAAAAACg                                      >  2:524912/1‑90 (MQ=255)
                                                    aTCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGAAGGTTTTTTTTAAAAAAAATCATGGTGTCGCTCTCCCGAAAAACg                                      >  1:1886056/1‑90 (MQ=255)
                                                         ggTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCAt                                 <  2:898485/90‑1 (MQ=255)
                                                          gtgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCAtt                                <  2:970072/90‑1 (MQ=255)
                                                          gtgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCAtt                                <  1:1152367/90‑1 (MQ=255)
                                                          gtgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCAtt                                <  1:1835528/90‑1 (MQ=255)
                                                          gtgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCAtt                                <  1:175011/90‑1 (MQ=255)
                                                          gtgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCAtt                                <  1:131882/90‑1 (MQ=255)
                                                          gtgtGACGTCCAGGCAGGGCATCGGCGGGGGAAACCATATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCAtt                                <  1:2020373/90‑1 (MQ=255)
                                                           tgtTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCTTTTTTTCAAATAAACTCATGGGGTCGCTCTCCCGAAAATCGGTCATTg                               >  2:1123426/1‑90 (MQ=255)
                                                           tgtTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCTTTTTTTCAAATAAAATCATGGTGTCGCTCTCCCGAAAAACGGGCATTg                               >  2:175011/1‑90 (MQ=255)
                                                           tgtTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTCAAATAAACTCATGGGGTCGCTCTCCCGAAAATCGGTCATTg                               >  2:273208/1‑90 (MQ=255)
                                                           tgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATTCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTg                               <  1:101253/90‑1 (MQ=255)
                                                           tgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCAtt                                <  2:1999346/90‑1 (MQ=255)
                                                           tgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTg                               <  2:169928/90‑1 (MQ=255)
                                                           tgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTg                               <  2:2009999/90‑1 (MQ=255)
                                                           tgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTg                               <  2:1823099/90‑1 (MQ=255)
                                                           tgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTg                               <  2:1549514/90‑1 (MQ=255)
                                                            gtTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGAAGGTTTTTTTTAAAAAAAATCATGGGGTCGCTCTCCCGAAAAACGGGCAtttg                              >  1:239948/1‑88 (MQ=255)
                                                               aCGTCCAGGCAGGGCATCGGCGGGGGAAAACAGATGCTTTTTTTCAAAAAAACTCATGGGGTCGCTCTCCCGAAAAACGGGCATTggggg                           >  2:1441156/1‑89 (MQ=255)
                                                                 gTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTa                         <  1:1018405/90‑1 (MQ=255)
                                                                 gTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTa                         <  2:85743/90‑1 (MQ=255)
                                                                 gTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTa                         <  1:1688003/90‑1 (MQ=255)
                                                                  tCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTTTAAAAAAAATCATGGGGTCGCTCTCCCCAAAAACGGTCATTTGGGTTa                         >  1:1999346/1‑90 (MQ=255)
                                                                    caggcaggGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAgg                      <  1:1425069/90‑1 (MQ=255)
                                                                     aggcaggGCATCGGCGGGGGGAAACAGATGCTTTTTTTTAAAAAAAATCATGGGGGGGCTCTCCCCGAAAACGGGCGTTTGGGTTAAggg                     >  1:912600/1‑89 (MQ=255)
                                                                     aggcaggGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGt                     <  1:848242/90‑1 (MQ=255)
                                                                      ggcaggGCATCGGCGGGGGGAACCAGAAGCTTTTTTTTAAAAAAAATCATGGGGTCGCTCTCCCCAAAAACGGTCATTTGGGGTAAgggt                    >  1:2002173/1‑88 (MQ=255)
                                                                      ggcaggGCATCGGCGGGGGGAAACAGATTCTTTTTTTTAAATAAAAGCATGGGGGGGCTCTCCCGAAAAACGGGCATTTGGGGTAAAGtt                    >  1:59583/1‑90 (MQ=255)
                                                                      ggcaggGCATCGGCGGGGGGAAACAGATGCTTTTTTTTAAATAAAATCATGGGGTGGCTCTCCCGAAAAACGGGCATTGGGGGTAAgggt                    >  1:747536/1‑88 (MQ=255)
                                                                       gcaggGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTg                   <  2:994577/90‑1 (MQ=255)
                                                                          gggCATCGGCGGGGGGAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTaa                >  2:599941/1‑90 (MQ=255)
                                                                          gggCATCGGCGGGGGGAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTaa                >  2:1800460/1‑90 (MQ=255)
                                                                          gggCATCGGCGGGGGGAAACAGATGCTTTTTTTTAAAAAAAATCATGGGGTCGCTCTCCCCAAAAACGGGCATTTGGGGTAAGGCTGtct                >  1:1597727/1‑88 (MQ=255)
                                                                          gggCATCGGCGGGGGGAAACAGATGCTTTTTTTCAAAAAAAATCATGGGGGCGCTCTCCCCAAAAACGGGCATTTGGGGTAAGGTTGTaa                >  1:1661538/1‑90 (MQ=255)
                                                                          gggCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAAATCATGGGGTCGCTCTCCCCAAAATCGGTCATTGGGGTTAAGGTTTTaa                >  2:2101618/1‑90 (MQ=255)
                                                                            gCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACa              <  2:1395014/90‑1 (MQ=255)
                                                                               tCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAaga           <  1:1328693/90‑1 (MQ=255)
                                                                                  gcggGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGc        <  1:871962/90‑1 (MQ=255)
                                                                                    ggggggAACCAGATGCTTTTTTTTAAAAAAAATCATGGTGTCGCTCTCCCGAAAAACGGGCATTGGGGGTAAGGTTGTAACAAGAGGCGt      >  1:1367911/1‑90 (MQ=255)
                                                                                    gggggAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGt      <  2:377623/90‑1 (MQ=255)
                                                                                    gggggAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGt      <  1:164373/90‑1 (MQ=255)
                                                                                       ggAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTAtt   <  1:673962/90‑1 (MQ=255)
                                                                                         aaaCCAGATGCTTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTAttt  >  1:825604/1‑90 (MQ=255)
                                                                                         |                                                                                         
AGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTT  >  NZ_CP009273/4431874‑4432051

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