Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 431,355 A→G 6.6% L320P (CTG→CCG)  PP_0354 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947431,3550AG 6.6% 223.3 / 4.4 76L320P (CTG→CCG) PP_0354CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base G (3/2);  total (38/38)
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.83e-01

GCCCGGGTCGCCACGTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATGGGCAATATGCCG  >  NC_002947/431210‑431503
                                                                                                                                                  |                                                                                                                                                    
gCCCGGGTCGCCACGTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGc                                                                                                                                                    >  5:155309/1‑148 (MQ=255)
gCCCGGGTCGCCACGTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGc                                                                                                                                                    >  2:234979/1‑148 (MQ=255)
    gggTCGCCACGTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTTATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAAc                                                                                                                                               <  4:233406/149‑1 (MQ=255)
    gggTCGCCACGTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAAc                                                                                                                                               <  1:234979/149‑1 (MQ=255)
              gTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATc                                                                                                                                      >  2:365314/1‑148 (MQ=255)
                 cGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAg                                                                                                                                   <  6:254796/148‑1 (MQ=255)
                 cGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAg                                                                                                                                   <  6:277048/148‑1 (MQ=255)
                  ggccagggccagTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGc                                                                                                                                 <  1:284327/149‑1 (MQ=255)
                   gccagggccagTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCg                                                                                                                                >  4:341063/1‑149 (MQ=255)
                   gccagggccagTTCGATGACCCGGCACAAGGTATGGTCGTTGAGAAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCg                                                                                                                                >  8:185039/1‑149 (MQ=255)
                        ggccagTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTAgcg                                                                                                                           <  7:188382/149‑1 (MQ=255)
                        ggccagTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTAgcg                                                                                                                           >  8:188382/1‑149 (MQ=255)
                             gTTCGATGACCCGGCACACGGTATGGGCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATgg                                                                                                                      <  7:163129/149‑1 (MQ=255)
                              ttCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGc                                                                                                                     >  6:345340/1‑149 (MQ=255)
                               tCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGc                                                                                                                     <  4:47890/148‑1 (MQ=255)
                                 gATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATggaga                                                                                                                   >  3:20200/1‑145 (MQ=255)
                                 gATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATg                                                                                                                       <  4:20200/144‑1 (MQ=255)
                                 gATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATg                                                                                                                       >  4:263135/1‑144 (MQ=255)
                                 gATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATg                                                                                                                       <  3:263135/144‑1 (MQ=255)
                                         ggCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTgcgc                                                                                                          <  6:192019/149‑1 (MQ=255)
                                          gCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTgcgcg                                                                                                         >  8:350055/1‑149 (MQ=255)
                                               cGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAgg                                                                                                    >  6:331030/1‑149 (MQ=255)
                                                ggTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGt                                                                                                   <  6:155309/149‑1 (MQ=255)
                                                      gTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAAcc                                                                                              <  7:185039/148‑1 (MQ=255)
                                                      gTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAAcc                                                                                              <  3:341063/148‑1 (MQ=255)
                                                      gTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAAc                                                                                               >  5:181336/1‑147 (MQ=255)
                                                        cGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGa                                                                                           <  7:221425/149‑1 (MQ=255)
                                                          ttGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATc                                                                                         <  8:315433/149‑1 (MQ=255)
                                                             aGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGAc                                                                                      <  1:365314/149‑1 (MQ=255)
                                                             aGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGAc                                                                                      >  4:136396/1‑149 (MQ=255)
                                                               cAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGAccc                                                                                    <  6:181336/149‑1 (MQ=255)
                                                                aaGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGGTTTCGCTGCGCGGCGAAATCAGcc                                                                                                                                          <  3:329858/94‑2 (MQ=38)
                                                                aaGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGGTTTCGCTGCGCGGCGAAATCAGcc                                                                                                                                          >  4:329858/1‑93 (MQ=38)
                                                                      aTGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAg                                                                             <  3:312032/149‑1 (MQ=255)
                                                                                   tCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTcg                                                                >  6:65253/1‑149 (MQ=255)
                                                                                   tCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTc                                                                 <  3:136396/148‑1 (MQ=255)
                                                                                    cTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTcgc                                                               >  7:53367/1‑149 (MQ=255)
                                                                                      gCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCGCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTcgcg                                                              >  1:280967/1‑149 (MQ=255)
                                                                                        gTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTc                                                           <  5:331030/149‑1 (MQ=255)
                                                                                         tCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCg                                                          <  5:345340/149‑1 (MQ=255)
                                                                                         tCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCg                                                          <  1:136623/149‑1 (MQ=255)
                                                                                            acgaTGTGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGaa                                                       <  5:81378/149‑1 (MQ=255)
                                                                                                tGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCa                                                   >  5:222590/1‑149 (MQ=255)
                                                                                                   gcCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCAcgc                                                >  5:199358/1‑149 (MQ=255)
                                                                                                     cATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCAcgccg                                              >  5:357821/1‑149 (MQ=255)
                                                                                                            cAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGaa                                                       <  2:186019/133‑1 (MQ=255)
                                                                                                            cAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGaa                                                       >  1:186019/1‑133 (MQ=255)
                                                                                                            cAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAg                                         >  6:154482/1‑147 (MQ=255)
                                                                                                             aGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCg                                      >  6:143537/1‑149 (MQ=255)
                                                                                                                aTGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTcc                                   <  6:222590/149‑1 (MQ=255)
                                                                                                                 tGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCg                                  >  3:198429/1‑149 (MQ=255)
                                                                                                                   cGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGtt                                <  8:53367/149‑1 (MQ=255)
                                                                                                                    gTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGttt                               <  2:280967/149‑1 (MQ=255)
                                                                                                                    gTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGttt                               >  3:88583/1‑149 (MQ=255)
                                                                                                                     ttCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGAc                                                                                      >  1:297537/1‑93 (MQ=255)
                                                                                                                     ttCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGAc                                                                                      <  2:297537/93‑1 (MQ=255)
                                                                                                                     ttCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCg                                  <  4:235132/145‑1 (MQ=255)
                                                                                                                     ttCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCg                                  >  3:235132/1‑145 (MQ=255)
                                                                                                                     ttCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGtttt                              <  4:198429/149‑1 (MQ=255)
                                                                                                                       cTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGttt                               >  4:318137/1‑146 (MQ=255)
                                                                                                                       cTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGttt                               <  3:318137/146‑1 (MQ=255)
                                                                                                                        tACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTcca                           >  3:30195/1‑149 (MQ=255)
                                                                                                                          cAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTccacc                         <  4:88583/149‑1 (MQ=255)
                                                                                                                           aaGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTccac                          >  3:122237/1‑147 (MQ=255)
                                                                                                                           aaGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTcaac                          >  6:114614/1‑147 (MQ=255)
                                                                                                                             gctggctgATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAgg                      >  5:352894/1‑149 (MQ=255)
                                                                                                                               tggctgATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAgg                      >  7:129182/1‑147 (MQ=255)
                                                                                                                               acgctgGTTTCGCTGCGCGGCGAAATCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAgg                      >  1:166118/3‑147 (MQ=255)
                                                                                                                                ggctgATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAgg                      >  4:185423/1‑146 (MQ=255)
                                                                                                                                 gctgGTTTCGCTGCGCGGCGAAATCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGc                                       <  3:342521/128‑1 (MQ=255)
                                                                                                                                 gctgGTTTCGCTGCGCGGCGAAATCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGc                                       >  4:342521/1‑128 (MQ=255)
                                                                                                                                   tgATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAAc                                                     >  3:130662/1‑112 (MQ=255)
                                                                                                                                   tgATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAAc                                                     <  4:130662/112‑1 (MQ=255)
                                                                                                                                      tttCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATgg             >  3:103833/1‑149 (MQ=255)
                                                                                                                                      tttCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCAGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATgg             <  7:350055/149‑1 (MQ=255)
                                                                                                                                               cgcAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATGGGCAATATg     <  6:352894/149‑1 (MQ=255)
                                                                                                                                                  aGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATGGGCAATATGCCg  >  6:223110/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCCGGGTCGCCACGTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATGGGCAATATGCCG  >  NC_002947/431210‑431503

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