Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,392,049 T→C 9.2% intergenic (+74/‑41) aspS → / → PP_1214 aspartate‑‑tRNA ligase/transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,392,0490TC 9.2% 239.5 / 12.0 87intergenic (+74/‑41)aspS/PP_1214aspartate‑‑tRNA ligase/transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (38/41);  new base C (3/5);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACGGTCGCTGC  >  NC_002947/1391902‑1392196
                                                                                                                                                   |                                                                                                                                                   
ccccTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATgt                                                                                                                                                   <  2:465345/150‑2 (MQ=255)
  ccTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCCGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATgggg                                                                                                                                                 >  1:622674/1‑150 (MQ=255)
          ggTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCATGACggggcgggacgg                                                                                                                                         <  1:428745/150‑5 (MQ=255)
          ggTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGCCCCGCCCCATGACggggcgggacgg                                                                                                                                         <  2:984497/150‑5 (MQ=255)
           gTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCATGACggggcgggacggg                                                                                                                                        <  1:1200707/150‑6 (MQ=255)
            tCGACGCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGtt                                                                                                                                       >  1:3565/1‑150 (MQ=255)
                 gCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGa                                                                                                                                  >  2:376754/1‑150 (MQ=255)
                        cccTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAg                                                                                                                           <  2:147749/150‑1 (MQ=255)
                         ccTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGa                                                                                                                                  <  2:1814323/142‑1 (MQ=255)
                         ccTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGa                                                                                                                                  >  1:1814323/1‑142 (MQ=255)
                          cTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGa                                                                                                                         <  2:1793723/150‑1 (MQ=255)
                              gcgAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCg                                                                                                                     <  2:696971/150‑1 (MQ=255)
                                 aaCTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAg                                                                                                                  <  1:1130062/150‑1 (MQ=255)
                                   ggtcacaTCCGACTGAGAGAACAGACCCCGGTAGCGCGAGCAGCCCACCGTGGCAGCCACAGGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTc                                                                                                                <  2:511039/147‑1 (MQ=255)
                                   cTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTc                                                                                                                <  2:1367705/150‑1 (MQ=255)
                                    tGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCg                                                                                                               <  2:1155247/150‑1 (MQ=255)
                                      cacaTCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCATTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTTGACGCAGGCTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAACGATTCGGAGTCGtt                                                                                                             >  1:1295213/1‑150 (MQ=255)
                                      cacaTCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGtt                                                                                                             >  1:697608/1‑150 (MQ=255)
                                        caTCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTAt                                                                                                           <  1:376754/150‑1 (MQ=255)
                                          tCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATgg                                                                                                         >  2:132267/1‑150 (MQ=255)
                                            cGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGCTACGAAGGCGTTGTGGCt                                                                                                       >  2:1320927/1‑150 (MQ=255)
                                             gACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTg                                                                                                      <  2:3565/150‑1 (MQ=255)
                                             gACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTg                                                                                                      <  2:1070492/150‑1 (MQ=255)
                                              aCTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCTGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTgg                                                                                                     <  2:1006260/150‑1 (MQ=255)
                                              aCTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTgg                                                                                                     <  1:132267/150‑1 (MQ=255)
                                                tGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCCTGGAGCGTGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTc                                                                                                   <  2:273656/150‑1 (MQ=255)
                                                   gcgAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGCCCCGCCCCATGACggggcgggacg                                                                                                                                          >  1:1788915/1‑105 (MQ=255)
                                                   gcgAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGCCCCGCCCCATGACggggcgggacg                                                                                                                                          <  2:1788915/108‑4 (MQ=255)
                                                   ccgatctgACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTAGAAGGATTCGGAGTCGTTATGGCTGGTCAtt                                                                                                <  1:1022569/143‑1 (MQ=255)
                                                    cgAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTc                                                                                               <  1:1826287/150‑1 (MQ=255)
                                                         aGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGt                                                                                          <  2:697608/150‑1 (MQ=255)
                                                           aCCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTgg                                                                                        >  2:1662851/1‑150 (MQ=255)
                                                           aCCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGCGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTcagatcgg                                                                                        >  2:1022569/1‑143 (MQ=255)
                                                             cAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTc                                                                                               >  1:380547/1‑141 (MQ=255)
                                                             cAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTc                                                                                               <  2:380547/141‑1 (MQ=255)
                                                             cAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGc                                                                                      <  2:348760/150‑1 (MQ=255)
                                                             cAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGc                                                                                      >  1:1121984/1‑150 (MQ=255)
                                                             cAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGc                                                                                      <  1:170426/150‑1 (MQ=255)
                                                              aaGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCg                                                                                     <  2:1121984/150‑1 (MQ=255)
                                                                 gTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTTCTTTTTTTTCTCGTCATTGGGCTGGGCCGTTTGTGAGTCAAAGGATTTGGAGTCGGTATGGCTGTTCATTCCAAGTGGGCGGAc                                                                                  >  2:565078/1‑150 (MQ=255)
                                                                  tCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACa                                                                                 >  2:531718/1‑150 (MQ=255)
                                                                     aGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCGGGTTTTTTCCCGGCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCa                                                                              >  2:1041289/1‑150 (MQ=255)
                                                                     aGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCa                                                                              <  2:1651461/150‑1 (MQ=255)
                                                                                cccGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGa                                                                                                                                  >  2:474752/1‑87 (MQ=255)
                                                                                cccGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGa                                                                                                                                  <  1:474752/87‑1 (MQ=255)
                                                                                cccGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCcaa                                                                   >  1:491750/1‑150 (MQ=255)
                                                                                 ccGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTg                                                                                                      >  1:1204196/1‑114 (MQ=255)
                                                                                 ccGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTg                                                                                                      <  2:1204196/114‑1 (MQ=255)
                                                                                     gcgcATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAgcg                                                              >  1:1056163/1‑150 (MQ=255)
                                                                                        cATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCa                                                           <  2:491750/150‑1 (MQ=255)
                                                                                          tCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAgg                                                         >  2:1820745/1‑150 (MQ=255)
                                                                                             acaTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTAt                                                                                                           >  2:1564211/1‑97 (MQ=255)
                                                                                             acaTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTAt                                                                                                           <  1:1564211/97‑1 (MQ=255)
                                                                                              tgtGGATGCGCCCTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGc                                                     >  2:1014024/3‑150 (MQ=255)
                                                                                               aTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCg                                                                                                                     <  2:832597/85‑1 (MQ=255)
                                                                                               aTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCg                                                                                                                     >  1:832597/1‑85 (MQ=255)
                                                                                                    tGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGc                                                     <  2:1275893/144‑1 (MQ=255)
                                                                                                    tGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGc                                                     >  1:1275893/1‑144 (MQ=255)
                                                                                                     gctcgtgtGCGTTTCAACGCTGGTTTTTTCCCGTCCCGCCCCGTTATGGGGCGGGGCCGTTTGTGATTAAAAGGATTTGGAGTCGTTATGGCTGGTCAGTCCATGTGGGCGAACATCAAGCCCCGCAAAGAGCGCCCCGCTGCCAagaga                                              <  2:1295213/143‑1 (MQ=255)
                                                                                                     gcgcCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagaga                                              <  2:622674/150‑1 (MQ=255)
                                                                                                     gcgcCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagaga                                              >  1:1534007/1‑150 (MQ=255)
                                                                                                        ccTTTGTGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGc                                           >  1:141884/1‑150 (MQ=255)
                                                                                                           ttGCGTTTCGACGCAGGTTTTTTTTTTTTTTTTCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTTGGCGTACATCAAGCACCGCAAAGATCGCCAGGATGCCAAGAGAGGCAAg                                        >  2:1196938/1‑150 (MQ=255)
                                                                                                           ttGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAg                                        <  2:1534007/150‑1 (MQ=255)
                                                                                                            tGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCACGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAgg                                       <  1:651545/150‑1 (MQ=255)
                                                                                                             gCGTTTCGACGCAGGTTTTTTCCCGTTCCGCCACGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTGGGAGTGGTTATGGCTGGTTATTCGAAGTGGGCGAACATCAAGCACCGGAAAGAGCGTTAGGATGCCAAGAGAGGTAAGGt                                      >  2:651545/1‑150 (MQ=255)
                                                                                                              cGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTc                                     <  1:1662851/150‑1 (MQ=255)
                                                                                                               gTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCt                                    >  1:1701529/1‑150 (MQ=255)
                                                                                                               gTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGGCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCt                                    >  2:1206139/1‑150 (MQ=255)
                                                                                                                tttCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCa                                                                          <  1:160256/111‑1 (MQ=255)
                                                                                                                tttCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCa                                                                          >  2:160256/1‑111 (MQ=255)
                                                                                                                         aGGTTTTTTCCCGTCCCACCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTg                          >  1:1112880/1‑150 (MQ=255)
                                                                                                                                  cccgtcccgCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGTCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGa                 >  2:1496373/1‑150 (MQ=255)
                                                                                                                                    ggccccgCCCCATGACGGGGCGGGACCGTTTGTGATTCAAAGGATTCGGAGTCGGTATGGCTGGTCATTCCAAGTGGGTGTAGTGGATGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTGAACAAGTGGATCCGTGAAc               >  1:1088763/4‑150 (MQ=255)
                                                                                                                                     gtcccgCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGAGAAAAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACt              >  1:1070630/1‑150 (MQ=255)
                                                                                                                                      tcccgCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTg             >  2:60516/1‑150 (MQ=255)
                                                                                                                                      ccccgCCCCATGACGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCa                                                           >  1:642919/2‑104 (MQ=255)
                                                                                                                                      ccccgCCCCATGACGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCa                                                           <  2:642919/103‑1 (MQ=255)
                                                                                                                                       cccgCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGa            <  1:1041289/150‑1 (MQ=255)
                                                                                                                                       cccgCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGGGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGa            >  1:924384/1‑150 (MQ=255)
                                                                                                                                         cgCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACg          >  2:2867/1‑150 (MQ=255)
                                                                                                                                         cgCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACg          <  1:2867/150‑1 (MQ=255)
                                                                                                                                          gCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACgg         <  2:1056163/150‑1 (MQ=255)
                                                                                                                                               gTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACGGTCGCt    <  1:60516/150‑1 (MQ=255)
                                                                                                                                                tCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACGGTCGCTg   >  1:1939702/1‑150 (MQ=255)
                                                                                                                                                 cATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACGGTCGCTGc  <  2:1701529/150‑1 (MQ=255)
                                                                                                                                                 cATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACGGTCGCTGc  <  2:141884/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGTGAACTGACGGTCGCTGC  >  NC_002947/1391902‑1392196

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: