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.100.0% 245.2 / NA 60intergenic (‑283/‑262)trxB/lrpthioredoxin reductase, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (29/31);  total (29/31)

GCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCT  >  NC_000913/932195‑932465
                                                                                                                                          |                                                                                                                                    
gCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaaga                                                                                                                                      >  1:829972/1‑139 (MQ=255)
  tACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaa                                                                                                                                    >  2:1734071/1‑139 (MQ=255)
   aCTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaa                                                                                                                                   >  2:1479253/1‑139 (MQ=255)
    cTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaa                                                                                                                                  >  1:303270/1‑139 (MQ=255)
      tAAGTTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  1:285351/139‑1 (MQ=255)
      tAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  2:1503887/139‑1 (MQ=255)
        aaTTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  1:1137894/139‑3 (MQ=255)
        aaTTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  2:104189/139‑3 (MQ=255)
        aaTTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  2:193972/139‑3 (MQ=255)
          tttCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacagt                                                                                                                            <  2:303270/139‑5 (MQ=255)
              ccGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTc                                                                                                                       >  2:498698/1‑139 (MQ=255)
               cGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCt                                                                                                                      >  1:368273/1‑139 (MQ=255)
                  tttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTtat                                                                                                                   <  1:307867/139‑1 (MQ=255)
                  tttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTtat                                                                                                                   <  1:1734071/139‑1 (MQ=255)
                           cGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATa                                                                                                          >  1:192282/1‑139 (MQ=255)
                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  1:251335/115‑1 (MQ=255)
                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                >  2:251335/1‑115 (MQ=255)
                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTtata                                                                                                                  <  2:1264527/128‑1 (MQ=255)
                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTtata                                                                                                                  >  1:1264527/1‑128 (MQ=255)
                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                                       >  1:880428/1‑139 (MQ=255)
                               tAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                      <  2:368273/139‑1 (MQ=255)
                                 aaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                                    <  1:1479253/139‑1 (MQ=255)
                                      tgaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGt                                                                                               <  1:840354/139‑1 (MQ=255)
                                        aaTAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                             <  2:880428/139‑1 (MQ=255)
                                        aaTAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTACGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                             >  1:1088773/1‑139 (MQ=255)
                                          tAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAAt                                                                                           <  2:1088773/139‑1 (MQ=255)
                                              cGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCa                                                                                       <  2:1592651/139‑1 (MQ=255)
                                                ccccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacagt                                                                                                                            >  2:582843/1‑97 (MQ=255)
                                                ccccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacagt                                                                                                                            <  1:582843/101‑5 (MQ=255)
                                                   cTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         <  2:862250/132‑1 (MQ=255)
                                                   cTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         >  1:862250/1‑132 (MQ=255)
                                                    tGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATgcg                                                                                                              >  1:1010213/1‑110 (MQ=255)
                                                    tGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATgcg                                                                                                              <  2:1010213/110‑1 (MQ=255)
                                                              tatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                                       >  1:810640/1‑139 (MQ=255)
                                                               atCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  2:899345/82‑1 (MQ=255)
                                                               atCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                >  1:899345/1‑82 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  2:1706930/81‑1 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                >  1:1706930/1‑81 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTtat                                                                                                                   <  2:1420961/93‑1 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTtat                                                                                                                   >  1:1420961/1‑93 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                 <  2:627184/127‑1 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                 >  1:627184/1‑127 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                     >  2:635389/1‑139 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                     >  2:631280/1‑139 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                     >  1:666123/1‑139 (MQ=255)
                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                     >  2:485500/1‑139 (MQ=255)
                                                                   ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              >  2:168674/1‑78 (MQ=255)
                                                                   ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  1:168674/80‑3 (MQ=255)
                                                                   ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         >  2:962440/1‑116 (MQ=255)
                                                                   ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         <  1:962440/116‑1 (MQ=255)
                                                                   ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTg                                                                            >  2:601114/1‑129 (MQ=255)
                                                                   ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTg                                                                            <  1:601114/129‑1 (MQ=255)
                                                                    gCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaa                                                                                                                                   <  1:373341/74‑1 (MQ=255)
                                                                    gCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaa                                                                                                                                   >  2:373341/1‑74 (MQ=255)
                                                                          tgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTa                                                           <  2:810640/139‑1 (MQ=255)
                                                                             aCTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTT‑‑GAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCg                                                      <  1:1345488/139‑1 (MQ=255)
                                                                              cTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         >  1:1754914/1‑105 (MQ=255)
                                                                              cTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         <  2:1754914/105‑1 (MQ=255)
                                                                                   aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  1:1613189/62‑1 (MQ=255)
                                                                                   aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                >  2:1613189/1‑62 (MQ=255)
                                                                                      cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACTGTGCTAGt                                                                       >  2:659154/1‑115 (MQ=255)
                                                                                      cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACTGTGCTAGt                                                                       <  1:659154/115‑1 (MQ=255)
                                                                                       gATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGa                                              <  1:1553989/139‑1 (MQ=255)
                                                                                        aTGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              >  2:892668/1‑57 (MQ=255)
                                                                                        aTGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  1:892668/59‑3 (MQ=255)
                                                                                         tGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                     <  1:624499/114‑1 (MQ=255)
                                                                                         tGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                     >  2:624499/1‑114 (MQ=255)
                                                                                                   tGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACcat                                                                                                     <  2:848864/72‑1 (MQ=255)
                                                                                                   tGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACcat                                                                                                     >  1:848864/1‑72 (MQ=255)
                                                                                                             cTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCa                                                                                                            >  1:685841/1‑55 (MQ=255)
                                                                                                             cTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCa                                                                                                            <  2:685841/55‑1 (MQ=255)
                                                                                                               ggTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTg                      >  2:847984/1‑139 (MQ=255)
                                                                                                                gtgtTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGc                     <  1:631280/139‑1 (MQ=255)
                                                                                                                   ttttGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTc                              >  1:464902/1‑127 (MQ=255)
                                                                                                                   ttttGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTc                              <  2:464902/127‑1 (MQ=255)
                                                                                                                   ttttGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCcac                  >  2:1360322/1‑139 (MQ=255)
                                                                                                                     ttGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACaa                >  1:1234310/1‑139 (MQ=255)
                                                                                                                          aaaaCATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCaa           >  1:752286/1‑139 (MQ=255)
                                                                                                                            aaCATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         >  1:1573352/1‑59 (MQ=255)
                                                                                                                            aaCATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         <  2:1573352/59‑1 (MQ=255)
                                                                                                                            aaCATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACg         <  2:1234310/139‑1 (MQ=255)
                                                                                                                             aCATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACgg        <  1:1360322/139‑1 (MQ=255)
                                                                                                                                  cGAGGAAGAAAAAAACAGTATTCTTATATGCGCATACCCATGCATGTAAATACCATGATTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTc                                      >  1:1482359/1‑104 (MQ=255)
                                                                                                                                  cGAGGAAGAAAAAAACAGTATGCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTc                                      <  2:1482359/104‑1 (MQ=255)
                                                                                                                                   gAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGtcttct  <  2:192282/139‑1 (MQ=255)
                                                                                                                                   gAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGtcttct  <  1:847984/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                    
GCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCT  >  NC_000913/932195‑932465

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: