Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,486,223 A→T 7.1% K36* (AAA→TAA)  crp → cAMP‑activated global transcription factor, mediator of catabolite repression

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,486,2230AT 7.1% 205.7 / 5.5 66K36* (AAA→TAA) crpcAMP‑activated global transcription factor, mediator of catabolite repression
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base T (3/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGC  >  NC_000913/3486075‑3486371
                                                                                                                                                    |                                                                                                                                                     
ttATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGaaa                                                                                                                                                       >  1:716080/1‑149 (MQ=255)
    cTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTTTTCACCCTGGTGAATAAGc                                                                                                                                                   <  1:470788/149‑1 (MQ=255)
    cTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGc                                                                                                                                                   <  1:749838/149‑1 (MQ=255)
    cTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGc                                                                                                                                                   <  1:483381/149‑1 (MQ=255)
    cTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGc                                                                                                                                                   <  1:199683/149‑1 (MQ=255)
     tGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCg                                                                                                                                                  <  2:649013/149‑1 (MQ=255)
       gCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTTTTCACCCTAGTGAATAAGCgtg                                                                                                                                                <  2:448504/149‑3 (MQ=255)
       gCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGa                                                                                                                                                <  2:59658/149‑1 (MQ=255)
                aaaGCTTATAACCGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGta                                                                                                                                       <  1:106691/149‑1 (MQ=255)
                       ataACAGAGGATAATCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGGAACGCTGTACTACATc                                                                                                                                <  1:491011/149‑1 (MQ=255)
                         aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTCCCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGt                                                                                                                              <  1:222927/149‑1 (MQ=255)
                         aaCAGAGGATAACCGCGCATGGCGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTTCCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGt                                                                                                                              <  1:222794/149‑1 (MQ=255)
                          aCAGCGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGGCCCGACTCTCGGATGGTTCTTGTCTCATTGCCACAATCATAAGTGCCCATCCAAGAGCACGCTTATTCACCAGGGAGAAAAGGCGGAAACGCTGTACTACATCGta                                                                                                                             >  1:86961/1‑148 (MQ=255)
                             gagGATAACCGCGCATGGTGCTTGGCAAGCCGCAAACAGCCCCGACCATCGAATGGTTCTTGTCTCATTGCCACATTCATAGTTACGCATCCAAGAGCACGCTTATTCACCAGGGGGAAAAAGCGGAAACGCTGTACTACATCGTTaaa                                                                                                                          <  2:133129/149‑1 (MQ=255)
                                    aCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCGGTGCTACATCGTTAAAGGCTCtg                                                                                                                   >  2:549615/1‑149 (MQ=255)
                                     ccGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCtgt                                                                                                                  <  1:353834/149‑1 (MQ=255)
                                       gcgcATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTgg                                                                                                                <  2:396688/149‑1 (MQ=255)
                                             tgggCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGc                                                                                                          <  1:206247/146‑1 (MQ=255)
                                             ggTGCGTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGc                                                                                                          <  1:348426/149‑1 (MQ=255)
                                                  ttGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATc                                                                                                     >  1:250179/1‑149 (MQ=255)
                                                      cAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCGGATCAAAg                                                                                                >  2:603469/1‑149 (MQ=255)
                                                              aaaCAGACCCGACTGTCGAATGGTTATTGTCTCGTTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAGAGCGGAATCGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAagag                                                                                        <  1:477811/149‑1 (MQ=255)
                                                              aaaCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTCCATCGTAAAAGGCTCTGTGGCAGTGCTGACCAAAGCCGGagag                                                                                        >  2:56876/1‑149 (MQ=255)
                                                              aaaCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAagag                                                                                        <  2:14033/149‑1 (MQ=255)
                                                                         aCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAAt                                                                             <  1:117408/149‑1 (MQ=255)
                                                                         aCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAAt                                                                             <  1:760945/149‑1 (MQ=255)
                                                                         aCTCTCGAATGGTTCTTGTCTCATTGCAACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAAt                                                                             <  2:571795/149‑1 (MQ=255)
                                                                             tcGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATc                                                                         <  1:580406/149‑1 (MQ=255)
                                                                               gAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACAAAGAGGGTAAAGAAATGATccg                                                                       >  1:600486/1‑148 (MQ=255)
                                                                                aaTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCctc                                                                      >  1:298076/1‑149 (MQ=255)
                                                                                  tggaccTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAGCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCtctc                                                                    <  1:97884/144‑1 (MQ=255)
                                                                                          gTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGa                                                            >  1:633900/1‑149 (MQ=255)
                                                                                           tctcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGaa                                                           >  1:749990/1‑149 (MQ=255)
                                                                                           tctcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGaa                                                           >  1:630321/1‑149 (MQ=255)
                                                                                           tctcATTGCCACATTCATAAGTACCCAGCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGaa                                                           <  1:763247/149‑1 (MQ=255)
                                                                                            ctcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAAt                                                          >  1:122186/1‑149 (MQ=255)
                                                                                             tcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTTTTTATTTTGAAAAAGCGGAAACGCCGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                                         >  2:777541/1‑149 (MQ=255)
                                                                                             tcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                                         >  2:143257/1‑149 (MQ=255)
                                                                                             tcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                                         >  2:692569/1‑149 (MQ=255)
                                                                                                 tGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAggg                                                     >  1:739434/1‑149 (MQ=255)
                                                                                                   ccACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTg                                                   >  1:540347/1‑149 (MQ=255)
                                                                                                     acaTTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAt                                                 >  1:128150/1‑149 (MQ=255)
                                                                                                     acaTTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAt                                                 >  1:630801/1‑149 (MQ=255)
                                                                                                     acaTTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAt                                                 >  1:91470/1‑149 (MQ=255)
                                                                                                      cagtCATAAGTACCCACCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCACTGCCGTCGTGCTGTCAATACAATTAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAtt                                                <  1:186238/146‑1 (MQ=255)
                                                                                                         tCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTa                                             >  1:559759/1‑149 (MQ=255)
                                                                                                           ctAAGCACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCCGAATCAGGGTGATTTTAtt                                           <  2:534009/148‑1 (MQ=255)
                                                                                                           aTAAGTACCCATCCAAGAGCACGCTTCTTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAACGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTAtt                                           <  2:515730/149‑1 (MQ=255)
                                                                                                           aTAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTAtt                                           <  2:413604/149‑1 (MQ=255)
                                                                                                                   ccATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGGGAACt                                   >  1:387027/1‑149 (MQ=255)
                                                                                                                     aTCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTgg                                 >  1:84372/1‑149 (MQ=255)
                                                                                                                     aTCCAAGAGCACGCTTATTCACCAGGGTGAAAAAACGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTgg                                 >  2:436503/1‑149 (MQ=255)
                                                                                                                       ccAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGc                               >  1:444869/1‑149 (MQ=255)
                                                                                                                             gCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGttt                         <  1:551107/149‑1 (MQ=255)
                                                                                                                                cgcgTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGaa                      <  2:569712/149‑1 (MQ=255)
                                                                                                                                cacTTATACACCAGGGGGAAAAAGCGGAAACGCTGTACTACATCGTGCAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGTGTGTAAAGAAATGATCCTCTCCTATCTGAATCAGGTTGATTTTATTGGCGAACTGGGCCTGTTTGaa                      <  2:356022/147‑1 (MQ=255)
                                                                                                                                 gCTTTTTCACCCTGGTGAATAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCCGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAGTTTATTGGCGAACTGGGCCTGTTTGCag                     >  1:147789/1‑149 (MQ=255)
                                                                                                                                 gCTTTTTCACCCTGGTGAATAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTCGATCAGGGTGATGTTATTGGCGAACTGGGCCTGTTTGAag                     >  1:18893/1‑149 (MQ=255)
                                                                                                                                   ttATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTGCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGCagcg                   >  1:103980/1‑147 (MQ=255)
                                                                                                                                   ttATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAagag                   >  1:621897/1‑149 (MQ=255)
                                                                                                                                         aCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAg             >  1:74128/1‑149 (MQ=255)
                                                                                                                                         aCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAg             >  1:365787/1‑149 (MQ=255)
                                                                                                                                         aCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAg             >  1:197817/1‑149 (MQ=255)
                                                                                                                                           cctggTGAATAAGCGGAAACGCTGGACTGCATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGCTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGa           >  1:585304/4‑149 (MQ=255)
                                                                                                                                           cctggTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGa           >  2:595168/4‑149 (MQ=255)
                                                                                                                                              ggTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAGCTGGGCCTGTTTCAAGAGGGCCAGGCACg        >  2:399048/1‑149 (MQ=255)
                                                                                                                                                tGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTGCTTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTa      >  1:339243/1‑149 (MQ=255)
                                                                                                                                                  aaaaaGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGCAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGACTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTGgc    >  1:306508/1‑149 (MQ=255)
                                                                                                                                                   aaaaGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAggg   >  2:665020/1‑147 (MQ=255)
                                                                                                                                                    aaaGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAgcgc  <  1:429009/149‑1 (MQ=255)
                                                                                                                                                    aaaGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAgcgc  <  1:724758/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGC  >  NC_000913/3486075‑3486371

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