Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,603,056 C→T 87.3% W142* (TGG→TGA PP_4924 ← subtilase family serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,603,0560CT87.3% 156.5 / 13.1 63W142* (TGG→TGAPP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base T (26/29);  total (30/33)
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.59e-01

GTGACTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAGG  >  NC_002947/5602911‑5603196
                                                                                                                                                 |                                                                                                                                             
gTGACTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCag                                                                                                                                            >  5:303307/1‑149 (MQ=255)
   aCTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCagta                                                                                                                                          <  4:287166/148‑1 (MQ=255)
   aCTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCagta                                                                                                                                          <  4:145136/148‑1 (MQ=255)
    cTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCagtagt                                                                                                                                        >  2:303320/1‑149 (MQ=255)
          ccTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCtctc                                                                                                                                  >  8:268813/1‑149 (MQ=255)
                  tCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTTATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTc                                                                                                                           <  3:296347/148‑1 (MQ=255)
                   ccACTCCTTAGCAATTAAACGGTTTTCGCCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGa                                                                                                                         >  1:210950/1‑149 (MQ=255)
                     aCTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGa                                                                                                                         >  2:226099/1‑147 (MQ=255)
                     aCTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACa                                                                                                                       <  7:81957/149‑1 (MQ=255)
                          ttAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAAt                                                                                                                  >  7:67331/1‑149 (MQ=255)
                            aGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAAt                                                                                                                  <  6:71501/147‑1 (MQ=255)
                             gCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATcac                                                                                                               <  2:113488/149‑1 (MQ=255)
                                                cggAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAAc                                                                                           <  2:14166/149‑1 (MQ=255)
                                                 ggAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACg                                                                                          >  4:303798/1‑149 (MQ=255)
                                                    agaaACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTc                                                                                        >  2:213299/1‑148 (MQ=255)
                                                    agaaACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTcc                                                                                       <  1:303320/149‑1 (MQ=255)
                                                     gaaACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCa                                                                                      >  3:156303/1‑149 (MQ=255)
                                                       aaCTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATg                                                                                    >  5:11721/1‑149 (MQ=255)
                                                         cgcatGATGGCAATGGCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCGGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGcc                                                                                  >  1:237466/5‑149 (MQ=255)
                                                          tactGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGccc                                                                                 >  2:335551/1‑149 (MQ=255)
                                                               aTGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCctct                                                                              >  4:189208/1‑147 (MQ=255)
                                                                  gcaatggcatcatGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGa                                                                         <  2:237466/142‑1 (MQ=255)
                                                                  ccattgccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAAACACCGGACTCAACGTCCATGCCCTCTCAgg                                                                          <  1:226099/149‑1 (MQ=255)
                                                                  ccattgccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAgg                                                                          <  4:156303/148‑1 (MQ=255)
                                                                   cattgccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCa                                                                            <  7:262487/145‑1 (MQ=255)
                                                                   cattgccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACAGGACTCAAAGTCCATGCACTCTCAgag                                                                         >  8:262487/1‑146 (MQ=255)
                                                                        ccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTc                                                                    <  6:11721/148‑1 (MQ=255)
                                                                        ccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACAAGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTc                                                                    >  4:5439/1‑148 (MQ=255)
                                                                         catcatGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTc                                                                    >  5:326818/1‑147 (MQ=255)
                                                                         catcatGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCaa                                                                  <  6:303307/149‑1 (MQ=255)
                                                                          atcatGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAATGTCCATGCCCTCTCAGGATATTCAAg                                                                 <  2:210950/149‑1 (MQ=255)
                                                                                cGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTc                                                           <  7:268813/149‑1 (MQ=255)
                                                                                 ggTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTc                                                           >  3:17298/1‑148 (MQ=255)
                                                                                     aGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGaaa                                                       <  4:17298/148‑1 (MQ=255)
                                                                                         tgAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTc                                                  >  1:362213/1‑149 (MQ=255)
                                                                                           aaaggaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCaaa                                                                                                        <  2:326996/93‑1 (MQ=255)
                                                                                           aaaggaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCaaa                                                                                                        >  1:326996/1‑93 (MQ=255)
                                                                                           aaaggaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTcc                                                 <  6:326818/148‑1 (MQ=255)
                                                                                           aaaggaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTcc                                                 <  3:303798/148‑1 (MQ=255)
                                                                                           aaaggaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTcc                                                 <  3:189208/148‑1 (MQ=255)
                                                                                               gaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCtt                                               <  1:211040/146‑1 (MQ=255)
                                                                                               gaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGAATCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCtt                                               >  2:211040/1‑146 (MQ=255)
                                                                                                      tAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAg                                                                                                                      >  5:183025/1‑69 (MQ=255)
                                                                                                      tAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAg                                                                                                                      <  6:183025/69‑1 (MQ=255)
                                                                                                           tCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGAc                                <  1:213299/149‑1 (MQ=255)
                                                                                                              tCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACctc                             <  8:67331/149‑1 (MQ=255)
                                                                                                                 gTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGAAATCACATCAAACACAGGACTCAACGTCCATGCCCTCTCAGGATATTAAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCttt                          >  7:41806/1‑149 (MQ=255)
                                                                                                                  ttAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCtttt                         >  2:14296/1‑149 (MQ=255)
                                                                                                                         ccccGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGc                     >  2:9115/1‑146 (MQ=255)
                                                                                                                         ccccGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGc                     <  1:9115/146‑1 (MQ=255)
                                                                                                                          cccGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTc                                                           >  6:253902/1‑107 (MQ=255)
                                                                                                                          cccGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTc                                                           <  5:253902/107‑1 (MQ=255)
                                                                                                                          cccGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGACCTCTCAGCATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACATCTTTTTGGCTTAc                 >  5:159501/1‑149 (MQ=255)
                                                                                                                                 aGGAAAGCACACGCCGTCAGTAGTCCTCTCGATTTTCGACATCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGAc            <  3:5439/147‑1 (MQ=255)
                                                                                                                                  ggAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCttt                          >  6:219837/1‑132 (MQ=255)
                                                                                                                                  ggAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCttt                          <  5:219837/132‑1 (MQ=255)
                                                                                                                                       gCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGcc                                                                                  >  7:208165/1‑71 (MQ=255)
                                                                                                                                       gCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGcc                                                                                  <  8:208165/71‑1 (MQ=255)
                                                                                                                                        cacaTGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAg   <  1:335551/149‑1 (MQ=255)
                                                                                                                                         acaTGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACTCCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAg   <  5:282697/148‑1 (MQ=255)
                                                                                                                                         acaTGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAg   >  6:370510/1‑148 (MQ=255)
                                                                                                                                         acaTGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAg   >  6:282697/1‑148 (MQ=255)
                                                                                                                                         acaTGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAgg  <  2:362213/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GTGACTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAGG  >  NC_002947/5602911‑5603196

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: