Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,511,139 G→C 7.3% W379C (TGG→TGC ydcR → putative DNA‑binding transcriptional regulator and putative aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,511,1390GC 7.3% 220.7 / 6.5 67W379C (TGG→TGCydcRputative DNA‑binding transcriptional regulator and putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (36/26);  new base C (3/2);  total (39/28)
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.81e-01

TAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGG  >  NC_000913/1510992‑1511285
                                                                                                                                                   |                                                                                                                                                  
tAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTggc                                                                                                                                                   >  2:513658/1‑149 (MQ=255)
tAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTggc                                                                                                                                                   >  1:20589/1‑149 (MQ=255)
tAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTggc                                                                                                                                                   >  1:38979/1‑149 (MQ=255)
tAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTggc                                                                                                                                                   >  4:140220/1‑149 (MQ=255)
tAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTggc                                                                                                                                                   >  1:534237/1‑149 (MQ=255)
 aaaaTTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTggca                                                                                                                                                  <  3:679431/149‑1 (MQ=255)
 aaaaTTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTggca                                                                                                                                                  >  4:679431/1‑149 (MQ=255)
       cAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGCCAGGCAcg                                                                                                                                            <  2:646632/149‑2 (MQ=255)
         aCGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCActgc                                                                                                                                          <  4:584687/149‑1 (MQ=255)
           gCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCActgctg                                                                                                                                        <  1:207876/149‑1 (MQ=255)
           gCTTGCAGTTGATGAGTACGCTTGCCACCATCGCGCCGATACAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACAACGCCCATCTTCGTCGCATGCGTCGCCAGCTTGCGTAACGTAAACAACGTTCCTGCCAGGCAcgttgt                                                                                                                                        <  2:760884/149‑6 (MQ=255)
            cTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACtgctgc                                                                                                                                       <  4:236341/149‑1 (MQ=255)
             ttGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCg                                                                                                                                      >  3:697200/1‑149 (MQ=255)
                     gatgaGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATctgc                                                                                                                              >  1:627213/1‑149 (MQ=255)
                     gatgaGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATctgc                                                                                                                              >  2:735857/1‑149 (MQ=255)
                       tgaGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATctgcct                                                                                                                            >  2:87188/1‑149 (MQ=255)
                       tgaGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATctgcct                                                                                                                            <  4:697200/149‑1 (MQ=255)
                          gTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCg                                                                                                                         >  2:259149/1‑149 (MQ=255)
                          gTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCg                                                                                                                         >  4:538584/1‑149 (MQ=255)
                          gTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCg                                                                                                                         >  2:517133/1‑149 (MQ=255)
                              gCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAg                                                                                                                     <  3:140220/149‑1 (MQ=255)
                              gCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTCCTGCGTTATCTGCCTGCGGAAg                                                                                                                     <  1:469965/149‑1 (MQ=255)
                                    caccaGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGaaaa                                                                                                               <  1:107223/149‑1 (MQ=255)
                                    caccaGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGaaaa                                                                                                               <  2:38979/149‑1 (MQ=255)
                                     accaGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCTAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAAt                                                                                                              >  2:583227/1‑149 (MQ=255)
                                      ccaGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAAtt                                                                                                             >  4:289401/1‑149 (MQ=255)
                                       caGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATtc                                                                                                            <  4:184289/149‑1 (MQ=255)
                                        aGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATtca                                                                                                           <  2:639095/149‑1 (MQ=255)
                                        aGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATtca                                                                                                           >  4:298697/1‑149 (MQ=255)
                                         gCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATtcat                                                                                                          <  4:371708/149‑1 (MQ=255)
                                           tCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATtcatca                                                                                                        >  3:127441/1‑149 (MQ=255)
                                             aCCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATa                                                                                                      >  1:466760/1‑149 (MQ=255)
                                                      aCTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTg                                                                                             >  2:418087/1‑149 (MQ=255)
                                                                 tggATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTtctc                                                                                  <  2:513028/149‑1 (MQ=255)
                                                                     ttACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGt                                                                              >  2:752634/1‑149 (MQ=255)
                                                                       aCCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTg                                                                            <  2:402660/149‑1 (MQ=255)
                                                                        ccTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTgg                                                                           >  4:409119/1‑149 (MQ=255)
                                                                         cTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGa                                                                          <  2:627213/149‑1 (MQ=255)
                                                                              cACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTcc                                                                     <  3:241364/149‑1 (MQ=255)
                                                                               aCGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTccc                                                                    >  2:586708/1‑149 (MQ=255)
                                                                               aCGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTccc                                                                    >  2:56408/1‑149 (MQ=255)
                                                                                    ccGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGc                                                               >  3:610375/1‑149 (MQ=255)
                                                                                      gATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCg                                                             >  1:487690/1‑149 (MQ=255)
                                                                                           gacgCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTctct                                                                                 >  3:706317/1‑124 (MQ=255)
                                                                                           gacgCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTctct                                                                                 <  4:706317/124‑1 (MQ=255)
                                                                                              gCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGc                                                     >  4:116802/1‑149 (MQ=255)
                                                                                                ccATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTtctc                                                                                  >  2:13203/1‑118 (MQ=255)
                                                                                                ccATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTtctc                                                                                  <  1:13203/118‑1 (MQ=255)
                                                                                                 cATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCgg                                                  >  3:451203/1‑149 (MQ=255)
                                                                                                    cTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGa                                               <  1:566535/149‑1 (MQ=255)
                                                                                                      tcgtcgCCTGCGTCGCCATCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCAGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAAt                                             <  4:610375/149‑1 (MQ=255)
                                                                                                         tcgCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTaa                                          >  3:641758/1‑149 (MQ=255)
                                                                                                          cgCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAg                                         <  3:279428/149‑1 (MQ=255)
                                                                                                          cgCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAg                                         >  4:279428/1‑149 (MQ=255)
                                                                                                          cgCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAg                                         <  2:466760/149‑1 (MQ=255)
                                                                                                          cgCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAg                                         <  1:24074/149‑1 (MQ=255)
                                                                                                                     ccAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCCTCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCAc                              <  3:516218/149‑1 (MQ=255)
                                                                                                                     ccAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGt                                                                              <  3:691463/101‑1 (MQ=255)
                                                                                                                     ccAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGt                                                                              >  4:691463/1‑101 (MQ=255)
                                                                                                                       aGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTg                            >  4:11495/1‑149 (MQ=255)
                                                                                                                         cTTGCGGAACGTAAACAACGTGCCTGGCAGGCGCTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGAc                          >  2:754370/1‑149 (MQ=255)
                                                                                                                            gCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGca                       >  1:107778/1‑149 (MQ=255)
                                                                                                                              ggAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGcatc                     <  3:344808/149‑1 (MQ=255)
                                                                                                                               gAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGcatga                    >  3:617907/1‑147 (MQ=255)
                                                                                                                               gAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGcatca                    >  1:403090/1‑149 (MQ=255)
                                                                                                                                 aCGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCata                  >  3:559204/1‑149 (MQ=255)
                                                                                                                                   gTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATc                >  1:263260/1‑149 (MQ=255)
                                                                                                                                      gcagcagTGCCTGCCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGt             >  4:613096/7‑149 (MQ=255)
                                                                                                                                      aacaacGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGt             >  2:752675/1‑149 (MQ=255)
                                                                                                                                         gcagTGCCTGCCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTAtt          >  1:624777/4‑149 (MQ=255)
                                                                                                                                           agTGCCTGCCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTgc        >  2:503403/2‑149 (MQ=255)
                                                                                                                                              gCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGcc     <  4:451203/149‑1 (MQ=255)
                                                                                                                                              gCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGcc     >  3:537865/1‑149 (MQ=255)
                                                                                                                                                 tGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCggg  >  2:409768/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGG  >  NC_000913/1510992‑1511285

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