Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,604,754 T→A 7.8% S199T (TCT→ACT)  bglJ → bgl operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,7540TA 7.8% 302.0 / 14.1 90S199T (TCT→ACT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (42/41);  new base A (5/2);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01

AGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGCCACATTG  >  NC_000913/4604607‑4604901
                                                                                                                                                    |                                                                                                                                                   
aGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTc                                                                                                                                                    <  2:257458/149‑1 (MQ=255)
   cGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                                                                                 <  4:454724/149‑1 (MQ=255)
     gACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACg                                                                                                                                               >  1:642456/1‑149 (MQ=255)
        gTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAg                                                                                                                                            >  1:453003/1‑149 (MQ=255)
           aaGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGAACTGACGcccagt                                                                                                                                         <  4:685995/149‑6 (MQ=255)
               ccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                                                                     <  2:169549/149‑1 (MQ=255)
                cGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                                                                                                                                    >  1:124549/1‑149 (MQ=255)
                  aCGGAGAGCGAAATGTTGCGCTTTATGTCGCCTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACCCTATATCAAAACGATTCGTGCTCATAAATTTAATCTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                                                                  <  3:468791/149‑1 (MQ=255)
                  aCGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                                                                  <  2:700362/149‑1 (MQ=255)
                  aCGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                                                                  <  4:165328/149‑1 (MQ=255)
                   cGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAg                                                                                                                                 <  1:307215/149‑1 (MQ=255)
                   cGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAg                                                                                                                                 <  4:566068/149‑1 (MQ=255)
                    ggAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                                                                                                                                >  2:8667/1‑149 (MQ=255)
                    ggAGAGAAAAATATTGCGGTTTAAGTCGCCGTGGTTATTAATGGCACAAAGTGCCGAGCAGGTTTAACGCAAAATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                                                                                                                                <  2:637303/149‑1 (MQ=255)
                         gagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCa                                                                                                                           >  1:193429/1‑149 (MQ=255)
                          agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAg                                                                                                                          >  2:64920/1‑149 (MQ=255)
                          agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAg                                                                                                                          >  2:167690/1‑149 (MQ=255)
                             aaTATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGata                                                                                                                       <  2:392196/149‑1 (MQ=255)
                               tatTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATAtt                                                                                                                     <  4:749972/149‑1 (MQ=255)
                                 ttGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCt                                                                                                                   >  4:230023/1‑149 (MQ=255)
                                 ttGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCt                                                                                                                   <  4:549989/149‑1 (MQ=255)
                                   gcgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                                                                                 <  2:660890/149‑1 (MQ=255)
                                   gcgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                                                                                 >  3:275832/1‑149 (MQ=255)
                                    cgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGtt                                                                                                                <  2:453003/149‑1 (MQ=255)
                                    cgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCGTAAACGCAATAACAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGtt                                                                                                                <  1:32739/149‑1 (MQ=255)
                                     gcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTta                                                                                                               >  2:346863/1‑149 (MQ=255)
                                      cTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat                                                                                                              <  4:365534/149‑1 (MQ=255)
                                                      cTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCg                                                                                               >  4:494988/1‑149 (MQ=255)
                                                       tACTCAATGACACAAATTGCCGAGCGACTAAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGAATTGCgt                                                                                              >  2:276712/1‑148 (MQ=255)
                                                       tACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGa                                                                                              >  2:387327/1‑149 (MQ=255)
                                                       tACTCAATGACACAAATTGCCGAGCAACTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGa                                                                                              >  4:118602/1‑149 (MQ=255)
                                                         cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                                                                  >  1:302910/1‑111 (MQ=255)
                                                         cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                                                                  <  2:302910/111‑1 (MQ=255)
                                                         cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaaa                                                                                            <  4:542691/149‑1 (MQ=255)
                                                            aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                                                                         >  4:570512/1‑149 (MQ=255)
                                                            aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                                                                         >  4:68400/1‑149 (MQ=255)
                                                            aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                                                                         >  4:117569/1‑149 (MQ=255)
                                                             aTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAg                                                                                        >  1:188532/1‑149 (MQ=255)
                                                              tGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGt                                                                                       <  2:142290/149‑1 (MQ=255)
                                                               gACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTa                                                                                      <  3:275467/149‑1 (MQ=255)
                                                                acacAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCGGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTaa                                                                                     <  1:631423/149‑1 (MQ=255)
                                                                  acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                                                                   >  3:390584/1‑149 (MQ=255)
                                                                  acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                                                                   >  3:498391/1‑149 (MQ=255)
                                                                     aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATgtgt                                                                                >  1:131081/1‑149 (MQ=255)
                                                                      aTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                                                                               >  1:275981/1‑149 (MQ=255)
                                                                      aTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                                                                               <  1:8667/149‑1 (MQ=255)
                                                                          ccGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCAt                                                                           <  2:131081/149‑1 (MQ=255)
                                                                           cGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATc                                                                          <  2:594879/149‑1 (MQ=255)
                                                                             agcagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATccc                                                                        <  3:102942/149‑1 (MQ=255)
                                                                             agcagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATccc                                                                        >  4:594367/1‑149 (MQ=255)
                                                                                agcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTAt                                                                     <  3:117569/149‑1 (MQ=255)
                                                                                     aaaCGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATc                                                                <  4:498391/149‑1 (MQ=255)
                                                                                        cGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGc                                                             <  3:394487/149‑1 (MQ=255)
                                                                                          cAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCat                                                           <  4:744240/149‑1 (MQ=255)
                                                                                           aaTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatg                                                          <  3:570512/149‑1 (MQ=255)
                                                                                               tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgt                                                                                  <  2:492974/121‑1 (MQ=255)
                                                                                               tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgt                                                                                  >  1:492974/1‑121 (MQ=255)
                                                                                                 aaaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg                                                    >  2:718082/1‑149 (MQ=255)
                                                                                                 aaaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg                                                    >  4:14617/1‑149 (MQ=255)
                                                                                                 aaaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg                                                    >  4:389946/1‑149 (MQ=255)
                                                                                                 aaaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg                                                    >  1:142751/1‑149 (MQ=255)
                                                                                                 aaaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg                                                    >  3:281567/1‑149 (MQ=255)
                                                                                                  aaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                                                                     >  1:564232/1‑67 (MQ=255)
                                                                                                  aaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                                                                     <  2:564232/67‑1 (MQ=255)
                                                                                                  aaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACGAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg                                                   >  4:315096/1‑149 (MQ=255)
                                                                                                  aaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg                                                   <  2:124549/149‑1 (MQ=255)
                                                                                                  aaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg                                                   >  4:220989/1‑149 (MQ=255)
                                                                                                     cGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGCTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGtt                                                <  1:226841/149‑1 (MQ=255)
                                                                                                           gTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGAACTGACGc                                                                                                                                              >  3:61454/1‑49 (MQ=255)
                                                                                                           gTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGAACTGACGc                                                                                                                                              <  4:61454/49‑1 (MQ=255)
                                                                                                             gCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTccc                                        <  3:315096/149‑1 (MQ=255)
                                                                                                               acaTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCg                                      >  1:378123/1‑149 (MQ=255)
                                                                                                                 aTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGc                                    >  4:345946/1‑149 (MQ=255)
                                                                                                                    aaTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGt                                 <  4:404432/149‑1 (MQ=255)
                                                                                                                     aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTg                                >  2:288679/1‑149 (MQ=255)
                                                                                                                      tttAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGATTACTTCCCCGGCAGTGc                               >  2:29525/1‑149 (MQ=255)
                                                                                                                      tttAATGGGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGc                               >  2:411982/1‑149 (MQ=255)
                                                                                                                         aaTGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCttt                            <  4:361460/149‑1 (MQ=255)
                                                                                                                             tgATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCtttcatt                        <  1:64920/149‑1 (MQ=255)
                                                                                                                                tGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCtttcatttca                     <  1:387327/149‑1 (MQ=255)
                                                                                                                                  tCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGc                   >  2:213432/1‑149 (MQ=255)
                                                                                                                                       accctgCGTCAGAACTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACa              >  1:617319/6‑149 (MQ=255)
                                                                                                                                         tGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGc                   <  4:144777/142‑1 (MQ=255)
                                                                                                                                         tGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGc                   >  3:144777/1‑142 (MQ=255)
                                                                                                                                         cctgCGTCAGAACTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAAt            >  4:335568/4‑149 (MQ=255)
                                                                                                                                         cctgCGTCAGAACTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAAt            >  1:44373/4‑149 (MQ=255)
                                                                                                                                          ctgCGTCAGAACTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATc           >  3:190979/3‑149 (MQ=255)
                                                                                                                                           ggCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCg          >  3:713027/1‑149 (MQ=255)
                                                                                                                                            gCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGc         <  1:758484/149‑1 (MQ=255)
                                                                                                                                               tCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGCcac      >  3:698000/1‑149 (MQ=255)
                                                                                                                                                   ttCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGCCACATtg  >  3:236686/1‑149 (MQ=255)
                                                                                                                                                   ttCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGCCACATtg  >  4:386627/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGCCACATTG  >  NC_000913/4604607‑4604901

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