Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,168,9660AC38.2% 104.4 / 17.8 69I85I (ATA→ATC) yegSphosphatidylglycerol kinase, metal‑dependent
Reads supporting (aligned to +/- strand):  ref base A (31/11);  new base C (0/26);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 2.02e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTGCGGGTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACA  >  NC_000913/2168818‑2169106
                                                                                                                                                    |                                                                                                                                             
aTGTGCGGGTCACCTGGGAGAAAGTCGATGCCGCACGATATGTAGAGGAGACCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGCTGGCACCATTAATGAAGTTTCTACGGCGTTGATTCCGTGTGAGGGGGATCACATAcc                                                                                                                                             >  2:831214/1‑151 (MQ=255)
 tgtgCGGGTCACCTGGGAGAAAAGCGATGCCGCACGACATGTAGAGGAGGCCCGGCAGTTTGGCGTCGCAACGGTGATTGACGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCCGAGTGAGGGGGATGACATAccc                                                                                                                                            >  1:1043413/1‑151 (MQ=255)
     cGGGTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCgcgc                                                                                                                                        >  1:1016078/1‑151 (MQ=255)
     cGGGTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCgcgc                                                                                                                                        >  2:671083/1‑151 (MQ=255)
       ggTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTg                                                                                                                                      <  2:172962/151‑1 (MQ=255)
       ggTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTg                                                                                                                                      <  1:320028/151‑1 (MQ=255)
       ggTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTg                                                                                                                                      <  2:316465/151‑1 (MQ=255)
                 gagaAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGc                                                                                                                           >  1:629059/1‑151 (MQ=255)
                          gATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGaa                                                                                                                  <  1:1022686/151‑1 (MQ=255)
                               cGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGcc                                                                                                             <  2:147518/151‑1 (MQ=255)
                                     atatGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGAt                                                                                                       >  1:1011627/1‑151 (MQ=255)
                                      tatGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGAtt                                                                                                      >  2:283724/1‑151 (MQ=255)
                                       atGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGGGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGAGAGGGGGATGACATACGCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGAt                                                                                                       >  1:1163359/1‑149 (MQ=255)
                                        tGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGAtttt                                                                                                    >  1:1167581/1‑151 (MQ=255)
                                        tGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGAtttt                                                                                                    >  1:1193289/1‑151 (MQ=255)
                                        tGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGAtttt                                                                                                    >  2:934028/1‑151 (MQ=255)
                                          tAGAGGAGGCCCAGAAGTTTGGCGTCGAAACGGTGATTGCGGGTGGTGGGGAGGGCCCCATTAATGAAGTTTCTCCGGGGTTGATTCAGTGTGAGGGGGATGACATCCCCGCGCGGGGAA‑TTTTGCCATTGGTAACCGCCAATGAtttt                                                                                                    <  1:983208/149‑1 (MQ=255)
                                            gaggagGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGc                                                                                                  >  1:442972/1‑149 (MQ=255)
                                              ggaggCCGGGAATTTGGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCTTTAATGAAGTTTCTACGGCTTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGccacc                                                                                              <  2:841186/151‑1 (MQ=255)
                                               ggggCCGGGAGGTTTGGCTTCGCACCGGTGTTTGCCGGTGGTGGCGGTGGCTCCATTAATGAAGTTTCTGCGGCGTTATTTCAGTGTGAGGGGGATGACATCCCCGCGCGGGGAA‑TTTTGCCATTAGAAACCGCCAATGATTTTGccacca                                                                                             <  2:935720/149‑1 (MQ=255)
                                                 ggCCCGGAAGTTTGGCGTCGCAACGGTGATTGCGGGTGGTGGCGATGGCCCCATTAATGAAGTTTCTCCGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAg                                                                                            <  2:580475/150‑1 (MQ=255)
                                                   cccGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAgtgt                                                                                         >  1:657869/1‑151 (MQ=255)
                                                        aaGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATCCCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGa                                                                                    <  2:251579/151‑1 (MQ=255)
                                                        aaGTTGGGCGTCGCAGCGGTGATTGCGGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAT‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGa                                                                                    <  1:255926/151‑1 (MQ=255)
                                                          gTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCACCAATGATTTTGCCACCAGTGTAGGGAtt                                                                                  >  1:748429/1‑151 (MQ=255)
                                                              ggCTTCGCAACGGTGATTCCCGGTGGTGGCGATGGCACCTTTAATGAATTTTCTACGGCTTTGATTCAGTGTGGGGGGGTTGACTTCCCCGCGCTGGAAA‑TTTGGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTg                                                                              <  2:769590/151‑1 (MQ=255)
                                                                  tCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAGGTTTCTAGGGGGTTGATTCAGTGTGGGGGGGATGACATCCCCGCGCTGGAAA‑TTTTGCCGTTAGGACCCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGa                                                                             <  2:1112741/148‑1 (MQ=255)
                                                                    ccAACGGTGTTTCCCGGTGGGGGCGTTGGAACCATAAATGAAGTTTCTACGGCTTGAATACAGGGTGGGGGGGAGGCCATCCCCGCGTGGGGAA‑TTTCGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCa                                                                         <  2:471371/149‑1 (MQ=255)
                                                                      aaCGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTg                                                                      >  1:134812/1‑151 (MQ=255)
                                                                            gATGGCGGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGAAA‑TTTTGCAATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATaa                                                                 <  2:1103122/150‑1 (MQ=255)
                                                                              ttGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGgc                                                              >  1:693101/1‑151 (MQ=255)
                                                                                gcgggtggtGGCGATGGCCCCATTAATGAATTTTCTACGGCTTTGATTCAGTGGGGGGGGGATGACATCCCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGgc                                                              <  1:1115445/146‑1 (MQ=255)
                                                                                 ccGGTGGTGGCGATGGCACCATAAATGAATTTTCTACGGCGTGGATTCAGTGTGAGGGGGAGGACATCCCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGgcg                                                             <  2:1184990/149‑1 (MQ=255)
                                                                                 ccGGTGGGGGCGATGGCACCATTAAGGAAGTTTCTCCGGTGTTGATTCAGTGGGAGGGGGATGACATCCCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGgcgc                                                            <  1:1204936/150‑1 (MQ=255)
                                                                                   ggtggtGGGGATGGCACCATTAATGAAGTTTCTCCGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGAAA‑TTTTGCCTTTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGa                                                         <  1:989509/151‑1 (MQ=255)
                                                                                      ggtggGGAGGGCCCCATTAATGAAGTTTCTCCGGCGTTGATTCAGTGTGAGGGGGAAGACATCCCCGCGGTGGGAA‑TTTTGCCATTGGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAAc                                                      <  2:124090/151‑1 (MQ=255)
                                                                                      ggtggCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAAc                                                      >  2:821058/1‑151 (MQ=255)
                                                                                         ggCGTTGGCCCAATTAATGAAGTTTCTCCGGCTTTTTTTCAGTGTGGGGGGGATGACATCCCCGGCCTGGGAT‑TTTTGCCATTAGGAGCCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTg                                                    <  1:996229/150‑1 (MQ=255)
                                                                                         ggCGATGGCACCATTAATGAAGTTTCTGCGGCGTTGATTCAGTGTGAGGGGGATGACATCCCCGCGCGGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACt                                                     <  1:630330/149‑1 (MQ=255)
                                                                                          gCGATGGCACCATTAATGAAGTTTCTACGGCGTTGTTTCGGTGTGAGGGGGATTACATCCCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGc                                                  <  2:134816/151‑1 (MQ=255)
                                                                                                gCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGc                                            >  1:760877/1‑151 (MQ=255)
                                                                                                 cACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGcc                                           >  1:411243/1‑151 (MQ=255)
                                                                                                  aCCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCg                                          >  1:574327/1‑151 (MQ=255)
                                                                                                    cATTAATGAATTTTCTCGGGCGTTGATTCAGTGTGAGGGGGATGACATCCCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGcc                                           <  1:355083/148‑1 (MQ=255)
                                                                                                      ttAATGAAGTTTCTACGGCGTTGATTCGGGGTGAGGGGGATGACATCCCCG‑GGTGGGAATTTTTGCCATTAGGAGCCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTg                                       <  2:716161/150‑1 (MQ=255)
                                                                                                      ttAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACGGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTg                                       >  1:1025366/1‑150 (MQ=255)
                                                                                                      ttAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTg                                       >  2:870430/1‑150 (MQ=255)
                                                                                                       tAATGAAGTTTCTACGGCGTTGTTTCAGTGGGGGGGGGACGCCATCCCCGGGCTGGGAA‑TTTTCCCAGTAGAAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGa                                      <  1:765828/150‑1 (MQ=255)
                                                                                                       tAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGAc                                     >  2:17966/1‑151 (MQ=255)
                                                                                                                 tCTGCGGCTTTGTTTCAGGGTGAGGGGGATGAAATCCCCGCGCTGGGAT‑TTTGGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTg                              <  2:846285/148‑1 (MQ=255)
                                                                                                                      ggggTTGATTCAGTGTGAGGGGGATGACATCCCCGCGCTGGGAT‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCgatagat                      <  2:1213929/148‑1 (MQ=255)
                                                                                                                      ggCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCgat                          >  2:1185782/1‑147 (MQ=255)
                                                                                                                      ggCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCgat                          >  2:780399/1‑147 (MQ=255)
                                                                                                                         gTTGATTCAGTGTGCGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATg                   >  2:95948/1‑151 (MQ=255)
                                                                                                                          ttGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATgg                  >  2:306358/1‑151 (MQ=255)
                                                                                                                          ttGATACAGTGGGGGGGGGATGACATCCCCGGGCTGGTAA‑TTTTGCCTTTAGGACCCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATg                   <  2:828706/150‑1 (MQ=255)
                                                                                                                              ttCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGgcgc               >  1:680898/1‑150 (MQ=255)
                                                                                                                                 aGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGt           >  1:1224559/1‑151 (MQ=255)
                                                                                                                                    gtgAGGGGGATGACATCCCCGCGCTGGGAT‑TTTTCCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTcaa        <  1:671058/151‑1 (MQ=255)
                                                                                                                                    gtgAGGGGGATGACATCCCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTca         <  2:218151/150‑1 (MQ=255)
                                                                                                                                    gggggggggATGACATCCCCGCGCGGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTcaa        <  1:476838/147‑1 (MQ=255)
                                                                                                                                     tgAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTcaac       >  1:1222786/1‑151 (MQ=255)
                                                                                                                                      gAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTcaac       >  2:540140/1‑150 (MQ=255)
                                                                                                                                        gggggATGACATCCCCGCGCGGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCaacaaa    <  1:169062/151‑1 (MQ=255)
                                                                                                                                        gggggATGACACGCCCGCGCTGGGAA‑TTTTGCCATTGGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTcaacaa     <  2:27565/150‑1 (MQ=255)
                                                                                                                                        gggggATAACATCCCCGCGCGGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTcaa        <  1:669580/147‑1 (MQ=255)
                                                                                                                                        gggggATAACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTcaacaa     <  1:932197/150‑1 (MQ=255)
                                                                                                                                          gggATGACATCCCCGCGTTGGGAT‑TTTTCCCATTGGTAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCaacaaaca  <  2:683807/151‑1 (MQ=255)
                                                                                                                                          gggATGACATCCCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCaacaaaca  <  1:165960/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ATGTGCGGGTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAA‑TTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACA  >  NC_000913/2168818‑2169106

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