Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 2,600,577 G→T intergenic (‑139/+728) SAUSA300_RS13365 ← / ← SAUSA300_RS13370 diaminopimelate epimerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077932,600,5770GT98.9% 275.5 / ‑4.6 90intergenic (‑139/+728)SAUSA300_RS13365/SAUSA300_RS13370diaminopimelate epimerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (6/83);  total (6/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTAGTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGTAATTTTTT  >  NC_007793/2600479‑2600673
                                                                                                  |                                                                                                
ttAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCtagg                                                                                                 <  1:1537716/100‑2 (MQ=255)
 tAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCtaggg                                                                                                <  1:1361232/100‑3 (MQ=255)
              ctctTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCt                                                                                   <  1:396217/100‑1 (MQ=25)
              ctctTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCt                                                                                   <  1:3888700/100‑1 (MQ=25)
                 ttATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGcggtg                                                                                <  1:4477538/100‑5 (MQ=255)
                      cAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCa                                                                           <  1:992227/100‑1 (MQ=255)
                       aGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCag                                                                          <  1:4623600/100‑2 (MQ=255)
                       aGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCaa                                                                          <  1:1601179/100‑1 (MQ=255)
                       aGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCaa                                                                          <  1:3555906/100‑1 (MQ=255)
                       aGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCaa                                                                          <  1:4271202/100‑1 (MQ=255)
                                      aacaacTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATgg                                                           <  1:630553/100‑1 (MQ=255)
                                      aacaacTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATga                                                           <  1:3226759/100‑2 (MQ=255)
                                       acaacTGTATTCCTCTAACCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATggt                                                          <  1:5434337/100‑2 (MQ=255)
                                       acaacTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCAtcta                                                          <  1:3495674/100‑4 (MQ=255)
                                       acaacTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATggt                                                          <  1:1926750/100‑2 (MQ=255)
                                       acaacTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATggg                                                          <  1:506683/100‑2 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAATCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:400099/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:3177225/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:4553974/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:4851536/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:4868140/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:488859/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:4896857/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:5046885/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:5175442/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:178636/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:783000/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTACGc                                                  <  1:4519456/100‑1 (MQ=255)
                                               aTTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTATATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGc                                                  <  1:1191942/100‑1 (MQ=255)
                                                 tCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAg                                                <  1:4306383/100‑1 (MQ=255)
                                                 tCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAg                                                <  1:2851816/100‑1 (MQ=255)
                                                 cgatctAAACCTGTAATTTTTTAGAAGATAAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAg                                                <  1:3886843/97‑1 (MQ=255)
                                                  cgatctaaCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:19141/94‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGACCCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:1537388/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:2453557/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:2301604/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:3530544/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:1636618/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:2150646/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:4910609/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:5346032/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:3110423/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:1482050/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgt                                               <  1:5391053/100‑1 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgc                                               <  1:887911/100‑2 (MQ=255)
                                                  ccTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAga                                               <  1:3766221/100‑2 (MQ=255)
                                                   ctctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgtt                                              <  1:3525656/100‑1 (MQ=255)
                                                   ctctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgtt                                              <  1:4712706/100‑1 (MQ=255)
                                                   ctctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgtt                                              <  1:1685986/100‑1 (MQ=255)
                                                   ctctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgtt                                              <  1:3002531/100‑1 (MQ=255)
                                                   ctctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgtt                                              <  1:5517122/100‑1 (MQ=255)
                                                   ctctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTAACTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgta                                              <  1:5268018/100‑2 (MQ=255)
                                                   atctAAGCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgtt                                              <  1:4793751/99‑1 (MQ=255)
                                                    tctAGACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgttg                                             <  1:3541588/100‑1 (MQ=255)
                                                    tctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgttg                                             <  1:1576768/100‑1 (MQ=255)
                                                    tctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgttg                                             <  1:2624492/100‑1 (MQ=255)
                                                    tctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgtta                                             <  1:5219008/100‑2 (MQ=255)
                                                    tctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgtta                                             <  1:220998/100‑2 (MQ=255)
                                                    tctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAggat                                             <  1:2512432/100‑4 (MQ=255)
                                                    atctcgccTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAgttg                                             <  1:3783482/94‑1 (MQ=255)
                                                     ctAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAGAACCATGGCTCTTATGCAgttgt                                            <  1:5424821/100‑1 (MQ=255)
                                                       atgccTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGttgttg                                          >  1:4261215/4‑100 (MQ=255)
                                                       aaaCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGttgttg                                          >  1:4538077/1‑100 (MQ=255)
                                                       aaaCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGttgttg                                          >  1:2256635/1‑100 (MQ=255)
                                                       aaaCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGttgttg                                          >  1:887731/1‑100 (MQ=255)
                                                       aaaCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGttgttg                                          <  1:943245/100‑1 (MQ=255)
                                                       aaaCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGttgttg                                          >  1:1440627/1‑100 (MQ=255)
                                                       aaaCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGttgttg                                          <  1:2177549/100‑1 (MQ=255)
                                                       aaaCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGttgttg                                          >  1:4085888/1‑100 (MQ=255)
                                                       aaaCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGttgttg                                          <  1:282077/100‑1 (MQ=255)
                                                        aaCCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGc                                         <  1:3785074/100‑1 (MQ=255)
                                                         atcTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGca                                        <  1:3227833/98‑2 (MQ=255)
                                                         atcTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCg                                        <  1:3405046/98‑1 (MQ=255)
                                                          ccTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGttgttgggt                                       <  1:2660811/100‑4 (MQ=255)
                                                          ccTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGa                                       <  1:2087759/100‑1 (MQ=255)
                                                                 ttttttAGAAGATGAACCTAATAGCTCTATCTATTACGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGaagcaagc                                <  1:2552366/100‑1 (MQ=255)
                                                                   ttttAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGaagccgggtg                              <  1:2647076/100‑7 (MQ=255)
                                                                      tAGAAGATGAACCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGaagcaagggggtt                           <  1:2076514/100‑7 (MQ=255)
                                                                               aaCCTAATAGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCctct                  <  1:2042644/100‑1 (MQ=255)
                                                                                aCCTAATAGCTCTATCTATTAGGTAAGTCATGCCAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTa                 <  1:321256/100‑1 (MQ=255)
                                                                                       aGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGt          <  1:622511/100‑1 (MQ=255)
                                                                                       aGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGt          <  1:973653/100‑1 (MQ=255)
                                                                                       aGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGt          <  1:4503559/100‑1 (MQ=255)
                                                                                       aGCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGt          <  1:2120784/100‑1 (MQ=255)
                                                                                        gCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGTa         <  1:1828766/100‑1 (MQ=255)
                                                                                        gCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGTa         <  1:218607/100‑1 (MQ=255)
                                                                                        gCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGTa         <  1:3848883/100‑1 (MQ=255)
                                                                                        gCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGTa         <  1:3501760/100‑1 (MQ=255)
                                                                                        gCTCTATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGTa         <  1:2729568/100‑1 (MQ=255)
                                                                                        gCTCTATCTAGTAGGTAAGTCATGCCAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGTa         <  1:3605030/100‑1 (MQ=255)
                                                                                         ctctATCTATTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGtcg        <  1:5003334/100‑3 (MQ=255)
                                                                                               ctaTTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGTAAtttttt  <  1:4814766/100‑1 (MQ=255)
                                                                                                  |                                                                                                
TTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTCTATCTAGTAGGTAAGTCATGCTAATTCGCAAAAGCGAAAACCATGGCTCTTATGCAGTTGTTGCGAAGCAAGCAACTGTATTCCTCTAAAACTGTAATTTTTT  >  NC_007793/2600479‑2600673

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