Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,799,214 A→G 16.5% intergenic (‑185/+31) PP_2452 ← / ← ansB hypothetical protein/glutaminase‑asparaginase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,799,2140AG16.5% 191.5 / 27.1 84intergenic (‑185/+31)PP_2452/ansBhypothetical protein/glutaminase‑asparaginase
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base G (7/7);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CTGTAAGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTC  >  NC_002947/2799072‑2799358
                                                                                                                                               |                                                                                                                                                
cTGTAAGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGgcg                                                                                                                                           >  2:1126090/1‑150 (MQ=255)
      gtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCgg                                                                                                                                     >  2:1142404/1‑150 (MQ=255)
            tCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGAAGCGCTTGGGTTGTAACCTCCTGGTGTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGa                                                                                                                               <  1:1206610/150‑1 (MQ=255)
             cGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAg                                                                                                                              <  1:1049977/150‑1 (MQ=255)
             cGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAg                                                                                                                              <  1:1757939/150‑1 (MQ=255)
              gCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAgg                                                                                                                             <  1:1126090/150‑1 (MQ=255)
              gCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAgg                                                                                                                             <  2:1812871/150‑1 (MQ=255)
               cAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAggg                                                                                                                            >  2:1210305/1‑150 (MQ=255)
                aaCAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAgggg                                                                                                                           <  1:1729154/150‑1 (MQ=255)
                    gCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAg                                                                                                                              <  1:1641136/143‑1 (MQ=255)
                    gCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAg                                                                                                                              >  2:1641136/1‑143 (MQ=255)
                      ttGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCa                                                                                                                     >  2:704372/1‑150 (MQ=255)
                          cAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATGGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCa                                                                                                                 <  1:598869/150‑1 (MQ=255)
                            ggAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTCCGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAg                                                                                                                              <  1:1074106/135‑1 (MQ=255)
                            ggAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAg                                                                                                                              >  2:1074106/1‑135 (MQ=255)
                             gagaAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTa                                                                                                              <  1:704372/150‑1 (MQ=255)
                                 aaTCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTcc                                                                                                          <  1:1210305/150‑1 (MQ=255)
                                 aaTCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTcc                                                                                                          <  2:564507/150‑1 (MQ=255)
                                   tCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTATACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCa                                                                                                        >  2:707674/1‑150 (MQ=255)
                                    ccAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCag                                                                                                       >  1:1051110/1‑150 (MQ=255)
                                             gCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGc                                                                                              <  1:219995/150‑1 (MQ=255)
                                             gCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGc                                                                                              >  2:1617872/1‑150 (MQ=255)
                                              ccTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACTCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAg                                                                                                                              <  1:1459005/117‑1 (MQ=255)
                                              ccTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACTCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAg                                                                                                                              >  2:1459005/1‑117 (MQ=255)
                                               cTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTg                                                                                            >  2:372199/1‑150 (MQ=255)
                                                 gTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCa                                                                                          <  1:1142404/150‑1 (MQ=255)
                                                  tAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAg                                                                                         >  1:838046/1‑150 (MQ=255)
                                                    aaCCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGct                                                                                       >  2:1299357/1‑150 (MQ=255)
                                                         tGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTg                                                                                  <  2:1051110/150‑1 (MQ=255)
                                                         tGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTg                                                                                  <  1:1299357/150‑1 (MQ=255)
                                                         tGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACGCAGTACTCCCAGAAGATCCGCTGCAGCTAttta                                                                                  >  2:804061/1‑149 (MQ=255)
                                                                ggTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCt                                                                           >  1:1128759/1‑150 (MQ=255)
                                                                  tCATAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTgg                                                                         <  1:675067/150‑1 (MQ=255)
                                                                  tCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTgg                                                                         <  2:1128759/150‑1 (MQ=255)
                                                                      aaaGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGATGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGctct                                                                                                                                 >  2:1479796/1‑87 (MQ=37)
                                                                      aaaGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGATGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGctct                                                                                                                                 <  1:1479796/90‑1 (MQ=37)
                                                                      aaaGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCt                                                                     >  1:1777025/1‑150 (MQ=255)
                                                                        aGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTg                                                                   <  2:1555783/150‑1 (MQ=255)
                                                                          gTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGCGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCAGCTGCAGCTCTTTTCTGTCCTGGGTCTTGGt                                                                 >  1:823603/1‑150 (MQ=255)
                                                                           tAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTg                                                                                                                                <  1:1129462/86‑1 (MQ=255)
                                                                           tAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTg                                                                                                                                >  2:1129462/1‑86 (MQ=255)
                                                                             aaCGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCAt                                                              >  1:561037/1‑150 (MQ=255)
                                                                                  tGTCGCCGGTTTGTATCCTTTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTa                                                         <  2:1414547/150‑1 (MQ=255)
                                                                                  tGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTcc                                                                                                          <  2:90004/101‑1 (MQ=255)
                                                                                  tGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTcc                                                                                                          >  1:90004/1‑101 (MQ=255)
                                                                                    tCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTggg                                                                        <  2:1662215/133‑1 (MQ=255)
                                                                                    tCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTggg                                                                        >  1:1662215/1‑133 (MQ=255)
                                                                                       ccGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAgg                                                                                                                             <  2:194747/77‑1 (MQ=255)
                                                                                       ccGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAgg                                                                                                                             >  1:194747/1‑77 (MQ=255)
                                                                                          ggTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTg                                                                                  <  1:688673/118‑1 (MQ=255)
                                                                                          ggTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTg                                                                                  >  2:688673/1‑118 (MQ=255)
                                                                                           gTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGcc                                                 >  1:1533341/1‑150 (MQ=255)
                                                                                            tttGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCAt                                                              <  2:513376/136‑1 (MQ=255)
                                                                                            tttGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTACCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCAt                                                              >  1:513376/1‑136 (MQ=255)
                                                                                                tAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGAt                                            >  1:865931/1‑150 (MQ=255)
                                                                                                    ctctGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCcgc                                        <  1:1617872/150‑1 (MQ=255)
                                                                                                         tgtCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGc                                                                                              <  2:116151/91‑1 (MQ=255)
                                                                                                         tgtCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGc                                                                                              >  1:116151/1‑91 (MQ=255)
                                                                                                              aCGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCa                                                                                                                 >  2:798519/1‑67 (MQ=255)
                                                                                                              aCGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCa                                                                                                                 <  1:798519/67‑1 (MQ=255)
                                                                                                                      gAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAgg                      <  1:707674/150‑1 (MQ=255)
                                                                                                                       aGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGt                     >  1:169045/1‑150 (MQ=255)
                                                                                                                       aGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGt                     >  1:1308311/1‑150 (MQ=255)
                                                                                                                        ggAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGc                               <  2:21881/139‑1 (MQ=255)
                                                                                                                        ggAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGc                               >  1:21881/1‑139 (MQ=255)
                                                                                                                         gAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGTCAGGATcc                                          <  2:271841/127‑1 (MQ=255)
                                                                                                                         gAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGTCAGGATcc                                          >  1:271841/1‑127 (MQ=255)
                                                                                                                         gAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCg                   >  1:519441/1‑150 (MQ=255)
                                                                                                                          aaTTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGt                  <  2:865931/150‑1 (MQ=255)
                                                                                                                            ttGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                <  2:838046/150‑1 (MQ=255)
                                                                                                                                cAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCa            <  2:519441/150‑1 (MQ=255)
                                                                                                                                    ttcccccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTgg                                                                         >  1:1409673/5‑85 (MQ=25)
                                                                                                                                    ttcccccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTgg                                                                         <  2:1409673/81‑1 (MQ=25)
                                                                                                                                      ggCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTg                                                                   <  1:1294400/89‑1 (MQ=255)
                                                                                                                                      ggCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTg                                                                   >  2:1294400/1‑89 (MQ=255)
                                                                                                                                        cccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTgg                                                                  >  1:1444705/1‑88 (MQ=37)
                                                                                                                                        cccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTgg                                                                  <  2:1444705/88‑1 (MQ=37)
                                                                                                                                        cccTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCgg                             >  1:1861579/1‑125 (MQ=255)
                                                                                                                                        cccTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCgg                             <  2:1861579/125‑1 (MQ=255)
                                                                                                                                         ccTCAGAGCCGGGGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGt   >  2:1050705/1‑150 (MQ=255)
                                                                                                                                          cTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                >  2:497550/1‑136 (MQ=255)
                                                                                                                                          cTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                <  1:497550/136‑1 (MQ=255)
                                                                                                                                          cTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGGTCAGGTCGTGAGCCACGACCCAGTc  >  2:1770979/1‑150 (MQ=255)
                                                                                                                                          cTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTc  >  1:588673/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CTGTAAGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTC  >  NC_002947/2799072‑2799358

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