Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,508,007 T→A 5.2% R409R (CGT→CGA ydcP → putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,508,0070TA 5.2% 261.0 / 3.6 88R409R (CGT→CGAydcPputative peptidase
Reads supporting (aligned to +/- strand):  ref base T (39/43);  new base A (2/3);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAA  >  NC_000913/1507860‑1508153
                                                                                                                                                   |                                                                                                                                                  
tACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTc                                                                                                                                                   >  7:373084/1‑149 (MQ=255)
tACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCGCCAGCGAACGTc                                                                                                                                                   >  3:272642/1‑149 (MQ=255)
  cGGTCGAGAAAAGCGGAGTAAATCAGTACCGCGTCTGGGCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAAGCGTAATCTTGATCATAACTGGGAGCAGGCACTGACAAAAACCTCCAGCGAACGTCgg                                                                                                                                                 <  7:11601/149‑1 (MQ=255)
    gTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGt                                                                                                                                               >  4:188695/1‑149 (MQ=255)
     tCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTg                                                                                                                                              <  2:273801/149‑1 (MQ=255)
       gagaAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGCAATGCCAGCAGATTTGCACAAAATTCGTCCACATCCCCCCCTAACCCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAACACCCTCCAGAGCACGGCGGGTGGc                                                                                                                                            >  3:276942/1‑149 (MQ=255)
       gagaAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTgg                                                                                                                                             <  7:20096/148‑1 (MQ=255)
        agaAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCCCCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCg                                                                                                                                           >  8:119175/1‑149 (MQ=255)
               cGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACa                                                                                                                                    <  1:5316/149‑1 (MQ=255)
               cGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACa                                                                                                                                    <  8:316681/149‑1 (MQ=255)
                 gagaAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACAtt                                                                                                                                  >  4:252955/1‑149 (MQ=255)
                    aaaaTCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGaa                                                                                                                               <  1:241971/149‑1 (MQ=255)
                    aaaaTCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGACGGGTGGCGGTAGACATTGaa                                                                                                                               <  2:358035/149‑1 (MQ=255)
                     aaaTCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAAc                                                                                                                              <  7:85430/149‑1 (MQ=255)
                     aaaTCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCAAAAAATTCGTCCACATCACCCACTATACCGTAATCTTGATCATAACTTGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCTGTAGACATTGAAc                                                                                                                              <  5:120816/149‑1 (MQ=255)
                              ccGCGTCTGGCCCAATGAAATGCCAGCATATTTGCACAACATTCGTCCACATCACCCACTAAACAGTAATCTTTATCATAACTGGCAGCAGGCACTTACAAAAACCTCCAGAGAAAGTCGGGGGGCGGTATACATGTAACTGGGCGGCt                                                                                                                     >  6:354453/1‑149 (MQ=255)
                              ccGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCt                                                                                                                     <  8:267231/149‑1 (MQ=255)
                              ccGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCt                                                                                                                     >  3:212546/1‑149 (MQ=255)
                              ccGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCt                                                                                                                     >  4:312926/1‑149 (MQ=255)
                               cgcgTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTg                                                                                                                    <  7:313402/149‑1 (MQ=255)
                               cgcgTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGACGGGTGGCGGTAGACATTGAACTGGGCGGCTg                                                                                                                    >  2:58693/1‑149 (MQ=255)
                               cgcgTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGACGGGTGGCGGTAGACATTGAACTGGGCGGCTg                                                                                                                    >  2:64467/1‑149 (MQ=255)
                                gcgTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTgg                                                                                                                   <  3:188695/149‑1 (MQ=255)
                                 cgTCTGGCCCAATGACATGCCAGCAGATTTGCACAAAATGCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGc                                                                                                                  <  5:354453/149‑1 (MQ=255)
                                  gTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCa                                                                                                                 <  8:45171/149‑1 (MQ=255)
                                  gTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTAGGGTGGCGGTAGACATTGAACTGGGCGGCTgg                                                                                                                   >  1:161305/1‑147 (MQ=255)
                                       gCCCAATGAAATGCCAGCAGATTTGCACAAACTTCGTCCACCTCACCCACTACACCATAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGa                                                                                                              <  2:351989/147‑1 (MQ=255)
                                       gCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGa                                                                                                              <  8:26247/147‑1 (MQ=255)
                                         ccAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGaacaa                                                                                                          >  3:330652/1‑149 (MQ=255)
                                         ccAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGaacaa                                                                                                          <  1:77602/149‑1 (MQ=255)
                                          cAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATACCTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGaacaac                                                                                                         >  8:302903/1‑149 (MQ=255)
                                          cAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGaacaac                                                                                                         <  1:410422/149‑1 (MQ=255)
                                             tGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGa                                                                                                      <  3:252955/149‑1 (MQ=255)
                                             tGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGa                                                                                                      <  8:373084/149‑1 (MQ=255)
                                               aaaTGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATAATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATAGAACTCGGCGGCTGGCAGGAACAACAGAt                                                                                                     >  4:113959/1‑148 (MQ=255)
                                                aaTGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTc                                                                                                   >  4:88807/1‑149 (MQ=255)
                                                aaTGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTAGCAGGAACAACTGATTc                                                                                                   <  2:299280/149‑1 (MQ=255)
                                                  tGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGACCGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGCTTCt                                                                                                  >  6:181517/1‑148 (MQ=255)
                                                    ccAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGAc                                                                                               <  4:212546/149‑1 (MQ=255)
                                                    ccAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGAc                                                                                               <  8:53765/149‑1 (MQ=255)
                                                     cagcagATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGAcc                                                                                              >  4:162237/1‑149 (MQ=255)
                                                      agcagATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGAccc                                                                                             >  7:114004/1‑149 (MQ=255)
                                                          gATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCAc                                                                                         <  3:187580/149‑1 (MQ=255)
                                                           aTTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCCAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCAc                                                                                         >  8:216076/1‑148 (MQ=255)
                                                                cacaAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTg                                                                                    >  5:226282/1‑148 (MQ=255)
                                                                   aaaaTTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAaga                                                                                <  8:405971/149‑1 (MQ=255)
                                                                          gTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGgtgt                                                                          >  2:13517/1‑148 (MQ=255)
                                                                          gTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGACGTGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTc                                                                         <  1:64467/149‑1 (MQ=255)
                                                                          gTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGACGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTc                                                                         <  1:58693/149‑1 (MQ=255)
                                                                           tCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGgtgt                                                                          >  5:400192/1‑147 (MQ=255)
                                                                             cacaTCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAg                                                                       >  7:94846/1‑148 (MQ=255)
                                                                                    cccACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGc                                                               <  6:400192/149‑1 (MQ=255)
                                                                                           aaCCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCt                                                        <  6:226282/149‑1 (MQ=255)
                                                                                            aCCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTg                                                       >  1:382883/1‑149 (MQ=255)
                                                                                             ccGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGaa                                                                                  <  5:100988/121‑1 (MQ=255)
                                                                                             ccGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGaa                                                                                  >  6:100988/1‑121 (MQ=255)
                                                                                             ccGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAATAGGGTGTCAGCATCACGCATACGCTg                                                       >  1:263140/1‑148 (MQ=255)
                                                                                              cGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGa                                                     >  7:178532/1‑149 (MQ=255)
                                                                                                 aaTCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACgg                                                  >  1:277851/1‑149 (MQ=255)
                                                                                                     ttGATCATAACTGGCAGCAGGCACTGCCAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGCACTGGGCGGCGGGCAGGAACAACTGATTATGACCCTCACAAGTGAAGAGGGTGTCAGCATCACGCATACTCTGGACGTGCAg                                              >  6:329927/1‑149 (MQ=255)
                                                                                                        aTCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTTGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTc                                           <  4:300321/149‑1 (MQ=255)
                                                                                                         tCATAACTGGCAGCGGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTc                                           >  3:197033/1‑148 (MQ=255)
                                                                                                                  gcagcaGGCACTGACAAAAACCTCCAGCGCACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGcc                                  >  1:26735/1‑148 (MQ=255)
                                                                                                                  gcagcaGGCACTGACAAAAACCTCCAGCGACCGTCGGTGGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGc                                   <  3:113959/147‑1 (MQ=255)
                                                                                                                   cagcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCaa                                <  3:29253/149‑1 (MQ=255)
                                                                                                                    agcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGGGTCAGCATCACGCATACGCTGGACGTGCAGTTCGACGAAGCCaa                                >  3:42652/1‑148 (MQ=255)
                                                                                                                     gcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAAt                               >  8:121044/1‑148 (MQ=255)
                                                                                                                     gcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATa                              <  2:382883/149‑1 (MQ=255)
                                                                                                                      cagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATaa                             <  3:6064/149‑1 (MQ=255)
                                                                                                                       agGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAAc                            <  3:312926/149‑1 (MQ=255)
                                                                                                                          cACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGcc                         >  1:250341/1‑149 (MQ=255)
                                                                                                                          cACTGACAAAAACCTCCAGCGAACGTCGGGGGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGcc                         <  4:272642/149‑1 (MQ=255)
                                                                                                                            atGACAAAAACCTCCAGCGAACGTAGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGAGCAGTACGACGAAGCCAAGAACGCCg                        >  8:304492/2‑148 (MQ=255)
                                                                                                                             tGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGCCGAAGCCAATAACGCCGa                       >  4:356110/1‑148 (MQ=255)
                                                                                                                               aCAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGaaa                     >  8:259786/1‑148 (MQ=255)
                                                                                                                                 aaaaaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAg                  >  5:310246/1‑149 (MQ=255)
                                                                                                                                 aaaaaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAg                  <  3:88807/149‑1 (MQ=255)
                                                                                                                                 aaaaaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCAGACGCTGGACGGGCAGTTCGACGAAGCCAAATACGCCGAAAAAg                  >  1:61309/1‑149 (MQ=255)
                                                                                                                                    aaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCaa               <  4:330652/149‑1 (MQ=255)
                                                                                                                                     aCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAAt              <  7:302903/149‑1 (MQ=255)
                                                                                                                                     aCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAAt              <  8:175437/149‑1 (MQ=255)
                                                                                                                                        tCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCa                                                              <  4:100465/98‑1 (MQ=255)
                                                                                                                                        tCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCa                                                              >  3:100465/1‑98 (MQ=255)
                                                                                                                                        tCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGaa           >  2:129222/1‑149 (MQ=255)
                                                                                                                                             cGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATc      <  2:250341/149‑1 (MQ=255)
                                                                                                                                                 cGTCGGGTGGCGGTAGACATTTAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGaa  <  8:282607/149‑1 (MQ=255)
                                                                                                                                                 cGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGaa  <  7:216076/149‑1 (MQ=255)
                                                                                                                                                 cGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGaa  <  8:94846/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAA  >  NC_000913/1507860‑1508153

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