Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,355,475 Δ1 bp 8.7% intergenic (‑93/‑1) nrdA ← / → PP_1180 ribonucleoside‑diphosphate reductase subunit alpha/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,355,4750A. 8.7% 204.6 / 14.8 69intergenic (‑93/‑1)nrdA/PP_1180ribonucleoside‑diphosphate reductase subunit alpha/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base . (3/3);  total (36/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.53e-01

TCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAA  >  NC_002947/1355327‑1355622
                                                                                                                                                      |                                                                                                                                                    
tCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCaa                                                                                                                                                     >  2:1158338/1‑150 (MQ=255)
     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                                                                                          <  2:1499242/141‑1 (MQ=255)
                                                                   tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAggg                                                                                          >  1:1499242/1‑141 (MQ=255)
                                                                   tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCt                                                                                 >  2:1076603/1‑150 (MQ=255)
                                                                   tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCt                                                                                 >  2:1165817/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                                                                                                               <  2:1166768/90‑1 (MQ=255)
                                                                                                 aTGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCg                                                                                                               >  1:1166768/1‑90 (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                                                               <  2:1322318/131‑1 (MQ=255)
                                                                                                        cccACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGaa                                                               >  1:1322318/1‑131 (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                                                                                                        <  1:1069831/58‑1 (MQ=255)
                                                                                                                                        gcgG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGa                                                                                                        >  2:1069831/1‑58 (MQ=255)
                                                                                                                                              ggATGCGCATGCATTGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctggc        >  2:586796/1‑150 (MQ=255)
                                                                                                                                                aTGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctggcct      <  1:1165817/150‑1 (MQ=255)
                                                                                                                                                aTGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctggcct      >  1:1066104/1‑150 (MQ=255)
                                                                                                                                                   cgcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGAt                        >  1:159829/1‑129 (MQ=255)
                                                                                                                                                   cgcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGAt                        <  2:159829/129‑1 (MQ=255)
                                                                                                                                                    gcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCaa  >  2:224218/1‑150 (MQ=255)
                                                                                                                                                      aaTGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGc                                                                   >  1:1494729/1‑83 (MQ=255)
                                                                                                                                                      aaTGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGc                                                                   <  2:1494729/83‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAA  >  NC_002947/1355327‑1355622

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