Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,632,128 G→T 8.0% E327* (GAG→TAG) ‡ PP_5590 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,632,1280GT 8.0% 221.1 / 3.7 70E327* (GAG→TAG) ‡PP_5590hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/35);  new base T (2/4);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGGCTG  >  NC_002947/3631979‑3632275
                                                                                                                                                     |                                                                                                                                                   
cacCGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCg                                                                                                                                                     >  1:631972/1‑150 (MQ=255)
   tgACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAgg                                                                                                                                                  >  1:1160791/2‑150 (MQ=255)
         cGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCgctgc                                                                                                                                            <  2:631972/150‑1 (MQ=255)
           gCCACCGCCGGCGTCATTAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGctgctg                                                                                                                                          <  1:556022/150‑1 (MQ=255)
               tggccgGGGTCGTGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGccc                                                                                                                                      <  2:411678/148‑1 (MQ=255)
               tggccgGGGTCGTGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGccc                                                                                                                                      <  2:411591/148‑1 (MQ=255)
                       gTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACACCGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCAc                                                                                                                              <  1:395889/150‑1 (MQ=255)
                                         cTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATc                                                                                                            >  1:1182710/1‑150 (MQ=255)
                                                  ttCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGTTGTTCACGGCCAGCAACGCCTAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAg                                                                                                   <  2:682384/150‑1 (MQ=255)
                                                    cAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCGGGAgaac                                                                                                <  2:171987/150‑4 (MQ=255)
                                                     aaCGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTg                                                                                                >  2:369570/1‑150 (MQ=255)
                                                     aaCGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTg                                                                                                >  2:1389458/1‑150 (MQ=255)
                                                     aaCGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTg                                                                                                >  1:58605/1‑150 (MQ=255)
                                                     aaCGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGGGCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTg                                                                                                >  1:1342993/1‑150 (MQ=255)
                                                          tCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCg                                                                                           >  2:975946/1‑150 (MQ=255)
                                                             ccTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgt                                                                                        <  2:1047730/150‑1 (MQ=255)
                                                             ccTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgt                                                                                        <  1:369570/150‑1 (MQ=255)
                                                              cTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtg                                                                                       <  2:784335/150‑1 (MQ=255)
                                                              cTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtg                                                                                       >  2:1294367/1‑150 (MQ=255)
                                                                 gACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtgt                                                                                      <  1:1315460/148‑1 (MQ=255)
                                                                 gACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtgt                                                                                      >  2:1315460/1‑148 (MQ=255)
                                                                     gCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTgtgggtgg                                                                                >  2:815316/1‑150 (MQ=255)
                                                                        tCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGtggtgg                                                                             <  1:975946/150‑1 (MQ=255)
                                                                         cGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGa                                                                            <  1:304466/150‑1 (MQ=255)
                                                                            gTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGcc                                                                                                               <  1:310044/112‑1 (MQ=255)
                                                                            gTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGcc                                                                                                               >  2:310044/1‑112 (MQ=255)
                                                                               gAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTGGcg                                                                                                                                       >  2:1364213/1‑84 (MQ=37)
                                                                               gAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTGGcg                                                                                                                                       <  1:1364213/85‑2 (MQ=37)
                                                                                 aGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCTATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGc                                                                    <  2:281277/150‑1 (MQ=255)
                                                                                       tggGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTgaaga                                                              >  2:871638/1‑150 (MQ=255)
                                                                                          gCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCTAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATg                                                           <  1:430966/150‑1 (MQ=255)
                                                                                          gCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATg                                                           <  1:871638/150‑1 (MQ=255)
                                                                                          gCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATg                                                           <  2:1182710/150‑1 (MQ=255)
                                                                                           cAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtg                                                                                         >  2:210880/1‑119 (MQ=255)
                                                                                           cAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtg                                                                                         <  1:210880/119‑1 (MQ=255)
                                                                                           cAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATgg                                                          <  2:109003/150‑1 (MQ=255)
                                                                                            aagaagGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGa                                                                                                    <  1:1514421/107‑1 (MQ=255)
                                                                                            aagaagGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGa                                                                                                    >  2:1514421/1‑107 (MQ=255)
                                                                                               aagGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAg                                                      >  2:1404499/1‑150 (MQ=255)
                                                                                                agGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtgt                                                                                      <  2:711753/117‑1 (MQ=255)
                                                                                                agGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtgt                                                                                      >  1:711753/1‑117 (MQ=255)
                                                                                                agGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAgc                                                     >  1:748088/1‑150 (MQ=255)
                                                                                                  gagaTGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAgcgc                                                   <  1:103207/150‑1 (MQ=255)
                                                                                                         aCATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtgt                                                                                      >  2:1205245/1‑108 (MQ=255)
                                                                                                         aCATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtgt                                                                                      <  1:1205245/108‑1 (MQ=255)
                                                                                                           aTCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATg                                          >  1:1300442/1‑150 (MQ=255)
                                                                                                               ccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGcc                                        >  2:641400/1‑148 (MQ=255)
                                                                                                               ccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGcc                                        <  1:641400/148‑1 (MQ=255)
                                                                                                                  aCAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGCGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCTATGTGATGCCGCGTa                                   >  2:1140400/1‑150 (MQ=255)
                                                                                                                       gCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAAt                                                            >  2:236978/1‑120 (MQ=255)
                                                                                                                       gCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAAt                                                            <  1:236978/120‑1 (MQ=255)
                                                                                                                             aaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCa                                <  2:1643775/142‑1 (MQ=255)
                                                                                                                             aaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCa                                >  1:1643775/1‑142 (MQ=255)
                                                                                                                               ggTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtgt                                                                                      <  1:1379815/86‑1 (MQ=255)
                                                                                                                               ggTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtgt                                                                                      >  2:1379815/1‑86 (MQ=255)
                                                                                                                                  gtTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGCGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTcga                   <  1:1140400/150‑1 (MQ=255)
                                                                                                                                   ttCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTcgac                  <  1:1389458/150‑1 (MQ=255)
                                                                                                                                     cACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTcgacga                <  2:1160791/150‑1 (MQ=255)
                                                                                                                                      aCCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgt                                                                                        <  1:679564/77‑1 (MQ=255)
                                                                                                                                      aCCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgt                                                                                        >  2:679564/1‑77 (MQ=255)
                                                                                                                                      aCCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGcc                                        >  2:1641798/1‑125 (MQ=255)
                                                                                                                                      aCCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGcc                                        <  1:1641798/125‑1 (MQ=255)
                                                                                                                                      aCCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGaa               <  2:1342993/150‑1 (MQ=255)
                                                                                                                                       ttcccAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAg              >  1:1135603/4‑150 (MQ=255)
                                                                                                                                         gCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCAc                               >  2:604240/1‑131 (MQ=255)
                                                                                                                                         gCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCAc                               <  1:604240/131‑1 (MQ=255)
                                                                                                                                          ccAGCAGCGCCTCGGCGTTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATg                                                           >  2:1404765/1‑102 (MQ=255)
                                                                                                                                          ccAGCAGCGCCTCGGCGTTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATg                                                           <  1:1404765/102‑1 (MQ=255)
                                                                                                                                              cAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGtggtgg                                                                             <  1:1181035/80‑1 (MQ=255)
                                                                                                                                              cAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGtggtgg                                                                             >  2:1181035/1‑80 (MQ=255)
                                                                                                                                                   ccGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGGCTg  <  1:815316/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACCGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGGCTG  >  NC_002947/3631979‑3632275

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