Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,355,480:1 +T 8.6% coding (5/264 nt) PP_1180 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,355,4801.T 8.6% 265.6 / 19.8 70coding (5/264 nt)PP_1180hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (34/30);  new base T (3/3);  total (37/33)
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.65e-01

GGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAACACCCAT  >  NC_002947/1355332‑1355629
                                                                                                                                                       |                                                                                                                                                     
ggcCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAatgcat                                                                                                                                                      <  1:1280287/150‑1 (MQ=255)
        gCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATcc                                                                                                                                              >  1:383688/1‑150 (MQ=255)
          aCCTGGGCCTGGGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTACCCAGGTCCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTTCAGTTGCTCGATAACGTCGCGA‑‑CGTTGCCCGATGCATGCGCagctcc                                                                                                                                            >  2:308453/1‑145 (MQ=255)
             tGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGCCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATccgcccg                                                                                                                                         >  1:818883/1‑150 (MQ=255)
              gggCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCATGCATTGCGCATccccgcga                                                                                                                                        <  1:40647/150‑7 (MQ=255)
                         ttCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTAcc                                                                                                                             >  2:1244521/1‑150 (MQ=255)
                            tcGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGAc                                                                                                                          >  1:629768/1‑150 (MQ=255)
                             cgcgAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACt                                                                                                                         <  1:1586432/150‑1 (MQ=255)
                              gcgAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCAAAGTTGCGATATTGTTCAGGTGCCTTATGGGCGCTCACCGTTCCGTGTACTTTCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGACGTACCGACtt                                                                                                                        >  2:663193/1‑150 (MQ=255)
                              gcgAGTCGTGGCGGTTTGCATGGGTGGGTGTCTCCCCAGTAACTATACTGTTCAGGCGCCTTATGGCAGCCCACCGTTACGTTGACTTGCTCGATAACGTAGCAG‑‑CGTTGTGGAGTGCAAGAGTATCCACACGCGCGGGGGTACCgtttt                                                                                                                        >  2:691193/1‑146 (MQ=255)
                                     gtgtCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACa                                                                                                                 <  1:1158338/150‑1 (MQ=255)
                                            tttGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAg                                                                                                          <  1:1244521/150‑1 (MQ=255)
                                              cgCATCGGTGGGTGTATCACCGGGTTATATATTGGACAGGCGCTTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGa                                                                                                        <  1:691193/149‑1 (MQ=255)
                                               gCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGaa                                                                                                       <  2:790271/150‑1 (MQ=255)
                                                cATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAAt                                                                                                      >  1:943216/1‑150 (MQ=255)
                                                  tgggtgggTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATaa                                                                                                    <  2:279254/150‑1 (MQ=255)
                                                   gggtgggtGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGTCTTCGGGATAGACTCAGGAATAAg                                                                                                   <  2:1335189/150‑1 (MQ=255)
                                                   gggtgggtGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAg                                                                                                   <  1:308453/150‑1 (MQ=255)
                                                      tgggtgTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAgggg                                                                                                <  2:629768/150‑1 (MQ=255)
                                                              tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAggg                                                                                                 >  1:1499242/1‑141 (MQ=255)
                                                              tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAggg                                                                                                 <  2:1499242/141‑1 (MQ=255)
                                                              tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCt                                                                                        >  2:1165817/1‑150 (MQ=255)
                                                              tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCt                                                                                        >  2:1076603/1‑150 (MQ=255)
                                                               ccACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCTCGATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCtt                                                                                       <  2:1621922/150‑1 (MQ=255)
                                                                 acaGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGc                                                                                     <  1:1076603/150‑1 (MQ=255)
                                                                 acaGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCGCGCGGGCGGATGCGCATGCATTGCGCa                                                                                                                                                 >  2:204445/1‑92 (MQ=25)
                                                                 acaGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCGCGCGGGCGGATGCGCATGCATTGCGCa                                                                                                                                                 <  1:204445/92‑1 (MQ=25)
                                                                               gTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGa                                                                       >  1:1461182/1‑150 (MQ=255)
                                                                                ttCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGaa                                                                      >  1:739586/1‑150 (MQ=255)
                                                                                   aGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGT‑‑GGATGCGCAATGCATGCGCATCAGCCCGCGCGGGCGTACCGACTTCGGGACAGAATGAGGAATAAGGGGCAATGGCGTGCCGGTGCCTGGCCGGAgac                                                                   >  1:491640/1‑148 (MQ=255)
                                                                                   aGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTc                                                                   >  2:1537549/1‑150 (MQ=255)
                                                                                     gcgcCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCaa                                                                 >  2:1301968/1‑150 (MQ=255)
                                                                                     gcgcCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCaa                                                                 >  2:1262578/1‑150 (MQ=255)
                                                                                      cgcCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCaaa                                                                <  1:1262578/150‑1 (MQ=255)
                                                                                        ccTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                                   <  1:514859/97‑1 (MQ=255)
                                                                                        ccTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                                   >  2:514859/1‑97 (MQ=255)
                                                                                            aTGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCg                                                                                                                      >  1:1166768/1‑90 (MQ=255)
                                                                                            aTGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCg                                                                                                                      <  2:1166768/90‑1 (MQ=255)
                                                                                            aTGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGa                                                                                                        <  1:581499/104‑1 (MQ=255)
                                                                                            aTGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGa                                                                                                        >  2:581499/1‑104 (MQ=255)
                                                                                              gggCGCCTACCGTTCCGTCCACTTGCTCGAGAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGCCTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTgg                                                        <  2:701252/150‑1 (MQ=255)
                                                                                              gggCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTgg                                                        >  2:616637/1‑150 (MQ=255)
                                                                                              gggCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCCGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTgg                                                        <  1:425689/150‑1 (MQ=255)
                                                                                               ggCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTggg                                                       >  1:1166037/1‑150 (MQ=255)
                                                                                                   cccACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGaa                                                                      >  1:1322318/1‑131 (MQ=255)
                                                                                                   cccACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGaa                                                                      <  2:1322318/131‑1 (MQ=255)
                                                                                                      aCCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATg                                                                                           <  1:1068269/107‑1 (MQ=255)
                                                                                                      aCCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATg                                                                                           >  2:1068269/1‑107 (MQ=255)
                                                                                                             cGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAAc                                         >  1:943930/1‑150 (MQ=255)
                                                                                                              gTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACt                                        >  1:791238/1‑150 (MQ=255)
                                                                                                               tCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTg                                       <  1:616637/150‑1 (MQ=255)
                                                                                                                  aCTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTg                                    >  2:1252736/1‑150 (MQ=255)
                                                                                                                  aCTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTg                                    <  2:791238/150‑1 (MQ=255)
                                                                                                                       cTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACGTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGAt                               <  1:1252736/150‑1 (MQ=255)
                                                                                                                        tgcaTAACCTCGCGG‑‑GGATGCGCAATGCATGCGGATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATg                              <  2:312978/147‑1 (MQ=255)
                                                                                                                          gATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGc                                                                                     <  2:117437/93‑1 (MQ=255)
                                                                                                                          gATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGc                                                                                     >  1:117437/1‑93 (MQ=255)
                                                                                                                           aTAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGt                           <  2:943930/150‑1 (MQ=255)
                                                                                                                               cctcaagg‑‑ggATGCGCAATGCATGCGCATCCGGCAGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCCCCCTGGCCGAAGTCAAGTGTTCTGGGCCGAGCCAGAAAACTGTTGCTGCTGGGTGCCa                       >  2:380995/7‑150 (MQ=255)
                                                                                                                                   gcgggcggATGCGCATGCATTGCGCATC‑‑CCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGt                                      >  1:771078/7‑131 (MQ=38)
                                                                                                                                   gcgggcggATGCGCATGCATTGCGCATC‑‑CCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGt                                      <  2:771078/125‑1 (MQ=38)
                                                                                                                                   gcgG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGa                                                                                                               >  2:1069831/1‑58 (MQ=255)
                                                                                                                                   gcgG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGa                                                                                                               <  1:1069831/58‑1 (MQ=255)
                                                                                                                                         ggATGCGCATGCATTGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctggc               >  2:586796/1‑150 (MQ=255)
                                                                                                                                           aTGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctggcct             >  1:1066104/1‑150 (MQ=255)
                                                                                                                                           aTGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctggcct             <  1:1165817/150‑1 (MQ=255)
                                                                                                                                              cgcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGAt                               <  2:159829/129‑1 (MQ=255)
                                                                                                                                              cgcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGAt                               >  1:159829/1‑129 (MQ=255)
                                                                                                                                               gcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCaa         >  2:224218/1‑150 (MQ=255)
                                                                                                                                                 aaTGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGc                                                                          <  2:1494729/83‑1 (MQ=255)
                                                                                                                                                 aaTGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGc                                                                          >  1:1494729/1‑83 (MQ=255)
                                                                                                                                                  atgcattgCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGa                                                                                                               >  1:1503674/1‑46 (MQ=255)
                                                                                                                                                  atgcattgCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGa                                                                                                               <  2:1503674/46‑1 (MQ=255)
                                                                                                                                                      atgcGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGc                                   >  1:1372697/1‑117 (MQ=255)
                                                                                                                                                      atgcGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGc                                   <  2:1372697/117‑1 (MQ=255)
                                                                                                                                                      atgcGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAACACCCAt  <  2:1461182/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAACACCCAT  >  NC_002947/1355332‑1355629

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