Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,060,804 A→G 5.3% I37I (ATT→ATC cbl ← ssuEADCB/tauABCD operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,060,8040AG 5.3% 240.8 / 3.2 76I37I (ATT→ATCcblssuEADCB/tauABCD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (45/27);  new base G (2/2);  total (47/29)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGTCTGCAAGCCGACGAACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAGTTCAGATTTCCCGTTATTGT  >  NC_000913/2060659‑2060950
                                                                                                                                                 |                                                                                                                                                  
ggTCTGCAAGCCGACGAACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGaa                                                                                                                                                   >  3:257421/1‑147 (MQ=255)
  tCTGCAAGCCGACGAACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAatat                                                                                                                                                >  1:338192/1‑148 (MQ=255)
    tGCAAGCCGACGAACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGa                                                                                                                                              >  8:60482/1‑148 (MQ=255)
     gCAAGCCGACGAACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTGCCGAATATGAc                                                                                                                                             >  8:82110/1‑148 (MQ=255)
     gCAAGCCGACGAACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGAc                                                                                                                                             >  4:67446/1‑148 (MQ=255)
      cAAGCCGACGAACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACg                                                                                                                                            <  1:212920/148‑1 (MQ=255)
      cAAGCCGACGAACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACg                                                                                                                                            >  6:86070/1‑148 (MQ=255)
       aaGCCGACGAACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACgg                                                                                                                                           >  2:121239/1‑148 (MQ=255)
            gacgaACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGAc                                                                                                                                      <  8:191491/148‑1 (MQ=255)
             acgaACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGCCg                                                                                                                                     >  2:24304/1‑148 (MQ=255)
               gaACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGcc                                                                                                                                   >  8:334987/1‑148 (MQ=255)
                 aCATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTg                                                                                                                                 <  7:60482/148‑1 (MQ=255)
                       cTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACtgt                                                                                                                            <  1:121239/147‑1 (MQ=255)
                             tcattcaGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGATCAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGCCTGTGAggga                                                                                                                      >  6:109100/1‑145 (MQ=255)
                               attcaGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTaaaa                                                                                                                   <  2:338192/148‑1 (MQ=255)
                                    gAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACtgtg                                                                                                                           >  8:211509/1‑135 (MQ=255)
                                    gAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACtgtg                                                                                                                           <  7:211509/135‑1 (MQ=255)
                                      aTACGTTCTGCAATGACAAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATACTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCatat                                                                                                            >  1:6960/1‑148 (MQ=255)
                                           ttCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGc                                                                                                        >  7:220613/1‑147 (MQ=255)
                                             cTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCaa                                                                                                      <  5:127073/147‑1 (MQ=255)
                                                cAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAgg                                                                                                                        <  3:215960/126‑1 (MQ=255)
                                                cAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAgg                                                                                                                        >  4:215960/1‑126 (MQ=255)
                                                  aTGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTg                                                                                                >  8:364888/1‑148 (MQ=255)
                                                   tGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGGATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGt                                                                                               >  7:356597/1‑148 (MQ=255)
                                                   tGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGGATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGt                                                                                               >  7:361687/1‑148 (MQ=255)
                                                    gACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCat                                                                                                              >  3:279385/1‑132 (MQ=255)
                                                    gACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCat                                                                                                              <  4:279385/132‑1 (MQ=255)
                                                      ccAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAg                                                                                            <  4:102092/148‑1 (MQ=255)
                                                           aaTGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTg                                                                                        >  5:251259/1‑147 (MQ=255)
                                                            aTGCTTTGCCCGGTTCAGTCATTCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTa                                                                                      <  6:251259/148‑1 (MQ=255)
                                                             tGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTaa                                                                                     >  6:218862/1‑148 (MQ=255)
                                                                  tGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATcc                                                                                  <  3:245132/146‑1 (MQ=255)
                                                                    cccGATTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGAc                                                                              >  7:90937/1‑148 (MQ=255)
                                                                      cGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGt                                                                            >  4:70266/1‑148 (MQ=255)
                                                                      cGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGt                                                                            >  4:358489/1‑148 (MQ=255)
                                                                      cGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGt                                                                            >  2:358583/1‑148 (MQ=255)
                                                                      cGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGt                                                                            <  4:257421/148‑1 (MQ=255)
                                                                         ttCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGCCGTg                                                                           >  5:360058/1‑146 (MQ=255)
                                                                              tCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTc                                                                    >  8:128259/1‑148 (MQ=255)
                                                                              tCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTc                                                                    <  3:67446/148‑1 (MQ=255)
                                                                               cATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTcg                                                                   >  7:173885/1‑148 (MQ=255)
                                                                                     cagcagTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGAt                                                              >  1:80968/1‑147 (MQ=255)
                                                                                         agTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCGCGCGGATTATCt                                                         >  4:309106/1‑148 (MQ=255)
                                                                                           tCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATcat                                                        >  6:46080/1‑145 (MQ=255)
                                                                                            cGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTGCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATGtt                                                        >  7:95568/1‑146 (MQ=255)
                                                                                            cGCTTACCTCCTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGATCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTa                                                      <  4:300201/148‑1 (MQ=255)
                                                                                             gCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTa                                                      <  3:70266/147‑1 (MQ=255)
                                                                                              cTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAgt                                                    >  2:260340/1‑148 (MQ=255)
                                                                                                  ccTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAgttgtt                                                >  7:115964/1‑148 (MQ=255)
                                                                                                  ccTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAgttgtt                                                <  7:19719/148‑1 (MQ=255)
                                                                                                  ccTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAgttgtt                                                <  1:260340/148‑1 (MQ=255)
                                                                                                        cgAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGt                                                                                       >  7:184338/1‑103 (MQ=255)
                                                                                                        cgAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGt                                                                                       <  8:184338/103‑1 (MQ=255)
                                                                                                        cgAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAAt                                          >  8:172313/1‑148 (MQ=255)
                                                                                                          aaCAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTc                                        <  8:90937/148‑1 (MQ=255)
                                                                                                            cAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCAc                                      >  8:139605/1‑148 (MQ=255)
                                                                                                             aaaTATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACg                                     >  3:32018/1‑148 (MQ=255)
                                                                                                              aaTATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACgg                                    >  2:359598/1‑148 (MQ=255)
                                                                                                                    tCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAgttgt                                                 >  3:344262/1‑129 (MQ=255)
                                                                                                                    tCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAgttgt                                                 <  4:344262/129‑1 (MQ=255)
                                                                                                                      gATGCCAAGTTCATCCTCCAGTTCCCGGATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTcc                            <  8:356597/148‑1 (MQ=255)
                                                                                                                      gATGCCAAGTTCATCCTCCAGTTCCCGGATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTcc                            <  8:361687/148‑1 (MQ=255)
                                                                                                                      gATGCCAAGTTCATCCTCCAGGTCCCGAATATGACGGCTGACGCCTGCCTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGGTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTcc                            <  5:88555/148‑1 (MQ=255)
                                                                                                                           cAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTa                                                                                      <  2:307597/85‑1 (MQ=255)
                                                                                                                           cAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTa                                                                                      >  1:307597/1‑85 (MQ=255)
                                                                                                                            aaGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAg                      >  7:308524/1‑148 (MQ=255)
                                                                                                                            aaGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAg                      >  2:16506/1‑148 (MQ=255)
                                                                                                                                tCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAGtt                    >  5:64399/1‑146 (MQ=255)
                                                                                                                                  aTCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAGTTCAg                 <  3:358489/147‑1 (MQ=255)
                                                                                                                                        cAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGt                                                                            <  3:144925/82‑1 (MQ=255)
                                                                                                                                        cAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGt                                                                            >  4:144925/1‑82 (MQ=255)
                                                                                                                                         aGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAGTTCAGATTTCCCg         <  8:308524/148‑1 (MQ=255)
                                                                                                                                         aGTTACCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTTTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAAATCACGGTAAACTCCGGGCAGTTCAGATTTccag         >  6:184521/1‑146 (MQ=255)
                                                                                                                                              cagAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAGTTCAGATTTCCCGTTAtt    <  5:190798/146‑1 (MQ=255)
                                                                                                                                               cGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAGTTCAGATTTCCCGTTATTg   >  4:32355/1‑148 (MQ=255)
                                                                                                                                                gAATATGACGGCTGACGCCTGACTGTGAGGTAAATAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGCCGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAGTTCAGATTTCCCGTTATTGt  >  1:360308/1‑148 (MQ=255)
                                                                                                                                                gAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAGTTCAGCTTTCCCGTTATTGt  >  2:256557/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGTCTGCAAGCCGACGAACATTACTGGCTTCATTCAGAATACGTTCTGCAATGACCAGTAATGCTTTGCCCGGTTCAGTCATGCCCAGCAGTCGCTTACCTCGTCGAACAAATATTTCGATGCCAAGTTCATCCTCCAGTTCCCGAATATGACGGCTGACGCCTGACTGTGAGGTAAAAAGCATATTCGCAACCTCTGTCAGGTTGTAATCCTGACGTGCAGCCTCGCGGATTATCTTTAGTTGTTGGAAATTCACGGTAAACTCCGGGCAGTTCAGATTTCCCGTTATTGT  >  NC_000913/2060659‑2060950

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