Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,798,262 A→G 8.9% intergenic (‑307/+67) gltB ← / ← PP_5077 L‑glutamate synthase subunit alpha/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,798,2620AG 8.9% 189.2 / 7.3 68intergenic (‑307/+67)gltB/PP_5077L‑glutamate synthase subunit alpha/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/24);  new base G (3/3);  total (41/27)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

TGAAAATTTATGTCGCAATTTTGAATGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTC  >  NC_002947/5798116‑5798401
                                                                                                                                                   |                                                                                                                                           
tGAAAATTTATGTCGCAATTTTGGATGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTATTAGCGTTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCgcagca                                                                                                                                          <  1:329151/150‑1 (MQ=255)
tGAAAATTTATGTCGCAATTTTGAATGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCTTGCCGCGCTCGGGCgcagca                                                                                                                                          <  2:376653/150‑1 (MQ=255)
tGAAAATTTATGTCGCAATTTTGAATGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCTGCGCCCGAGCgcggca                                                                                                                                          <  1:172791/150‑1 (MQ=255)
 gAAAATTTATGTCGCAATTTTGAATGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGcagcag                                                                                                                                         >  1:582186/1‑150 (MQ=255)
            tCGCAATTTTGAATGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGaaa                                                                                                                              >  1:1112376/1‑150 (MQ=255)
             cGCAATTTTGAATGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCGGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAg                                                                                                                             <  1:1288362/150‑1 (MQ=255)
              gCAATTTTGAATGTTTTTGCGCCATGCGCAGGGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCGCGGGCACAGAAGGCTAATCCAAAGc                                                                                                                            >  2:479331/1‑150 (MQ=255)
               cAATTTTGAATGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCt                                                                                                                           >  2:1920971/1‑150 (MQ=255)
                   tttGAATGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGaa                                                                                                                       <  2:582186/150‑1 (MQ=255)
                       aaTGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCg                                                                                                                   <  2:1873581/150‑1 (MQ=255)
                         tGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCCGCCGGCTCGTCGAAAGCTAGAATTCGCt                                                                                                                 >  2:572212/1‑150 (MQ=255)
                             tttGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGTCTTGTGGCCGGCAGTCCGGGACAGAAAGGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGc                                                                                                             <  1:1170391/150‑1 (MQ=255)
                                   ccATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTgg                                                                                                       >  1:755606/1‑150 (MQ=255)
                                          ccGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGcggcgg                                                                                                <  1:370987/150‑1 (MQ=255)
                                           cGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTc                                                                                                                                  >  2:806667/1‑115 (MQ=255)
                                           cGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTc                                                                                                                                  <  1:806667/115‑1 (MQ=255)
                                                    ggCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATgccg                                                                                      >  1:1649923/1‑150 (MQ=255)
                                                     gCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATgccgc                                                                                     >  2:908494/1‑150 (MQ=255)
                                                           cTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGaa                                                                               <  2:1354313/150‑1 (MQ=255)
                                                            tagAGGCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAg                                                                              >  2:469864/3‑150 (MQ=255)
                                                            tCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGAt                                                                                          <  1:1597964/138‑1 (MQ=255)
                                                            tCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGAt                                                                                          >  2:1597964/1‑138 (MQ=255)
                                                            tCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAg                                                                              >  1:777494/1‑150 (MQ=255)
                                                                 tCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATc                                                                         >  2:1774080/1‑150 (MQ=255)
                                                                  cTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCa                                                                        <  2:1393585/150‑1 (MQ=255)
                                                                   tGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAg                                                                       >  1:564479/1‑150 (MQ=255)
                                                                   tGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAg                                                                       >  1:1701230/1‑150 (MQ=255)
                                                                   tGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAg                                                                       <  1:621809/150‑1 (MQ=255)
                                                                    gAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGc                                                                      >  2:579996/1‑150 (MQ=255)
                                                                      aGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTGCAGCGTTCGAGAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCgg                                                                    <  2:529013/150‑1 (MQ=255)
                                                                            ttCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGt                                                               <  1:908494/150‑1 (MQ=255)
                                                                               aGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGc                                                            >  2:1551342/1‑150 (MQ=255)
                                                                                gCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGcc                                                           >  2:201410/1‑150 (MQ=255)
                                                                                gCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGcc                                                           >  2:1750819/1‑150 (MQ=255)
                                                                                 ccATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCa                                                          >  1:1705745/1‑150 (MQ=255)
                                                                                   aTGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCTGCGCCCGAGCGCGGCAGGCtttt                                                                                                                                   <  2:132256/75‑4 (MQ=21)
                                                                                   aTGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCTGCGCCCGAGCGCGGCAGGCtttt                                                                                                                                   >  1:132256/1‑72 (MQ=21)
                                                                                   aTGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGc                                                        <  1:1303159/150‑1 (MQ=255)
                                                                                        tGCAACTACGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCGGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGAGGCCAGCTCttg                                                   >  2:1058488/1‑150 (MQ=255)
                                                                                        tGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttg                                                   >  2:244153/1‑150 (MQ=255)
                                                                                        tGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttg                                                   >  2:746125/1‑150 (MQ=255)
                                                                                         gCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttgt                                                  <  2:1071884/150‑1 (MQ=255)
                                                                                           aaCTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttgttg                                                >  2:1596670/1‑150 (MQ=255)
                                                                                           aaCTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGATGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttgttg                                                >  1:1851370/1‑150 (MQ=255)
                                                                                            aCTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTgg                                               <  2:1649923/150‑1 (MQ=255)
                                                                                              tAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGAc                                             <  2:755606/150‑1 (MQ=255)
                                                                                                 gggCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATgccgc                                                                                     >  1:1041348/1‑107 (MQ=255)
                                                                                                 gggCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATgccgc                                                                                     <  2:1041348/107‑1 (MQ=255)
                                                                                                  ggCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTg                                         >  2:1319195/1‑150 (MQ=255)
                                                                                                    gttgttgCCTGCAGTCCGGGACAGAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCg                                      >  2:320178/2‑150 (MQ=255)
                                                                                                     ttGTGGCCGGCAGTCCGGGGCAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCg                                      >  1:1288466/1‑150 (MQ=255)
                                                                                                     ttGNGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCg                                      >  2:1191029/1‑150 (MQ=255)
                                                                                                          gCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGt                                 >  2:458544/1‑150 (MQ=255)
                                                                                                           ccGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTa                                >  1:1130289/1‑150 (MQ=255)
                                                                                                               cAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTg                            <  1:335952/150‑1 (MQ=255)
                                                                                                                    cGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCaa                       <  2:1701230/150‑1 (MQ=255)
                                                                                                                       gACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCggg                                                                                               >  2:1272413/1‑75 (MQ=255)
                                                                                                                       gACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCggg                                                                                               <  1:1272413/75‑1 (MQ=255)
                                                                                                                           gAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAc                <  1:112148/150‑1 (MQ=255)
                                                                                                                           gAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTAGATTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAc                >  1:135534/1‑150 (MQ=255)
                                                                                                                            aaaaGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGAGATGCGGCGTTGATGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAcc               >  2:1096925/1‑150 (MQ=255)
                                                                                                                             aaaGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCa              <  2:1705745/150‑1 (MQ=255)
                                                                                                                                 ccTGCTGCGCCCGAGCGCGGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCggg                                                                                               <  1:1917195/65‑1 (MQ=25)
                                                                                                                                 ccTGCTGCGCCCGAGCGCGGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCggg                                                                                               >  2:1917195/1‑65 (MQ=25)
                                                                                                                                  ctgctgCGCCCGAGCGCGGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCt          <  1:773154/149‑1 (MQ=255)
                                                                                                                                  ctgctgCGCCCGAGCGCGGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCt          >  2:773154/1‑149 (MQ=255)
                                                                                                                                   tGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTgg        <  1:1596670/150‑1 (MQ=255)
                                                                                                                                    gCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGa       <  1:1135565/150‑1 (MQ=255)
                                                                                                                                     ccGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACAAGCCTGGAc      >  2:921654/1‑150 (MQ=255)
                                                                                                                                       gcgcTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCt    >  2:145644/1‑150 (MQ=255)
                                                                                                                                       gcgcTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCt    >  2:936041/1‑150 (MQ=255)
                                                                                                                                         gcTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTc  >  2:1007819/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
TGAAAATTTATGTCGCAATTTTGAATGTTTTTGCGCCATGCGCCGTGATAGCGGCTTGGCTCGAGTCTGAAGCGTTTCAGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTC  >  NC_002947/5798116‑5798401

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: