Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,176,799 T→C 6.1% E12G (GAG→GGG)  cpdA ← 3',5' cAMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,176,7990TC 6.1% 250.4 / 4.5 80E12G (GAG→GGG) cpdA3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base C (3/2);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

GACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGC  >  NC_000913/3176654‑3176946
                                                                                                                                                 |                                                                                                                                                   
gACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAc                                                                                                                                                  <  4:358807/149‑1 (MQ=255)
gACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAc                                                                                                                                                  <  2:135863/149‑1 (MQ=255)
gACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAc                                                                                                                                                  >  1:440604/1‑149 (MQ=255)
     tCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAccagcc                                                                                                                                             <  2:540345/149‑1 (MQ=255)
       aGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAccagccag                                                                                                                                           >  3:494693/1‑149 (MQ=255)
       aGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAccagccag                                                                                                                                           >  2:383635/1‑149 (MQ=255)
         gTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCagag                                                                                                                                         >  1:120530/1‑149 (MQ=255)
         gTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTCCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCagag                                                                                                                                         <  1:95763/149‑1 (MQ=255)
          tCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAgg                                                                                                                                        >  1:352418/1‑149 (MQ=255)
                    gCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAAc                                                                                                                              <  3:397188/149‑1 (MQ=255)
                       ggTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGCCTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAgg                                                                                                                           <  4:273268/149‑1 (MQ=255)
                          gtgGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCtt                                                                                                                        <  4:566133/149‑1 (MQ=255)
                           tgGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCttt                                                                                                                       <  2:333024/149‑1 (MQ=255)
                              ccGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCa                                                                                                                    >  1:170867/1‑149 (MQ=255)
                                gAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTc                                                                                                                      <  2:74121/145‑1 (MQ=255)
                                gAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTc                                                                                                                      >  1:74121/1‑145 (MQ=255)
                                gAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCaaa                                                                                                                  <  3:193416/149‑1 (MQ=255)
                                 aaTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAAt                                                                                                                 <  2:415596/149‑1 (MQ=255)
                                  aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCCCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATg                                                                                                                >  1:46601/1‑149 (MQ=255)
                                  aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCCCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATg                                                                                                                >  1:49782/1‑149 (MQ=255)
                                   tCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATgg                                                                                                               <  1:176006/149‑1 (MQ=255)
                                    cGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGgt                                                                                                              <  4:494693/149‑1 (MQ=255)
                                    cGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGgt                                                                                                              <  2:300224/149‑1 (MQ=255)
                                     gCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGgtg                                                                                                             <  1:532365/149‑1 (MQ=255)
                                       cTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTc                                                                                                           <  3:94652/149‑1 (MQ=255)
                                       cTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTc                                                                                                           <  3:127123/149‑1 (MQ=255)
                                          cAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCCCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCtt                                                                                                        <  2:49782/149‑1 (MQ=255)
                                          cAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCCCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCtt                                                                                                        <  2:46601/149‑1 (MQ=255)
                                              aCCGCCTCGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGt                                                                                                    <  1:216394/149‑1 (MQ=255)
                                               ccGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGtt                                                                                                   <  2:120530/149‑1 (MQ=255)
                                                   cTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCac                                                                                               <  4:524336/149‑1 (MQ=255)
                                                     ggTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacga                                                                                             >  4:242549/1‑149 (MQ=255)
                                                      gTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacgac                                                                                            >  3:438770/1‑149 (MQ=255)
                                                      gTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacgac                                                                                            <  3:219619/149‑1 (MQ=255)
                                                           tctcCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCtaa                                                                                       >  1:369120/1‑149 (MQ=255)
                                                           tctcCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCtaa                                                                                       >  4:544867/1‑149 (MQ=255)
                                                            ctcCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCtaat                                                                                      >  3:282817/1‑149 (MQ=255)
                                                            ctcCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCtaat                                                                                      >  3:189430/1‑149 (MQ=255)
                                                               ccAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGT‑AAATCCTTGCTCTTGCCCTAACCGCCAGGGGAAAGGGTAACAGGCCTTCCAAATTGTGGCCTTAAGTTTTCCGACGCCAATaaaa                                                                                  >  1:174153/1‑148 (MQ=255)
                                                                 agttgtTTACCCCTACCAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGGCTCTGGCCTCACCAGCCAGAGGAAGTGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAAcc                                                                                 <  2:362934/146‑1 (MQ=255)
                                                                  ggTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCg                                                                                <  2:368268/149‑1 (MQ=255)
                                                                      tttACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAAt                                                                            >  1:398464/1‑149 (MQ=255)
                                                                           cccTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGc                                                                         >  1:325465/1‑147 (MQ=255)
                                                                           cccTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGc                                                                         <  2:325465/147‑1 (MQ=255)
                                                                           cccTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCa                                                                       >  3:529552/1‑149 (MQ=255)
                                                                           cccTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCa                                                                       <  4:153181/149‑1 (MQ=255)
                                                                           cccTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCa                                                                       >  4:525517/1‑149 (MQ=255)
                                                                             cTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATc                                                                     >  2:88014/1‑149 (MQ=255)
                                                                             cTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATc                                                                     <  4:8933/149‑1 (MQ=255)
                                                                               aaCAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGc                                                                   >  1:299021/1‑149 (MQ=255)
                                                                                  aGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTcc                                                                >  2:265395/1‑149 (MQ=255)
                                                                                    ggCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCAt                                                              <  2:170867/149‑1 (MQ=255)
                                                                                      cTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCAtgt                                                            <  3:544867/149‑1 (MQ=255)
                                                                                        tCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCATCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGc                                                          <  3:392829/149‑1 (MQ=255)
                                                                                              ttttGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCaaa                                                                                                                  <  3:261057/87‑1 (MQ=255)
                                                                                              ttttGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCaaa                                                                                                                  >  4:261057/1‑87 (MQ=255)
                                                                                              ttttGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                    >  2:55210/1‑149 (MQ=255)
                                                                                              ttttGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                    <  2:139735/149‑1 (MQ=255)
                                                                                                   tgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATc                                               <  4:334387/149‑1 (MQ=255)
                                                                                                      aaaCAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCa                                            >  3:482645/1‑149 (MQ=255)
                                                                                                      aaaCAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGAGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCa                                            >  4:479616/1‑149 (MQ=255)
                                                                                                        aCAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAAc                                          <  2:172929/149‑1 (MQ=255)
                                                                                                           ggTGAGTGTCGGTAATTAGTAAAATCCTGACTCTGGCCTCTCCAGACAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAg                                       <  2:225762/149‑1 (MQ=255)
                                                                                                             tGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTc                                     >  4:209260/1‑149 (MQ=255)
                                                                                                                    cGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATc                                               <  1:382425/132‑1 (MQ=255)
                                                                                                                    cGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATc                                               >  2:382425/1‑132 (MQ=255)
                                                                                                                      gTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTaaaaa                            >  3:561123/1‑149 (MQ=255)
                                                                                                                       tAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAAc                           <  4:316441/149‑1 (MQ=255)
                                                                                                                       tAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAAc                           >  1:242584/1‑149 (MQ=255)
                                                                                                                         aTTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTg                         >  4:151269/1‑149 (MQ=255)
                                                                                                                         aTTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTg                         >  4:10520/1‑149 (MQ=255)
                                                                                                                         aTTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTg                         >  1:499867/1‑149 (MQ=255)
                                                                                                                               aaaaTCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAAt                   >  4:362411/1‑149 (MQ=255)
                                                                                                                                 aaTCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAAttt                 <  3:36388/149‑1 (MQ=255)
                                                                                                                                   tCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGa               <  1:88014/149‑1 (MQ=255)
                                                                                                                                    ccTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCGTCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGAt              <  2:373033/149‑1 (MQ=255)
                                                                                                                                      tGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGc            >  3:269726/1‑149 (MQ=255)
                                                                                                                                           ctGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCATTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtc       <  3:479616/149‑1 (MQ=255)
                                                                                                                                               ccTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtcgtcg   >  1:46415/1‑149 (MQ=255)
                                                                                                                                                cTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTcgc  >  1:206746/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGC  >  NC_000913/3176654‑3176946

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