Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,389,595 T→A 6.8% Q28H (CAA→CAT)  PP_2991 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,389,5950TA 6.8% 208.9 / 3.6 73Q28H (CAA→CAT) PP_2991hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/29);  new base A (3/2);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CCGGCAGCTCTGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTAC  >  NC_002947/3389447‑3389738
                                                                                                                                                    |                                                                                                                                               
ccGGCAGCTCTGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCt                                                                                                                                                 >  4:135832/1‑149 (MQ=255)
  ggCAGCTCTTCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCttg                                                                                                                                               <  3:237818/149‑1 (MQ=255)
    cAGGTCTGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCttggt                                                                                                                                             <  1:171335/149‑1 (MQ=255)
     aGCTCTGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCttggtt                                                                                                                                            <  4:144615/149‑1 (MQ=255)
      gCTCTGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCttggttg                                                                                                                                           <  8:148258/149‑1 (MQ=255)
       ctctGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCttggttgg                                                                                                                                          <  4:114957/149‑1 (MQ=255)
        agagctgcAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCttggttgg                                                                                                                                          >  3:311202/4‑148 (MQ=255)
            ctgcAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAAt                                                                                                                                      <  6:293682/148‑1 (MQ=255)
                 aTCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACCACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCt                                                                                                                                <  1:309808/149‑1 (MQ=255)
                  tCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCt                                                                                                                                >  7:111742/1‑148 (MQ=255)
                    aCCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGaa                                                                                                                             >  3:28979/1‑149 (MQ=255)
                     ccGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGaa                                                                                                                             >  1:172823/1‑148 (MQ=255)
                      cGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGaa                                                                                                                             >  5:134036/1‑147 (MQ=255)
                      cGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACa                                                                                                                           <  6:264631/149‑1 (MQ=255)
                            gggCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATc                                                                                                                     >  4:300105/1‑149 (MQ=255)
                            gggCGGGTTCGAGCGATAACGGTCGTACAGCCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGAGTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATc                                                                                                                     >  2:106324/1‑149 (MQ=255)
                             ggCGGGTTAGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACAGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGACGACGTGATCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATcc                                                                                                                    >  8:84654/1‑149 (MQ=255)
                               cGGGTTCGCGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCGTTGGTTGGCAATTGAGCTGAACAGCAATCCTg                                                                                                                  >  4:138326/1‑149 (MQ=255)
                                gggTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGt                                                                                                                 <  6:183024/149‑1 (MQ=255)
                                   ttCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTaa                                                                                                              >  8:152607/1‑149 (MQ=255)
                                        gCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCgg                                                                                                         >  2:296185/1‑149 (MQ=255)
                                         cGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGt                                                                                                        >  1:206373/1‑149 (MQ=255)
                                          gATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCAAGGGTTGGcgg                                                                                                                                       <  2:2781/117‑3 (MQ=255)
                                          gATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCAAGGGTTGGcgg                                                                                                                                       >  1:2781/1‑115 (MQ=255)
                                             aaCGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCtt                                                                                                    >  5:356662/1‑149 (MQ=255)
                                                 gTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCaa                                                                                                <  2:99223/149‑1 (MQ=255)
                                                 ctCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGGCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCATGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCaa                                                                                                <  2:94214/148‑1 (MQ=255)
                                                  tCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAAt                                                                                                             <  5:26994/135‑1 (MQ=255)
                                                  tCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAAt                                                                                                             >  6:26994/1‑135 (MQ=255)
                                                     tACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGtt                                                                                             >  3:105589/1‑148 (MQ=255)
                                                        aGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCaa                                                                                         >  1:358708/1‑149 (MQ=255)
                                                                     cgacCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGt                                                                                    >  3:127737/1‑141 (MQ=255)
                                                                     cgacCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGt                                                                                    <  4:127737/141‑1 (MQ=255)
                                                                      gacCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCcg                                                                           >  1:377175/1‑149 (MQ=255)
                                                                       acCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCcgc                                                                          <  3:300105/149‑1 (MQ=255)
                                                                             aCCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTgg                                                                    >  5:106645/1‑149 (MQ=255)
                                                                             aCCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTgg                                                                    <  4:105589/149‑1 (MQ=255)
                                                                             aCCGGGGGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTgg                                                                    <  7:152607/149‑1 (MQ=255)
                                                                                ggggggCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGaa                                                                 <  7:13428/149‑1 (MQ=255)
                                                                                gggTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGaa                                                                 <  3:135832/149‑1 (MQ=255)
                                                                                gggTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGaa                                                                 >  2:228457/1‑149 (MQ=255)
                                                                                      ccTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGc                                                                                                  >  5:117691/1‑110 (MQ=255)
                                                                                      ccTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGc                                                                                                  <  6:117691/110‑1 (MQ=255)
                                                                                      ccTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACtt                                                           >  3:298087/1‑149 (MQ=255)
                                                                                      ccTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACtt                                                           >  5:378712/1‑149 (MQ=255)
                                                                                       cTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACtt                                                           >  7:185656/1‑148 (MQ=255)
                                                                                        tGCGATTTCTATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAAACCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGTGTGGGAACAACtt                                                           >  2:210623/1‑147 (MQ=255)
                                                                                        tGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGCTCTTTCTCCGAGGGTGGGAACAACtt                                                           >  4:91339/1‑147 (MQ=255)
                                                                                         gCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAg                                                        >  6:294070/1‑149 (MQ=255)
                                                                                           gATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAggg                                                      <  8:185656/149‑1 (MQ=255)
                                                                                                    tgtgTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGa                                              <  2:172823/148‑1 (MQ=255)
                                                                                                     gtgtTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAg                                             <  2:377175/148‑1 (MQ=255)
                                                                                                     gtgtTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAg                                             <  4:227664/148‑1 (MQ=255)
                                                                                                     gtgtTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGc                                            >  7:60317/1‑149 (MQ=255)
                                                                                                      tgtTGGTTTGAGTCGACGTGAGGATTTCAGGCGTCCCAGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGc                                            <  4:28979/148‑1 (MQ=255)
                                                                                                       gtTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCg                                                                                      <  5:589/105‑1 (MQ=255)
                                                                                                       gtTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCg                                                                                      >  6:589/1‑105 (MQ=255)
                                                                                                        ttGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGaacaa                                                              >  8:11908/1‑128 (MQ=255)
                                                                                                        ttGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGaacaa                                                              <  7:11908/128‑1 (MQ=255)
                                                                                                           gTTTCAGTCAACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGtt                                      <  2:206373/149‑1 (MQ=255)
                                                                                                               cAGTCGACGTGAGCATTTCATGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGATTCAATCGGTCTTTCTCCGCGGGTGTGAACAACTTGAGGGTTCAGCGAGCTACGTTGTc                                   <  3:138326/148‑1 (MQ=255)
                                                                                                               cAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTc                                   <  4:311202/148‑1 (MQ=255)
                                                                                                                  tCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGc                                                                                                      >  5:24841/1‑78 (MQ=255)
                                                                                                                  tCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGc                                                                                                      <  6:24841/78‑1 (MQ=255)
                                                                                                                    gACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCaga                             <  6:120408/149‑1 (MQ=255)
                                                                                                                        tGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCagagag                          >  1:66379/1‑148 (MQ=255)
                                                                                                                         gAGCATTTCAGGCGTCCCCGCCAACCCTTGGTGGGCAATTGAGCTGAACAGCAAACCTGTTAATGCGGTGCTTGCCAGTTTCCATCGGTCTGTCTCCGCGGGTGGGAACAACTAGAGGGTTCCGCGCGCTCCGATGGCGCCAGAGCGCa                        >  4:257130/1‑149 (MQ=255)
                                                                                                                            cATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGCACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTAAGTTGTCGCCAGAGAGCACCg                     >  5:57151/1‑149 (MQ=255)
                                                                                                                            cATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCg                     <  2:358708/149‑1 (MQ=255)
                                                                                                                              tttCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGc                    >  7:272072/1‑148 (MQ=255)
                                                                                                                                  aGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGCTCTTTCTCCGCGGCTGGGAACAACTTGTGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCa               >  2:243263/1‑149 (MQ=255)
                                                                                                                                    gCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCAcc             >  2:275898/1‑149 (MQ=255)
                                                                                                                                         cccGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTc        <  6:134036/149‑1 (MQ=255)
                                                                                                                                         cacgacAACCCTTTGTTGGCAATTGAGCTTAACAGAAATGCTTTTAATGCGTTGCTTGCAAGTTTCTATCGGTCTTTCTCCGCGTGTGGGAACAACTTGATTGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGACATCACCGCTTc        <  8:60317/144‑1 (MQ=255)
                                                                                                                                          ttgCCAACCAAGGGTTGGCAATTGAGCTGAACCGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGATACGTTGTCGCAAGAGAGCACCGCCATCACCGCtt         >  5:200075/3‑147 (MQ=255)
                                                                                                                                          ttgCCAACCAAGGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCtt         >  3:184709/3‑147 (MQ=255)
                                                                                                                                             ccAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGc     >  8:79745/1‑148 (MQ=255)
                                                                                                                                               aaCCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTAc  >  5:74882/1‑149 (MQ=255)
                                                                                                                                                 cccTTGTTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGc                                                >  3:286991/1‑101 (MQ=255)
                                                                                                                                                 cccTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGc                                                <  4:286991/101‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCGGCAGCTCTGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTACAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTAC  >  NC_002947/3389447‑3389738

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: ‑