Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,763,446 A→T 6.2% M216L (ATG→TTG)  PP_3326 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,763,4460AT 6.2% 195.3 / 3.4 64M216L (ATG→TTG) PP_3326hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base T (2/2);  total (32/32)
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.49e-01

GCAAACCGGCTGCTGATATCCCGCTGCAGGGCGATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCT  >  NC_002947/3763300‑3763590
                                                                                                                                                  |                                                                                                                                                
gCAAACCGGCTGCTGATATCCCGCTGCAGGGCGATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATgg                                                                                                                                               >  2:208032/1‑150 (MQ=255)
     ccGGCTGCTGATATCCCGCTGCAGGGCGATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGAc                                                                                                                                          >  2:1141043/1‑150 (MQ=255)
        gctgctGATATCCCGCTGCAGGGCGATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGc                                                                                                                                       >  2:829639/1‑150 (MQ=255)
        gctgctGATATCCCGCTGCAGGGCGATTACCGAGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGc                                                                                                                                       >  2:16521/1‑150 (MQ=255)
         ctgctgATATCCCGCTGCAGGGCGATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCa                                                                                                                                      <  1:241812/150‑1 (MQ=255)
          tgctgATATCCCGCTGCAGGGCGATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAg                                                                                                                                     >  1:199772/1‑150 (MQ=255)
                      gCTGCAGGGCGATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCca                                                                                                                         <  1:71978/150‑1 (MQ=255)
                             ggCGATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCATTGGCAgggta                                                                                                                                        >  1:737675/1‑124 (MQ=255)
                             ggCGATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCATTGGCAgggta                                                                                                                                        <  2:737675/128‑5 (MQ=255)
                                gATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCATTGGCAgggtag                                                                                                                                       >  2:848491/1‑121 (MQ=255)
                                gATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCATTGGCAgggtag                                                                                                                                       <  1:848491/126‑6 (MQ=255)
                                     ccGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTgttcag                                                                                                          <  1:113878/150‑1 (MQ=255)
                                        gggAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGc                                                                                                                                              >  1:1176419/1‑111 (MQ=255)
                                        gggAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGc                                                                                                                                              <  2:1176419/111‑1 (MQ=255)
                                             ccccTCACCAGGTCAGCGCCACCACCGACAAGCACGGCCGGGCACGTCTTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCCATGGCAGACGGCAGCCCCGTGGCCCAGCAGAACCTGTTCAGCTCACTGt                                                                                                  >  2:54298/1‑150 (MQ=255)
                                               ccccACGAGGTCAGCGCCACCACCGACGAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCg                                                                                                >  2:1013308/4‑150 (MQ=255)
                                                   aCCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTcc                                                                                            <  1:829639/150‑1 (MQ=255)
                                                       ggTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACTCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTggg                                                                                        <  1:824175/150‑1 (MQ=255)
                                                       ggTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTggg                                                                                        >  1:539270/1‑150 (MQ=255)
                                                                 accaccGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTgttcagtt                                                                                                        >  1:1142397/1‑124 (MQ=255)
                                                                 accaccGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTgttcagtt                                                                                                        <  2:1142397/124‑1 (MQ=255)
                                                                   caccGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGcacca                                                                            >  1:141269/1‑150 (MQ=255)
                                                                    accGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGcaccac                                                                           >  2:178998/1‑150 (MQ=255)
                                                                      cGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGa                                                                         >  2:455105/1‑150 (MQ=255)
                                                                      cGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGa                                                                         >  2:570059/1‑150 (MQ=255)
                                                                           agcaagGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTaaa                                                                    <  1:453352/150‑1 (MQ=255)
                                                                             caagGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACt                                                                  <  2:539270/150‑1 (MQ=255)
                                                                              aagGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTg                                                                                                                              >  2:1149441/1‑89 (MQ=255)
                                                                              aagGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTg                                                                                                                              <  1:1149441/89‑1 (MQ=255)
                                                                                ggCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAg                                                               >  2:376657/1‑150 (MQ=255)
                                                                                        aCGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATAc                                                       <  1:1013308/150‑1 (MQ=255)
                                                                                        aCGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATAc                                                       >  2:170735/1‑150 (MQ=255)
                                                                                        aCGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATAc                                                       >  2:797043/1‑150 (MQ=255)
                                                                                           tATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTgg                                                                                         <  1:278494/113‑1 (MQ=255)
                                                                                           tATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTgg                                                                                         >  2:278494/1‑113 (MQ=255)
                                                                                           tATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTc                                                    >  1:476825/1‑150 (MQ=255)
                                                                                                     cGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACg                                          <  1:208032/150‑1 (MQ=255)
                                                                                                         aTGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCg                                      <  1:16521/150‑1 (MQ=255)
                                                                                                          tGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTg                                                                                                                              >  2:1183426/1‑61 (MQ=255)
                                                                                                          tGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTg                                                                                                                              <  1:1183426/61‑1 (MQ=255)
                                                                                                           gAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCGTCATGGCCGAg                                    <  2:1058816/150‑1 (MQ=255)
                                                                                                            aaGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCGTCATGGCCGAgg                                   <  2:476825/150‑1 (MQ=255)
                                                                                                             aGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGt                                  >  2:1000787/1‑150 (MQ=255)
                                                                                                                 cTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGc                                                              <  2:923959/118‑1 (MQ=255)
                                                                                                                 cTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGc                                                              >  1:923959/1‑118 (MQ=255)
                                                                                                                  tCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTaaa                                                                    >  1:1035925/1‑111 (MQ=255)
                                                                                                                  tCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTaaa                                                                    <  2:1035925/111‑1 (MQ=255)
                                                                                                                    aaCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGa                                            >  1:730561/1‑133 (MQ=255)
                                                                                                                    aaCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGa                                            <  2:730561/133‑1 (MQ=255)
                                                                                                                     aCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGt                          >  1:593190/1‑150 (MQ=255)
                                                                                                                      catcatCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGtt                         <  1:1141043/150‑1 (MQ=255)
                                                                                                                       atcatcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCATCATAGCTTGTGCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTg                        <  2:495457/150‑1 (MQ=255)
                                                                                                                         catcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTTTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGa                      >  2:941768/1‑150 (MQ=255)
                                                                                                                          atcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                     <  1:570059/150‑1 (MQ=255)
                                                                                                                                gCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTc               <  1:455105/150‑1 (MQ=255)
                                                                                                                                 cccAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCg              <  1:170735/150‑1 (MQ=255)
                                                                                                                                         aCCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGt      <  2:99481/150‑1 (MQ=255)
                                                                                                                                           ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGt      <  1:1076642/148‑1 (MQ=255)
                                                                                                                                           ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGt      >  2:1076642/1‑148 (MQ=255)
                                                                                                                                           ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTCAGTATGCCATACCGCTAGTAAACTAAGCCACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGc    <  1:797043/150‑1 (MQ=255)
                                                                                                                                             tGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt  <  2:823517/150‑1 (MQ=255)
                                                                                                                                             tGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt  >  1:823517/1‑150 (MQ=255)
                                                                                                                                             tGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt  <  2:257400/150‑1 (MQ=255)
                                                                                                                                             tGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGACGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt  >  1:269937/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCAAACCGGCTGCTGATATCCCGCTGCAGGGCGATTACCGGGGAGCCCCTCACCAGGTCAGCGCCACCACCGACAAGCAAGGCCGGGCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCT  >  NC_002947/3763300‑3763590

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