Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,880,851 T→A 8.1% intergenic (‑92/+35) leuE ← / ← dmlR leucine efflux protein/transcriptional activator of dmlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,880,8510TA 8.1% 292.3 / 13.1 86intergenic (‑92/+35)leuE/dmlRleucine efflux protein/transcriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base T (41/38);  new base A (3/4);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCC  >  NC_000913/1880707‑1880998
                                                                                                                                                |                                                                                                                                                   
aTGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAAc                                                                                                                                                 <  4:340667/149‑1 (MQ=255)
    ccccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGa                                                                                                                                             <  3:330864/149‑1 (MQ=255)
    ccccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGa                                                                                                                                             >  4:251962/1‑149 (MQ=255)
    ccccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTGCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTGGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGa                                                                                                                                             <  4:855639/149‑1 (MQ=255)
     cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGat                                                                                                                                            <  3:865389/149‑1 (MQ=255)
      ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGatt                                                                                                                                           <  3:501948/149‑2 (MQ=255)
        aaCCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGatttt                                                                                                                                         <  3:279201/149‑4 (MQ=255)
         aCCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGata                                                                                                                                           <  4:830691/146‑1 (MQ=255)
         aCCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGata                                                                                                                                           >  3:830691/1‑146 (MQ=255)
          ccAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCa                                                                                                                                        <  1:544202/148‑1 (MQ=255)
           cAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATCCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCaa                                                                                                                                       <  2:745150/148‑1 (MQ=255)
                   tCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAAt                                                                                                                              >  3:61999/1‑149 (MQ=255)
                    ccAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAAt                                                                                                                              >  2:381154/1‑148 (MQ=255)
                    ccAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCCCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAAt                                                                                                                              <  2:783731/148‑1 (MQ=255)
                      aGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTg                                                                                                                            <  2:636712/148‑1 (MQ=255)
                       gTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGa                                                                                                                           <  1:836961/148‑1 (MQ=255)
                        tAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAAt                                                                                                                         <  3:251962/149‑1 (MQ=255)
                           ttCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                      >  3:585542/1‑149 (MQ=255)
                           ttCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                      >  3:149292/1‑149 (MQ=255)
                           ttCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                      >  3:144832/1‑149 (MQ=255)
                            tCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTACTTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                      >  1:303783/1‑148 (MQ=255)
                               gAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTaa                                                                                                                   >  2:313331/1‑148 (MQ=255)
                               gAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAAt                                                                                                                  >  3:642089/1‑149 (MQ=255)
                                   cccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAgg                                                                                                              >  4:684887/1‑149 (MQ=255)
                                    ccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAgg                                                                                                              >  1:366151/1‑148 (MQ=255)
                                    ccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGc                                                                                                             <  4:149292/149‑1 (MQ=255)
                                     cGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGc                                                                                                             >  1:219234/1‑148 (MQ=255)
                                            aGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCa                                                                                                     <  4:404212/149‑1 (MQ=255)
                                                aaCACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTg                                                                                                  >  2:332412/1‑148 (MQ=255)
                                                      ttAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAAt                                                                                            >  1:341852/1‑148 (MQ=255)
                                                       tAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTc                                                                                          <  4:144832/149‑1 (MQ=255)
                                                         aCCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCgg                                                                                        >  4:508959/1‑149 (MQ=255)
                                                             ctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCa                                                                                                                <  2:358907/121‑1 (MQ=255)
                                                             ctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCa                                                                                                                >  1:358907/1‑121 (MQ=255)
                                                             ctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTa                                                                                     >  2:261360/1‑148 (MQ=255)
                                                               ttAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGt                                                                                   >  2:692238/1‑148 (MQ=255)
                                                                       cTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCaa                                                                          >  4:654482/1‑149 (MQ=255)
                                                                        tttCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCaaa                                                                         >  4:440610/1‑149 (MQ=255)
                                                                        tttCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCaaa                                                                         >  3:880044/1‑149 (MQ=255)
                                                                        tttCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCaaa                                                                         >  4:545287/1‑149 (MQ=255)
                                                                         ttCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCaaaa                                                                        <  3:807593/149‑1 (MQ=255)
                                                                             aTCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACt                                                                    <  4:58925/149‑1 (MQ=255)
                                                                             aTCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACt                                                                    <  3:684887/149‑1 (MQ=255)
                                                                              tCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCa                                                                                                                <  1:542538/104‑1 (MQ=255)
                                                                              tCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCa                                                                                                                >  2:542538/1‑104 (MQ=255)
                                                                                       aTTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTaaaa                                                           <  2:715936/148‑1 (MQ=255)
                                                                                       aTTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTaaaa                                                           >  2:48999/1‑148 (MQ=255)
                                                                                            aGGGGAATATGCTAACCAAGCGGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTa                                                     >  4:456397/1‑149 (MQ=255)
                                                                                             ggTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTAc                                                    <  4:61999/149‑1 (MQ=255)
                                                                                              gTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAGATACTGGCGTAAAAATTCTAcc                                                   >  3:554807/1‑149 (MQ=255)
                                                                                               tGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCg                                                  >  3:309867/1‑149 (MQ=255)
                                                                                                   tatGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTa                                               <  1:663421/148‑1 (MQ=255)
                                                                                                    atGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTAt                                              <  2:341852/148‑1 (MQ=255)
                                                                                                    atGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATc                                             <  3:363735/149‑1 (MQ=255)
                                                                                                    atGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATc                                             <  3:545287/149‑1 (MQ=255)
                                                                                                           ccAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACt                                       >  1:494627/1‑148 (MQ=255)
                                                                                                            cAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGc                                          <  1:163785/144‑1 (MQ=255)
                                                                                                            cAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGc                                          >  2:163785/1‑144 (MQ=255)
                                                                                                             aaTCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACttt                                     >  1:724792/1‑148 (MQ=255)
                                                                                                                cTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGc                                  <  2:366151/148‑1 (MQ=255)
                                                                                                                  gTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGa                               <  4:309867/149‑1 (MQ=255)
                                                                                                                   tAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCg                                                                                         <  2:79466/90‑1 (MQ=255)
                                                                                                                   tAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCg                                                                                         >  1:79466/1‑90 (MQ=255)
                                                                                                                      cGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGt                            >  2:638532/1‑148 (MQ=255)
                                                                                                                         cTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAACAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGc                         >  2:628341/1‑148 (MQ=255)
                                                                                                                         cTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGcc                        >  4:159626/1‑149 (MQ=255)
                                                                                                                               aaGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCtt                  >  4:561350/1‑149 (MQ=255)
                                                                                                                               aaGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCtt                  <  4:880044/149‑1 (MQ=255)
                                                                                                                                   agaAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGaa               >  1:878906/1‑148 (MQ=255)
                                                                                                                                   agaAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGaa               <  2:834165/148‑1 (MQ=255)
                                                                                                                                   agaAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGaa               <  2:724792/148‑1 (MQ=255)
                                                                                                                                    gaAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAAc             >  4:583357/1‑149 (MQ=255)
                                                                                                                                     aaaaTCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACa            <  3:440610/149‑1 (MQ=255)
                                                                                                                                       aaTCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTa                                                     >  3:434166/1‑106 (MQ=255)
                                                                                                                                       aaTCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTa                                                     <  4:434166/106‑1 (MQ=255)
                                                                                                                                       aaTCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACAt           <  1:692238/148‑1 (MQ=255)
                                                                                                                                        aTCAGGTTAAAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCg                                   >  2:453977/1‑123 (MQ=255)
                                                                                                                                        aTCAGGTTAAAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCg                                   <  1:453977/123‑1 (MQ=255)
                                                                                                                                         tCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATaaa        >  3:865547/1‑149 (MQ=255)
                                                                                                                                         tCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATaaa        <  3:456397/149‑1 (MQ=255)
                                                                                                                                         tCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATaa         <  2:219234/148‑1 (MQ=255)
                                                                                                                                          cAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATaaa        >  2:69438/1‑148 (MQ=255)
                                                                                                                                          cAGGTTAAAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATAc                                                                     >  1:326420/1‑87 (MQ=255)
                                                                                                                                          cAGGTTAAAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTg                                <  1:761958/124‑1 (MQ=255)
                                                                                                                                          cAGGTTAAAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTg                                >  2:761958/1‑124 (MQ=255)
                                                                                                                                            ggTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACgg     <  3:654482/149‑1 (MQ=255)
                                                                                                                                              ttttAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGc    <  2:303783/148‑1 (MQ=255)
                                                                                                                                               tttAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGccc  <  3:508959/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCC  >  NC_000913/1880707‑1880998

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