Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,880,852 T→A 7.2% intergenic (‑93/+34) 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,8520TA 7.2% 241.1 / 8.7 70intergenic (‑93/+34)leuE/dmlRleucine efflux protein/transcriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base A (3/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

AAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCC  >  NC_000913/1880704‑1880997
                                                                                                                                                    |                                                                                                                                                 
aaaaTGGCCCCAACCAGATAGGTCCATTAATTCAGAACCCCGTATTCAGCGAACACATTAGCCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGtttt                                                                                                                                                   <  3:288135/149‑1 (MQ=255)
aaaaTGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGtttt                                                                                                                                                   <  1:448666/149‑1 (MQ=255)
   aTGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAAc                                                                                                                                                >  3:409699/1‑149 (MQ=255)
   aTGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAAc                                                                                                                                                <  1:273924/149‑1 (MQ=255)
        cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGat                                                                                                                                           >  4:150153/1‑149 (MQ=255)
        cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGat                                                                                                                                           >  3:94355/1‑149 (MQ=255)
        cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGat                                                                                                                                           >  2:548404/1‑149 (MQ=255)
        cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGat                                                                                                                                           >  1:135749/1‑149 (MQ=255)
         ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGatt                                                                                                                                          <  2:135749/149‑2 (MQ=255)
           aaCCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATc                                                                                                                                        >  2:209399/1‑149 (MQ=255)
             ccAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCaa                                                                                                                                      <  1:548404/149‑1 (MQ=255)
               aGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAcc                                                                                                                                    >  1:795617/1‑149 (MQ=255)
                  tAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGa                                                                                                                                 <  3:130399/149‑1 (MQ=255)
                   aGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGAt                                                                                                                                <  2:795617/149‑1 (MQ=255)
                   aGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGAt                                                                                                                                <  2:270524/149‑1 (MQ=255)
                     gTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATaa                                                                                                                              <  2:56079/149‑1 (MQ=255)
                       ccAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAAtt                                                                                                                            <  1:376598/149‑1 (MQ=255)
                         aGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGattttc                                                                                                                                       <  1:116614/136‑5 (MQ=255)
                         aGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGattttc                                                                                                                                       >  2:116614/1‑132 (MQ=255)
                          gTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGaa                                                                                                                         >  2:611151/1‑149 (MQ=255)
                              ttCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                     <  2:299306/149‑1 (MQ=255)
                              ttCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                     >  1:135291/1‑149 (MQ=255)
                               tCAGAGCCCCGTATTCAGCGAACACATTAACCTGTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTGAACCTGATATCAACCCGATAATTGAATCAtt                                                                                                                    >  1:231357/1‑149 (MQ=255)
                                 aGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTaa                                                                                                                  >  2:658087/1‑149 (MQ=255)
                                  gAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATCGAATCATTAAt                                                                                                                 <  3:633132/149‑1 (MQ=255)
                                      cccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAgg                                                                                                             <  3:543771/149‑1 (MQ=255)
                                      cccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAgg                                                                                                             <  1:209399/149‑1 (MQ=255)
                                       gcGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATGTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGc                                                                                                            >  4:518434/2‑149 (MQ=255)
                                       ccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCa                                                                                                                      <  1:39035/139‑1 (MQ=255)
                                       ccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCa                                                                                                                      >  2:39035/1‑139 (MQ=255)
                                       ccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGc                                                                                                            >  1:327407/1‑149 (MQ=255)
                                               aGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCa                                                                                                    <  4:409699/149‑1 (MQ=255)
                                                 cGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGt                                                                                                  <  1:464393/149‑1 (MQ=255)
                                                     cacaTTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGaaa                                                                                              <  1:191600/149‑1 (MQ=255)
                                                            aCCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCgg                                                                                       <  1:658087/149‑1 (MQ=255)
                                                            aCCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCgg                                                                                       <  4:262556/149‑1 (MQ=255)
                                                               tctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                     <  1:659673/116‑1 (MQ=255)
                                                               tctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                     >  2:659673/1‑116 (MQ=255)
                                                               tctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTa                                                                                    <  2:231357/149‑1 (MQ=255)
                                                               tctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTa                                                                                    >  3:180972/1‑149 (MQ=255)
                                                                  ttAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTc                                                                                         >  3:586958/1‑141 (MQ=255)
                                                                  ttAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTc                                                                                         <  4:586958/141‑1 (MQ=255)
                                                                    aaTTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCt                                                                               >  2:716273/1‑149 (MQ=255)
                                                                     aTTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAGTTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTc                                                                              >  3:742013/1‑149 (MQ=255)
                                                                     aTTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTc                                                                              >  1:709760/1‑149 (MQ=255)
                                                                                 tCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTg                                                                  <  1:611151/149‑1 (MQ=255)
                                                                                              aaGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCt                                                     >  3:222958/1‑149 (MQ=255)
                                                                                               aGGTGAATATGGTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTa                                                    <  3:518434/149‑1 (MQ=255)
                                                                                               aGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTa                                                    <  2:135291/149‑1 (MQ=255)
                                                                                                 gTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTaaaaa                                                         >  4:32229/1‑142 (MQ=255)
                                                                                                 gTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTaaaaa                                                         <  3:32229/142‑1 (MQ=255)
                                                                                                    aaTATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGtt                                               <  4:327068/149‑1 (MQ=255)
                                                                                                        tGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATcc                                           <  4:94355/149‑1 (MQ=255)
                                                                                                           tAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCa                                        >  3:53582/1‑149 (MQ=255)
                                                                                                                           gCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGc                        >  2:331394/1‑149 (MQ=255)
                                                                                                                                gAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTAc                                                   <  1:46497/117‑1 (MQ=255)
                                                                                                                                gAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTAc                                                   >  2:46497/1‑117 (MQ=255)
                                                                                                                                 aaaGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTCATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGCAAAATAATTGCGTAAAAAATCTACCCTTTTCCGGACTTTTGCTGACGTTCCCAGGCtt                 >  2:696460/1‑149 (MQ=255)
                                                                                                                                  aaGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATAc                                                                    >  2:197060/1‑98 (MQ=255)
                                                                                                                                  aaGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATAc                                                                    <  1:197060/98‑1 (MQ=255)
                                                                                                                                  aaGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGa                              >  3:128218/1‑136 (MQ=255)
                                                                                                                                  aaGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGa                              <  4:128218/136‑1 (MQ=255)
                                                                                                                                    ggAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGa               >  2:749159/1‑149 (MQ=255)
                                                                                                                                    ggAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGa               >  2:512988/1‑149 (MQ=255)
                                                                                                                                     gttgataTCAGGTTAAAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGaa              >  2:170690/7‑149 (MQ=255)
                                                                                                                                          aaTCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATa         >  2:22800/1‑149 (MQ=255)
                                                                                                                                          aaTCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTCGaaa             <  3:793735/145‑1 (MQ=255)
                                                                                                                                          aaTCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTCGaaa             >  4:793735/1‑145 (MQ=255)
                                                                                                                                           aTCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATaa        <  1:716273/149‑1 (MQ=255)
                                                                                                                                            tCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATaaa       <  4:222958/149‑1 (MQ=255)
                                                                                                                                               ggTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACgg    >  1:783420/1‑149 (MQ=255)
                                                                                                                                                gTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGc   >  1:468654/1‑149 (MQ=255)
                                                                                                                                                 ttttAACCTGATATCAACCCGATAGTTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGcc  <  4:742013/149‑1 (MQ=255)
                                                                                                                                                 ttttAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGcc  <  2:709760/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCC  >  NC_000913/1880704‑1880997

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: