Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,176,319 C→A 54.1% C172F (TGT→TTT)  cpdA ← 3',5' cAMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,176,3190CA54.1% 8.7 / 142.4 100C172F (TGT→TTT) cpdA3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base A (23/31);  total (48/52)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 102 total aligned reads displayed.

GGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAT  >  NC_000913/3176172‑3176464
                                                                                                                                                   |                                                                                                                                                 
ggCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAc                                                                                                                                                   >  3:395094/1‑148 (MQ=255)
 gCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTaaaa                                                                                                                                                 <  3:320796/149‑1 (MQ=255)
   gTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAACCTTCGCCAGCACGGGATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCACCTAAAAcc                                                                                                                                               >  5:308817/1‑149 (MQ=255)
       ccAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTACTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCa                                                                                                                                           <  8:485889/149‑1 (MQ=255)
        cagcagGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAg                                                                                                                                          >  4:291326/1‑149 (MQ=255)
        cagcagGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAg                                                                                                                                          <  4:395094/149‑1 (MQ=255)
        cagcagGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAg                                                                                                                                          <  2:32231/149‑1 (MQ=255)
        cagcagGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAg                                                                                                                                          <  5:463889/149‑1 (MQ=255)
         agcagGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAgg                                                                                                                                         >  6:315902/1‑149 (MQ=255)
            agGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAg                                                                                                                                      >  8:146341/1‑149 (MQ=255)
              gcggcgACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCg                                                                                                                                    <  2:323363/149‑1 (MQ=255)
                ggcgACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGc                                                                                                                                     <  4:437542/146‑1 (MQ=255)
                ggcgACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGc                                                                                                                                     >  3:437542/1‑146 (MQ=255)
                  cgACAATTCCAGTCGAGATCCAGCTACTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCAAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACACGCAGGTAGCTgat                                                                                                                                 >  8:418527/1‑148 (MQ=255)
                       aTTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGgatgatgat                                                                                                                           <  7:42624/149‑1 (MQ=255)
                        ttCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGatgatgatg                                                                                                                          <  6:192867/149‑1 (MQ=255)
                        ttCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGatgatgatg                                                                                                                          <  1:135281/149‑1 (MQ=255)
                            aGTCGAGATCCATCTCCTATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGc                                                                                                                                     <  1:136267/133‑1 (MQ=255)
                            aGTCGAGATCCAGCTCCTATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCTTTACGCAGAATGTGTTGATCTAGCAAACTACAACCCGCAGGTAGc                                                                                                                                     >  2:136267/1‑133 (MQ=255)
                               cGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAgg                                                                                                                                         <  7:199085/127‑1 (MQ=255)
                               cGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAgg                                                                                                                                         >  8:199085/1‑127 (MQ=255)
                                gagaTCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcag                                                                                                                    >  3:188701/1‑147 (MQ=255)
                                gagaTCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATgcagcag                                                                                                                    >  6:376721/1‑147 (MQ=255)
                                  gaTCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcag                                                                                                                 >  3:298366/1‑148 (MQ=255)
                                  gaTCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcag                                                                                                                 >  8:273965/1‑148 (MQ=255)
                                  gaTCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATgcagcagcag                                                                                                                 >  2:243669/1‑148 (MQ=255)
                                   aTCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACAACAACACGCAGGTATCGGATGATGATgcagcagcagc                                                                                                                >  4:159688/1‑148 (MQ=255)
                                      cAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCGCGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACg                                                                                                            >  2:78110/1‑149 (MQ=255)
                                       aGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTTATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACg                                                                                                            >  8:448657/1‑148 (MQ=255)
                                           ccTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCAAACTACAACCCGCAGAAAGCGGATGATGATGCAGCAGCAGCAACGTATgg                                                                                                       >  4:429593/1‑149 (MQ=255)
                                            cTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCGCGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATgg                                                                                                       <  1:78110/148‑1 (MQ=255)
                                            cTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATgg                                                                                                       <  7:171198/148‑1 (MQ=255)
                                            cTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATg                                                                                                        <  1:243669/147‑1 (MQ=255)
                                            cTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGc                                                                                                      <  3:334052/149‑1 (MQ=255)
                                             tgatgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCg                                                                                                     >  7:3353/1‑149 (MQ=255)
                                              gatgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcag                                                                                                                    >  5:169610/1‑133 (MQ=255)
                                              gatgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcag                                                                                                                    <  6:169610/133‑1 (MQ=255)
                                               atgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGtt                                                                                                   >  7:129292/1‑149 (MQ=255)
                                               atgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGtt                                                                                                   >  1:300107/1‑149 (MQ=255)
                                                tgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCg                                                                                                     >  7:166071/1‑146 (MQ=255)
                                                tgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCg                                                                                                     <  8:166071/146‑1 (MQ=255)
                                                 gaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCt                                                                                                 >  1:120859/1‑149 (MQ=255)
                                                   atatGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTgg                                                                                               <  8:129292/149‑1 (MQ=255)
                                                       gACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCa                                                                                           >  8:379772/1‑149 (MQ=255)
                                                        aCCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCAt                                                                                          <  4:188701/149‑1 (MQ=255)
                                                        aCCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCAt                                                                                          >  1:367355/1‑149 (MQ=255)
                                                         ccGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATc                                                                                         <  1:380053/149‑1 (MQ=255)
                                                         ccGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATc                                                                                         >  7:498864/1‑149 (MQ=255)
                                                           gCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAAc                                                                                                             >  3:99636/1‑127 (MQ=255)
                                                           gCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCgg                                                                                       <  8:100741/149‑1 (MQ=255)
                                                            cacaGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGc                                                                                      <  3:291326/149‑1 (MQ=255)
                                                             acaGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGcc                                                                                     <  8:3353/149‑1 (MQ=255)
                                                              caGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCa                                                                                    <  2:29463/149‑1 (MQ=255)
                                                                 cAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGgc                                                                                              >  5:276415/1‑136 (MQ=255)
                                                                 cAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGgc                                                                                              <  6:276415/136‑1 (MQ=255)
                                                                 cAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGtt                                                                                 <  7:146341/149‑1 (MQ=255)
                                                                  aaGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCg                                                                                        >  2:398811/1‑141 (MQ=255)
                                                                  aaGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCg                                                                                        <  1:398811/141‑1 (MQ=255)
                                                                   aGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTa                                                                               <  5:376721/149‑1 (MQ=255)
                                                                     tATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATAGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATTCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCGAGTTTACg                                                                             <  2:84413/149‑1 (MQ=255)
                                                                           aCGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCa                                                                        <  7:273965/148‑1 (MQ=255)
                                                                           aCGCGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACGGTGTTGAGCGAGCCACCTACACCCCGCAGGTAGCGGCTGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCa                                                                        <  6:308817/148‑1 (MQ=255)
                                                                              tGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCa                                                                    <  5:315902/149‑1 (MQ=255)
                                                                              tGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCa                                                                    <  2:367355/149‑1 (MQ=255)
                                                                                 ggAAACTTCGCTAGCACTGTATCCAGTTCGCCCGCTTTACGCAGACTTTGTTGATAGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATTGCGTTCTGGCGCATCGTCCAGTTTACGTTCCAGCCACTc                                                                 <  7:418527/149‑1 (MQ=255)
                                                                                  gAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGAGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTc                                                                 >  6:89987/1‑148 (MQ=255)
                                                                                           ccAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAAc                                                       <  4:298366/149‑1 (MQ=255)
                                                                                             aGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTc                                                     >  1:288385/1‑149 (MQ=255)
                                                                                                 cGGTATCAAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCATCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctc                                                 <  3:429593/149‑1 (MQ=255)
                                                                                                  ggTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgct                                                  >  6:100452/1‑147 (MQ=255)
                                                                                                  ggTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGACAACTACAACCCGCAGGTAGCGGATGATGATGAAGAAGCAGCAACGTATGGCGTTCTGTCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctca                                                >  1:124549/1‑149 (MQ=255)
                                                                                                         aGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCCTCGGCCAGTTTACGTTCCCGCCACTCAAGCTGAAACTCGCTCAGCTcccc                                         >  7:215776/1‑149 (MQ=255)
                                                                                                         aGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCAcc                                         >  1:57458/1‑149 (MQ=255)
                                                                                                         aGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCAcc                                         <  8:498864/149‑1 (MQ=255)
                                                                                                          gTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCg                                        <  2:300107/149‑1 (MQ=255)
                                                                                                              gCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCg                                    <  7:379772/149‑1 (MQ=255)
                                                                                                              gCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCg                                    >  8:252668/1‑149 (MQ=255)
                                                                                                              gCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGt                                                                                  <  1:147477/103‑1 (MQ=255)
                                                                                                              gCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGt                                                                                  >  2:147477/1‑103 (MQ=255)
                                                                                                              gCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGc                                     <  6:169491/148‑1 (MQ=255)
                                                                                                              gCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCg                                    <  4:446331/149‑1 (MQ=255)
                                                                                                                ccGCGTTACGCAGACTGTGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCgg                                   >  7:488304/1‑148 (MQ=255)
                                                                                                                  gcgTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTg                                                                                                <  7:197332/85‑1 (MQ=255)
                                                                                                                  gcgTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTg                                                                                                >  8:197332/1‑85 (MQ=255)
                                                                                                                         gCAGACTGAGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAac                         >  4:371878/1‑149 (MQ=255)
                                                                                                                             aCTGAGTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCt                     <  5:89987/149‑1 (MQ=255)
                                                                                                                              cTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTg                    <  2:240259/149‑1 (MQ=255)
                                                                                                                              cTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTg                    <  7:66471/149‑1 (MQ=255)
                                                                                                                               tgtgTTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTgg                   <  7:448657/149‑1 (MQ=255)
                                                                                                                                gtgtTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCa                                                                                    <  2:254094/83‑1 (MQ=255)
                                                                                                                                gtgtTGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCa                                                                                    >  1:254094/1‑83 (MQ=255)
                                                                                                                                 tgtTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTgg                   <  2:57458/147‑1 (MQ=255)
                                                                                                                                   ttGATCGAGCCAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTAt               >  8:250760/1‑149 (MQ=255)
                                                                                                                                     gATCGAGCCAACTACAACCCGCAGGTAGCGCAACAACATCCAGCAGCACCAACGTAAAGCGTACTGCCGCATCGCCCAGAATACCAACCAGCCACACAAGCTAAAACTCCCACACCACACCCACCGCCACCCCAAACACCACGCTAAcc             >  6:273515/1‑149 (MQ=255)
                                                                                                                                      aTCGAGCCAACTACAACCCGCAGGTAGCTGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATAGGCCAGTTTACGTTCCAGAAACTCAAGCTGAAACTCgc                                                   >  6:177781/1‑110 (MQ=255)
                                                                                                                                            ccAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACa      <  5:121073/149‑1 (MQ=255)
                                                                                                                                            ccAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTGCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACa      <  2:120859/149‑1 (MQ=255)
                                                                                                                                             cAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAg     >  1:406183/1‑149 (MQ=255)
                                                                                                                                             cAACTAAAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAg     >  2:411221/1‑149 (MQ=255)
                                                                                                                                                cTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAt  >  2:158145/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAT  >  NC_000913/3176172‑3176464

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: