Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 564,380:1 +A 9.6% intergenic (+29/‑144) rplQ → / → katA 50S ribosomal protein L17/catalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947564,3801.A 9.6% 338.6 / 31.2 94intergenic (+29/‑144)rplQ/katA50S ribosomal protein L17/catalase
Reads supporting (aligned to +/- strand):  ref base . (42/43);  new base A (4/5);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

CCGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAG  >  NC_002947/564231‑564528
                                                                                                                                                       |                                                                                                                                                    
ccGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGc                                                                                                                                                       <  2:782739/150‑1 (MQ=255)
        gcggcgACCTGTCTATCCTGAACTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGGTCGTCCGGTCGGCGGTGCTGTAGAAGCTTCTGAGTAAGGCGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGcc                                                                                                                                              <  2:763988/150‑1 (MQ=255)
          ggcTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCg                                                                                                                                            >  1:866568/1‑150 (MQ=255)
             tACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCACTAGGCCCGGtt                                                                                                                                         <  1:28953/150‑1 (MQ=255)
             tACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGtt                                                                                                                                         >  1:1086115/1‑150 (MQ=255)
                cTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGttttt                                                                                                                                      <  2:177986/150‑1 (MQ=255)
                cTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAGGCCCGGttttc                                                                                                                                      >  1:1934733/1‑149 (MQ=255)
                  gcgcATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTg                                                                                                                                    <  2:294873/150‑1 (MQ=255)
                     cATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGttt                                                                                                                                 <  2:2326512/150‑1 (MQ=255)
                      aTCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTta                                                                                                                                >  2:2260962/1‑150 (MQ=255)
                       tCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAGGCCCGGTTTTCTGtagca                                                                                                                               <  2:1934733/150‑5 (MQ=38)
                        ccTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGtt                                                                                                                                  <  1:938961/146‑1 (MQ=255)
                        ccTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGGGCCCGGTTTTTTGtt                                                                                                                                  >  2:938961/1‑146 (MQ=255)
                          tGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTta                                                                                                                                <  1:1381198/146‑1 (MQ=255)
                          tGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTta                                                                                                                                >  2:1381198/1‑146 (MQ=255)
                          tGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACa                                                                                                                            >  2:350763/1‑150 (MQ=255)
                           gAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAACCCTGGCCGAGGGCGCGGTTGTTTGGGTATACAt                                                                                                                           <  2:2233126/150‑1 (MQ=255)
                              gTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTa                                                                                                                        >  1:1707608/1‑150 (MQ=255)
                                gCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTAtt                                                                                                                      <  1:2260962/150‑1 (MQ=255)
                                  ggTTCCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGCAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGt                                                                                                                    >  2:1172243/1‑150 (MQ=255)
                                     ttCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTAATCGTCAGAAAACTGATCCTTTTTCTCGTGTTTTTGTTTATACATCTATTGTTTa                                                                                                                 <  1:299351/150‑1 (MQ=255)
                                      tCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAt                                                                                                                <  2:540128/150‑1 (MQ=255)
                                         gcgcTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTc                                                                                                             >  1:515017/1‑150 (MQ=255)
                                               gggACAACGCGCCTGTGGGGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTGCAGCACACCGTGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGa                                                                                                       <  2:332234/148‑1 (MQ=255)
                                                 gACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATg                                                                                                     >  2:90084/1‑150 (MQ=255)
                                                    aaCNCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTaa                                                                                                  >  2:1498428/1‑150 (MQ=255)
                                                    aaCGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTaa                                                                                                  <  1:763439/150‑1 (MQ=255)
                                                    aaCGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTaa                                                                                                  >  1:1111118/1‑150 (MQ=255)
                                                       gcgcCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTg                                                                                                                     <  2:937265/128‑1 (MQ=255)
                                                       gcgcCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTg                                                                                                                     >  1:937265/1‑128 (MQ=255)
                                                       gcgcCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTc                                                                                               >  2:1153673/1‑150 (MQ=255)
                                                            tATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGtt                                                                                                                                  <  1:1026156/110‑1 (MQ=255)
                                                            tATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGtt                                                                                                                                  >  2:1026156/1‑110 (MQ=255)
                                                            tATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAtt                                                                                                         <  1:578775/135‑1 (MQ=255)
                                                            tATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAtt                                                                                                         >  2:578775/1‑135 (MQ=255)
                                                             aTGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTaa                                                                                         <  1:90084/150‑1 (MQ=255)
                                                                 cGTGTGTCGAGCTGGTTGACCGTCCGGTCGGAGGCGCTGTGCCTCGAGTGTTTCGGCTCTTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCCGTTTTTTGTTTATACCTCTATTCTTTATTTCTATTGATGTAAGTCATGTAATCaa                                                                                     <  2:1243471/150‑1 (MQ=255)
                                                                       tCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTc                                                                               <  2:339438/150‑1 (MQ=255)
                                                                          aGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCt                                                                            <  1:1459774/150‑1 (MQ=255)
                                                                              ggTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTaa                                                                        <  1:2071225/150‑1 (MQ=255)
                                                                               gTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAAt                                                                       >  2:418002/1‑150 (MQ=255)
                                                                               gTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAAt                                                                       >  2:419039/1‑150 (MQ=255)
                                                                                ttGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGt                                                                                                                                                 >  1:71526/1‑75 (MQ=255)
                                                                                ttGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGt                                                                                                                                                 <  2:71526/75‑1 (MQ=255)
                                                                                  gATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCtt                                                                    >  2:1988765/1‑150 (MQ=255)
                                                                                   gtCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTg                                                                   >  2:1150329/2‑150 (MQ=255)
                                                                                    tCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGc                                                                  >  2:2001146/1‑150 (MQ=255)
                                                                                     cGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCt                                                                 <  1:1684748/150‑1 (MQ=255)
                                                                                          ggTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGCAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGCATTTGTCCCTGTAATCTTGCTCGTCa                                                            >  2:778509/1‑150 (MQ=255)
                                                                                             cggcggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATa                                                         <  1:418002/150‑1 (MQ=255)
                                                                                             cggcggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATa                                                         <  2:515017/150‑1 (MQ=255)
                                                                                             cggcggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATa                                                         <  1:419039/150‑1 (MQ=255)
                                                                                                cggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActc                                                      <  1:2370947/150‑1 (MQ=255)
                                                                                                 ggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActct                                                     >  1:1183830/1‑150 (MQ=255)
                                                                                                    gCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctc                                                  <  1:1498428/150‑1 (MQ=255)
                                                                                                    gCTGTAGAAGCTGCTGAGTAAGACGTTCTGCTGGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctc                                                  <  1:1451624/150‑1 (MQ=255)
                                                                                                       gTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTAc                                               >  1:1222383/1‑150 (MQ=255)
                                                                                                        tAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACa                                              >  1:1713576/1‑150 (MQ=255)
                                                                                                           aaGCTGCTGAGTAAGACGTTCTGTCTGCCACATAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGcc                                           <  1:2087374/150‑1 (MQ=255)
                                                                                                               tgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGAtt                                       >  1:2319442/1‑150 (MQ=255)
                                                                                                               tgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGAtt                                       >  1:2318914/1‑150 (MQ=255)
                                                                                                                gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTa                                                                                                                                                   <  2:945272/41‑1 (MQ=255)
                                                                                                                gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTa                                                                                                                                                   >  1:945272/1‑41 (MQ=255)
                                                                                                                 ctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAg                                     >  1:1132791/1‑150 (MQ=255)
                                                                                                                 ctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGCGCCCGGCTTTCTGTTTCTACATCTATTGTTTATTTCTATTGATGTAAGCCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCACTACCCTCTCTCCAGCCGACTac                                     >  2:2165669/1‑149 (MQ=255)
                                                                                                                  tgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGt                                                                                                                                   >  2:875621/1‑55 (MQ=255)
                                                                                                                  tgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGt                                                                                                                                   <  1:875621/55‑1 (MQ=255)
                                                                                                                  tgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGc                                    >  1:909817/1‑150 (MQ=255)
                                                                                                                    aGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCg                                  <  2:1183830/150‑1 (MQ=255)
                                                                                                                      tAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGc                                >  2:588206/1‑150 (MQ=255)
                                                                                                                       aaGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGcc                               <  1:2373272/150‑1 (MQ=255)
                                                                                                                         gACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCt                             >  2:839887/1‑150 (MQ=255)
                                                                                                                           aCGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCtt                            <  2:1654264/150‑1 (MQ=255)
                                                                                                                             gTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGt                                                                                           <  2:99467/85‑1 (MQ=255)
                                                                                                                             gTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGACATGt                                                                                           >  1:99467/1‑85 (MQ=255)
                                                                                                                              ttCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGaa                         <  2:1111118/150‑1 (MQ=255)
                                                                                                                                ctgtctgCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGa                       <  1:1153673/150‑1 (MQ=255)
                                                                                                                                        tACAGAAAACCGGGCCTAGTGCCCTGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAg                                             >  2:1620766/1‑120 (MQ=255)
                                                                                                                                        tACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAg                                             <  1:1620766/120‑1 (MQ=255)
                                                                                                                                          cAGAAAACCGGGCCTAGTGCCCGGTTTTTTCTTTATACATCTATTGTTTATTTCTAATGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGga             <  2:866568/150‑1 (MQ=255)
                                                                                                                                           aaaaaaCCGGGCACTAGGCCCGGTTTTCTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGACCTTGAAGACCGAAATGgag            >  1:1304715/3‑150 (MQ=255)
                                                                                                                                            gAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGgaga           >  2:450275/1‑150 (MQ=255)
                                                                                                                                              aaaCCGGGCACTAGGCCCGGTTTTCTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGc         <  2:1304715/150‑1 (MQ=255)
                                                                                                                                               aaCCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTg                                                                   <  2:596447/91‑1 (MQ=255)
                                                                                                                                               aaCCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTg                                                                   >  1:596447/1‑91 (MQ=255)
                                                                                                                                               aaCCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACa                                              >  2:1389226/1‑112 (MQ=255)
                                                                                                                                               aaCCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACa                                              <  1:1389226/112‑1 (MQ=255)
                                                                                                                                               aaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAg                                     >  1:1267719/1‑121 (MQ=255)
                                                                                                                                               aaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAg                                     <  2:1267719/121‑1 (MQ=255)
                                                                                                                                               aaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGcc        <  1:1891618/150‑1 (MQ=255)
                                                                                                                                               aaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGcc        >  2:1534596/1‑150 (MQ=255)
                                                                                                                                                  cGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCAt     >  1:2294941/1‑150 (MQ=255)
                                                                                                                                                   gggCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATg    <  1:1899308/150‑1 (MQ=255)
                                                                                                                                                   gggCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATg    <  1:1150329/150‑1 (MQ=255)
                                                                                                                                                     gCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAg  >  2:1418526/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CCGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAG  >  NC_002947/564231‑564528

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