Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,676,800 A→T 5.7% F30Y (TTC→TAC)  PP_3241 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,676,8000AT 5.7% 214.7 / 3.7 71F30Y (TTC→TAC) PP_3241hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base T (2/2);  total (37/34)
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.79e-01

GGGCAATGTCGCCGAGTTTCTGGTTGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCG  >  NC_002947/3676651‑3676948
                                                                                                                                                     |                                                                                                                                                    
gggCAATGTCGCCGAGTTTCTGGTTGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATcgacga                                                                                                                                                      >  2:2234324/1‑150 (MQ=255)
  gCAATGTCGCCGAGTTTCTGGTTGGCCTCGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGaaa                                                                                                                                                    >  2:1252892/1‑150 (MQ=255)
      tGTCGCCGAGTTTCTGGTTGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAatcat                                                                                                                                                <  1:2172926/150‑1 (MQ=255)
      tGTCGCCGAGTTTCTGGTTGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAatcat                                                                                                                                                <  2:916356/150‑1 (MQ=255)
                tttCTGGTTGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGaa                                                                                                                                                     >  1:591584/1‑135 (MQ=255)
                tttCTGGTTGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGaa                                                                                                                                                     <  2:591584/135‑1 (MQ=255)
                  tCTGGTTGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCgg                                                                                                                                    <  1:770757/150‑1 (MQ=255)
                    tGGTTGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAatcatc                                                                                                                                               <  1:506520/137‑1 (MQ=255)
                    tGGTTGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAatcatc                                                                                                                                               >  2:506520/1‑137 (MQ=255)
                       ttGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCAtt                                                                                                                               >  2:343320/1‑150 (MQ=255)
                                 tgtgGTCGGTGACCATTTGCTGTGCGAATGACTTGATCTCGGCAAAACCTGTATTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCg                                                                                                                     >  1:2275572/1‑150 (MQ=255)
                                  gtgGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGa                                                                                                                    <  2:1803227/150‑1 (MQ=255)
                                   tgGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGACGTCGCGGCCATTGCGTGCACCGAt                                                                                                                   <  1:1298004/150‑1 (MQ=255)
                                        ggTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAg                                                                                                              <  1:343320/150‑1 (MQ=255)
                                         gTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCTACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAgg                                                                                                             >  1:144387/1‑150 (MQ=255)
                                         gTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAgg                                                                                                             >  1:1158421/1‑150 (MQ=255)
                                           gACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGcc                                                                                                           <  2:1158421/150‑1 (MQ=255)
                                              cATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCaa                                                                                                        >  1:2028915/1‑150 (MQ=255)
                                               aTTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGAt                                                                                                                   <  2:482698/138‑1 (MQ=255)
                                               aTTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGAt                                                                                                                   >  1:482698/1‑138 (MQ=255)
                                                  tGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGCCCATTGCGTGCACCGATGCCAGGCCCAACAcc                                                                                                    <  2:1344510/150‑1 (MQ=255)
                                                   gCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCa                                                                                                   >  1:1057910/1‑150 (MQ=255)
                                                     tGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATg                                                                                                 >  1:1775587/1‑150 (MQ=255)
                                                      gggCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATg                                                                                                 >  1:55678/1‑149 (MQ=255)
                                                      gggCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATg                                                                                                 <  2:55678/149‑1 (MQ=255)
                                                       ggCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCAtt                                                                                                                               >  1:1699497/1‑118 (MQ=255)
                                                       ggCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCAtt                                                                                                                               <  2:1699497/118‑1 (MQ=255)
                                                         cGAATGACTTGATCTCGGCGTGCCGGCCTATCTCCTGCGCCAGGTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGCAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGaaa                                                                                             <  2:2144153/150‑1 (MQ=255)
                                                          gAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCACGGaaaa                                                                                            >  1:209744/1‑150 (MQ=255)
                                                          gAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAACTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGaaaa                                                                                            <  1:1128695/150‑1 (MQ=255)
                                                           aaTGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGaaaaa                                                                                           <  1:2026130/150‑1 (MQ=255)
                                                           aaTGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCACAACACCATGGaaaaa                                                                                           <  2:144387/150‑1 (MQ=255)
                                                           aaTGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATACCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCAt                                                                                                                                >  1:1178538/1‑113 (MQ=255)
                                                           aaTGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATACCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCAt                                                                                                                                <  2:1178538/113‑1 (MQ=255)
                                                            aTGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAAc                                                                                          <  1:1504167/150‑1 (MQ=255)
                                                                  tGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAAcgc                                                                                    >  1:1465952/1‑150 (MQ=255)
                                                                       tcGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTcgcg                                                                                                                                     >  1:1411884/1‑96 (MQ=255)
                                                                       tcGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTcgcg                                                                                                                                     <  2:1411884/96‑1 (MQ=255)
                                                                       tcGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCtg                                                                               >  1:72192/1‑150 (MQ=255)
                                                                       tcGGCGCTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGCAATCATCGGAAGGCGCGGCCATTGCGTGCACCGAAGCCAGGCCCCACACCATTGAACAACCAAcgctgctt                                                                               >  2:877534/1‑145 (MQ=255)
                                                                         ggCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCtgat                                                                             >  2:229581/1‑150 (MQ=255)
                                                                         ggCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCtgat                                                                             >  2:2153732/1‑150 (MQ=255)
                                                                            tttCTGGGTTTTCTCCTGGGCCAGCTTGCCGGTCACCACTTCGGCCATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGaa                                                                          <  2:111719/149‑1 (MQ=255)
                                                                             ttCNGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGACATCGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAag                                                                         >  2:1518056/1‑150 (MQ=255)
                                                                              tCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAaga                                                                        <  1:1252892/150‑1 (MQ=255)
                                                                                    ttttCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTAt                                                                  <  2:1775587/150‑1 (MQ=255)
                                                                                    ttttCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTAGTGATCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTAt                                                                  >  1:61715/1‑150 (MQ=255)
                                                                                      ttCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGTAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTc                                                                >  2:937609/1‑150 (MQ=255)
                                                                                      ttCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGTAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTc                                                                >  2:2171626/1‑150 (MQ=255)
                                                                                        ctcCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCAt                                                              <  1:2153732/150‑1 (MQ=255)
                                                                                          ccTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGCCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGc                                                            <  2:1661932/150‑1 (MQ=255)
                                                                                             gAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTc                                                         <  1:357724/150‑1 (MQ=255)
                                                                                                ccAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACt                                                      >  1:2294443/1‑150 (MQ=255)
                                                                                                    cTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAACCCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCtt                                                  <  2:2000141/150‑1 (MQ=255)
                                                                                                        ccGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTaaa                                              >  2:670922/1‑150 (MQ=255)
                                                                                                        ccGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTaaa                                              <  2:1057910/150‑1 (MQ=255)
                                                                                                         cGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTc                                                     >  2:222967/1‑142 (MQ=255)
                                                                                                         cGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTc                                                     <  1:222967/142‑1 (MQ=255)
                                                                                                             caccacTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTg                                         >  1:465566/1‑150 (MQ=255)
                                                                                                             caccacTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGACCCAACACCATGGAAAAACCAACGCGCCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTg                                         >  2:1120591/1‑150 (MQ=255)
                                                                                                                cacTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGgcg                                      >  1:851796/1‑150 (MQ=255)
                                                                                                                cacTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGgc                                       >  2:234867/1‑149 (MQ=255)
                                                                                                                cacTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGgc                                       <  1:234867/149‑1 (MQ=255)
                                                                                                                 acTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGTAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGgcgc                                     <  1:937609/150‑1 (MQ=255)
                                                                                                                 acTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGTAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGgcgc                                     <  1:2171626/150‑1 (MQ=255)
                                                                                                                       gCAATGTCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTc                               >  1:1299292/1‑150 (MQ=255)
                                                                                                                       gCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCa                                                        >  2:1195629/1‑125 (MQ=255)
                                                                                                                       gCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCa                                                        <  1:1195629/125‑1 (MQ=255)
                                                                                                                        cAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGGAAAAATTTATTCTTCACTCCGTGAAACCGTGGCGCGACCTCa                              >  2:1840940/1‑150 (MQ=255)
                                                                                                                                gCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATg                                                             <  2:1459632/111‑1 (MQ=255)
                                                                                                                                gCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATg                                                             >  1:1459632/1‑111 (MQ=255)
                                                                                                                                gCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTcgcg                      >  2:589475/1‑150 (MQ=255)
                                                                                                                                   tCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCACGGAAAAACTAACGCGTCTGATGAAAGGACTGGTAGATCGGCAGTGCTTTAAACCGTGGCGCGACCTCAATGTCGCGTAc                   >  1:1116379/1‑150 (MQ=255)
                                                                                                                                         gCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGCTCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGa             <  1:1840940/150‑1 (MQ=255)
                                                                                                                                           tGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAgg           >  2:13686/1‑150 (MQ=255)
                                                                                                                                                   cgaAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGc   <  2:72192/150‑1 (MQ=255)
                                                                                                                                                    gaAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCg  >  2:286907/1‑150 (MQ=255)
                                                                                                                                                    gaAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCg  >  1:631811/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGCAATGTCGCCGAGTTTCTGGTTGGCCTGGGTGTGGTCGGTGACCATTTGCTGGGCGAATGACTTGATCTCGGCGTTCTGGCTTTTCTCCTGAGCCAGCTTGCCGGTCACCACTTCGGCAATGCCGGCCTCGGTGGCTGCATCGACGAAATCATCGGAAGTCGCGGCCATTGCGTGCACCGATGCCAGGCCCAACACCATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCG  >  NC_002947/3676651‑3676948

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