Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 932,340 (A)8→7 intergenic (‑290/‑255) trxB ← / → lrp thioredoxin reductase, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913932,3330A.96.5% 233.1 / 1.6 58intergenic (‑283/‑262)trxB/lrpthioredoxin reductase, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (27/29);  minor base G (1/1);  total (28/30)
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

GCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAAC  >  NC_000913/932195‑932457
                                                                                                                                          |                                                                                                                            
gCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaaga                                                                                                                              >  1:1237/1‑139 (MQ=255)
gCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaaga                                                                                                                              >  2:609081/1‑139 (MQ=255)
 cTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTGGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaagaa                                                                                                                             >  2:379719/1‑139 (MQ=255)
      tAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                        <  2:1079335/139‑1 (MQ=255)
          tttCCCGCTTTATAAGCCGATTAAATGATGAGTAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacagt                                                                                                                    <  2:1057433/139‑5 (MQ=255)
              ccGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTc                                                                                                               >  1:986989/1‑139 (MQ=255)
              ccGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTc                                                                                                               >  2:1080904/1‑139 (MQ=255)
                  tttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTtat                                                                                                           <  2:1237/139‑1 (MQ=255)
                          ccGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      >  2:1510845/1‑119 (MQ=255)
                          ccGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      <  1:1510845/121‑3 (MQ=255)
                           cGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATa                                                                                                  >  2:1147527/1‑139 (MQ=255)
                             aTTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                >  2:1273416/1‑139 (MQ=255)
                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                               >  1:483822/1‑139 (MQ=255)
                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                               >  2:285898/1‑139 (MQ=255)
                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                               >  2:326618/1‑139 (MQ=255)
                               tAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACca                                                                                              <  1:326618/139‑1 (MQ=255)
                                 aaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                            <  2:1439295/139‑1 (MQ=255)
                                     atgaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACcatgcatg                                                                                        <  2:963825/139‑1 (MQ=255)
                                      tgaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGt                                                                                       <  1:550017/139‑1 (MQ=255)
                                       gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                        >  2:746853/1‑106 (MQ=255)
                                       gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                        <  1:746853/106‑1 (MQ=255)
                                        aaTAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                     >  1:1007904/1‑139 (MQ=255)
                                            aaCGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                 <  1:285898/139‑1 (MQ=255)
                                            aaCGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                 <  1:1273416/139‑1 (MQ=255)
                                               gCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      >  1:1126130/1‑98 (MQ=255)
                                               gCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      <  2:1126130/100‑3 (MQ=255)
                                                 cccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGt                                                                            >  1:63427/1‑139 (MQ=255)
                                                        aaTGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGt                                                                     >  1:519220/1‑139 (MQ=255)
                                                         aTGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                <  2:695292/111‑1 (MQ=255)
                                                         aTGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                >  1:695292/1‑111 (MQ=255)
                                                         aTGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTg                                                                    <  2:63427/139‑1 (MQ=255)
                                                         aTGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAGGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTg                                                                    >  2:1439199/1‑139 (MQ=255)
                                                             atatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACACTCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAg                                                                <  1:49256/139‑1 (MQ=255)
                                                              tatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAGGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                               <  1:1439199/139‑1 (MQ=255)
                                                               atCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTg                                                              <  2:519220/139‑1 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACCTTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTATTGa                                                             <  2:986989/139‑1 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      >  1:1198406/1‑81 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      >  2:174820/1‑81 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      <  1:174820/83‑3 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      <  2:1198406/83‑3 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                            <  1:656536/108‑1 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                            >  2:656536/1‑108 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                             >  1:907729/1‑139 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                             >  1:1283931/1‑139 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                             >  1:478718/1‑139 (MQ=255)
                                                                  tcgcctgttgtACTAAAAATCGATGTTTTGCTTTGAACAGCACCTTGTGTTTTGCGAAAACAGTCGAGGAAGAAAAAAACAGCATTCTTATATGCGCATAATCATGCATGTAAATACCATGTTTACCGTGCTAGTGaaa                                                           <  2:1335707/134‑1 (MQ=255)
                                                                    gCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATc                                                         <  1:1147527/139‑1 (MQ=255)
                                                                    gCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATc                                                         <  2:1491025/139‑1 (MQ=255)
                                                                         ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                    <  2:483822/139‑1 (MQ=255)
                                                                              cTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTAtt                                                                                                                <  2:812502/74‑1 (MQ=255)
                                                                              cTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTAtt                                                                                                                >  1:812502/1‑74 (MQ=255)
                                                                              cTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGc                                               >  1:794976/1‑139 (MQ=255)
                                                                               tAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCg                                              <  2:478718/139‑1 (MQ=255)
                                                                                 aaaaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      >  2:1012087/1‑64 (MQ=255)
                                                                                 aaaaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      <  1:1012087/66‑3 (MQ=255)
                                                                                 aaaaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTg                                            >  2:599093/1‑139 (MQ=255)
                                                                                      cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTAATCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATATACGt                                                    >  2:330875/1‑126 (MQ=255)
                                                                                      cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATACGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                    <  1:330875/126‑1 (MQ=255)
                                                                                       gATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATaa                                                                                                 <  1:67110/80‑1 (MQ=255)
                                                                                       gATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATaa                                                                                                 >  2:67110/1‑80 (MQ=255)
                                                                                         tGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTtat                                                                                                           <  2:1293217/68‑1 (MQ=255)
                                                                                         tGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTtat                                                                                                           >  1:1293217/1‑68 (MQ=255)
                                                                                         tGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACg                                    <  1:599093/139‑1 (MQ=255)
                                                                                             ttGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTAc                                                                        >  1:933276/1‑99 (MQ=255)
                                                                                             ttGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTAc                                                                        <  2:933276/99‑1 (MQ=255)
                                                                                              tGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCAtt                               >  1:241944/1‑139 (MQ=255)
                                                                                                  ttGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTg                           <  1:609081/139‑1 (MQ=255)
                                                                                                   tGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGa                          >  2:889203/1‑139 (MQ=255)
                                                                                                     aCAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATg                        <  2:1007904/139‑1 (MQ=255)
                                                                                                               ggTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGt                                                                                       >  2:1363621/1‑66 (MQ=255)
                                                                                                               ggTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGt                                                                                       <  1:1363621/66‑1 (MQ=255)
                                                                                                               ggTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                     >  2:1216659/1‑68 (MQ=255)
                                                                                                               ggTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                     <  1:1216659/68‑1 (MQ=255)
                                                                                                               ggTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTg              >  1:1293723/1‑139 (MQ=255)
                                                                                                                           aaaCATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAAc  <  2:794976/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                            
GCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAAC  >  NC_000913/932195‑932457

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