Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,031,922 T→A 6.3% L5M (TTG→ATG) ‡ PP_0889 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,031,9220TA 6.3% 350.7 / 12.8 114L5M (TTG→ATG) ‡PP_0889hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (50/55);  new base A (4/3);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

Only 100 of 115 total aligned reads displayed.

GCAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTC  >  NC_002947/1031773‑1032062
                                                                                                                                                     |                                                                                                                                             
gCAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCCCAACCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCt                                                                                                                                               <  1:1478140/150‑1 (MQ=255)
   atatCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGc                                                                                                                                           >  1:596230/1‑150 (MQ=255)
      tctcGGCCACCGCACGGGCGCTGGGCATGCACAGACGGCCGTGCCTTCGCAGTGTGTAGAAGCGGCCGGTTCGCTGTTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcg                                                                                                                                        <  1:435322/150‑1 (MQ=255)
        tcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAACCCGGGCTGGCCcgcg                                                                                                                                      >  2:354618/1‑150 (MQ=255)
          ggCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGaa                                                                                                                                    <  1:898345/150‑1 (MQ=255)
           gCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAg                                                                                                                                   <  1:75812/150‑1 (MQ=255)
             cACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAggg                                                                                                                                 <  2:772824/150‑1 (MQ=255)
               ccGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCt                                                                                                                               >  1:1706995/1‑150 (MQ=255)
                    cGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCTCGAAGGGCTGCaaa                                                                                                                          <  2:1734882/150‑1 (MQ=255)
                       gcgcTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCa                                                                                                                       <  2:1571211/150‑1 (MQ=255)
                           tGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTcc                                                                                                                   >  1:1209435/1‑150 (MQ=255)
                             ggCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGCTGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCa                                                                                                                 >  2:92782/1‑150 (MQ=255)
                              gCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCaa                                                                                                                <  2:596230/150‑1 (MQ=255)
                                   cACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                           >  2:1410177/1‑150 (MQ=255)
                                   cACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                           <  2:326107/150‑1 (MQ=255)
                                     ccGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGc                                                                                                         >  1:402756/1‑150 (MQ=255)
                                     ccGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgctcaca                                                                                                                                <  1:1104407/127‑7 (MQ=255)
                                       gacggacgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGa                                                                                                       >  1:532860/1‑150 (MQ=255)
                                              gTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGAt                                                                                                <  2:92068/150‑1 (MQ=255)
                                                 gCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTTTTGCAACTCTGTGGCCGCCTGGGTTGGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGaa                                                                                             >  1:717500/1‑150 (MQ=255)
                                                  cAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                            <  2:1209435/150‑1 (MQ=255)
                                                    gcgcAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGa                                                                                                 >  2:1700786/1‑143 (MQ=255)
                                                    gcgcAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGa                                                                                                 <  1:1700786/143‑1 (MQ=255)
                                                          gCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                            <  2:192751/142‑1 (MQ=255)
                                                          gCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGtt                                                                                    >  1:1246256/1‑150 (MQ=255)
                                                          gCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgc                                                                                                                                     <  2:1734028/101‑2 (MQ=255)
                                                          gCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgc                                                                                                                                     >  1:1734028/1‑100 (MQ=255)
                                                           cTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                   <  2:532860/150‑1 (MQ=255)
                                                            tGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTgg                                                                                  <  2:1510687/150‑1 (MQ=255)
                                                             gCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgct                                                                                                                                    <  1:943396/99‑3 (MQ=255)
                                                             gCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgct                                                                                                                                    >  2:943396/1‑97 (MQ=255)
                                                                      gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                        >  1:1860713/1‑150 (MQ=255)
                                                                        cGGTTCGGCGCTGATTCCTTATCGATGATGTTTGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCCTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgcgc                                                                      <  1:815823/150‑1 (MQ=255)
                                                                        cGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgcg                                                                       >  2:1080548/1‑149 (MQ=255)
                                                                        cGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgcg                                                                       <  1:1080548/149‑1 (MQ=255)
                                                                          gTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTg                                                                    >  1:684075/1‑150 (MQ=255)
                                                                             cGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGtt                                                                 <  2:402756/150‑1 (MQ=255)
                                                                              ggCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTAATCCTTGCGCTAGtgt                                                                >  2:1379197/1‑148 (MQ=255)
                                                                                cgcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCa                                                              <  2:1753769/150‑1 (MQ=255)
                                                                                 gcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATcc                                                                            <  1:344590/135‑1 (MQ=255)
                                                                                 gcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATcc                                                                            >  2:344590/1‑135 (MQ=255)
                                                                                 gcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                             <  1:102393/150‑1 (MQ=255)
                                                                                 gcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                             <  2:684075/150‑1 (MQ=255)
                                                                                  cTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCACTCTCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAgg                                                            <  1:726871/150‑1 (MQ=255)
                                                                                    gATTCCTTATCGATGAGGTTGGCCTGTTCGCGGGGTAATCCGCGCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTa                                                          >  1:1341050/1‑150 (MQ=255)
                                                                                     aTTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGGAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAg                                                         >  1:1794969/1‑150 (MQ=255)
                                                                                     aTTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATcc                                                                            >  2:1225022/1‑131 (MQ=255)
                                                                                     aTTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATcc                                                                            <  1:1225022/131‑1 (MQ=255)
                                                                                        ccTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgcgc                                                      <  1:1870304/150‑1 (MQ=255)
                                                                                         cTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGtt                                                                                    <  1:650033/119‑1 (MQ=255)
                                                                                         cTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGtt                                                                                    >  2:650033/1‑119 (MQ=255)
                                                                                         cTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTc                                                                                >  2:342906/1‑123 (MQ=255)
                                                                                         cTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTc                                                                                <  1:342906/123‑1 (MQ=255)
                                                                                              cGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGACTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGc                                                <  2:38955/150‑1 (MQ=255)
                                                                                               gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCt                                                                                                                               >  2:484093/1‑70 (MQ=255)
                                                                                               gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCt                                                                                                                               <  1:484093/70‑1 (MQ=255)
                                                                                               gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGACGGTACTGTTGGTCATCCTTGCGCTGGTTGCa                                                              <  1:1807848/135‑1 (MQ=255)
                                                                                               gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCtt                                                                          <  1:1043673/123‑1 (MQ=255)
                                                                                               gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCtt                                                                          >  2:1043673/1‑123 (MQ=255)
                                                                                               gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCa                                                              >  2:1807848/1‑135 (MQ=255)
                                                                                                atgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTg                                              >  1:1667052/1‑150 (MQ=255)
                                                                                                 tgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGtt                                                                 >  1:1561680/1‑130 (MQ=255)
                                                                                                 tgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGtt                                                                 <  2:1561680/130‑1 (MQ=255)
                                                                                                  gatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGCTGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGaa                                            <  1:354618/150‑1 (MQ=255)
                                                                                                     gTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATc                                                                                                              >  1:1898661/1‑81 (MQ=255)
                                                                                                      ttGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACt                                        <  2:530133/150‑1 (MQ=255)
                                                                                                         gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                     <  2:719012/150‑1 (MQ=255)
                                                                                                         gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                     <  2:1860713/150‑1 (MQ=255)
                                                                                                          ccTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTgg                                                                   <  1:110154/119‑1 (MQ=255)
                                                                                                          ccTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTgg                                                                   >  2:110154/1‑119 (MQ=255)
                                                                                                           cTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCt                                   >  2:805281/1‑150 (MQ=255)
                                                                                                             gTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGa                                 <  1:1410177/150‑1 (MQ=255)
                                                                                                                   gggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                             <  1:870398/116‑1 (MQ=255)
                                                                                                                   gggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                             >  2:870398/1‑116 (MQ=255)
                                                                                                                    ggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAg                          >  2:1638029/1‑150 (MQ=255)
                                                                                                                       gAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCCGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTgg                       >  1:242575/1‑150 (MQ=255)
                                                                                                                       gAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgcgc                                                                      <  2:1909921/103‑1 (MQ=255)
                                                                                                                       gAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgcgc                                                                      >  1:1909921/1‑103 (MQ=255)
                                                                                                                        agtCCGCTCCTGCAACTGTGAGGGCGGGCTTGCCCGCGAAGGACTTCAAAGCAGTGCCAATCAAGGCGAGCTCGATGAAGGTACTGTTTGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCGGACCCGCAGTGGc                      <  2:242575/148‑1 (MQ=255)
                                                                                                                          tCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAg                    >  2:421488/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAg                    >  1:999537/1‑150 (MQ=255)
                                                                                                                             gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                        <  2:1123783/95‑1 (MQ=255)
                                                                                                                             gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                        >  1:1123783/1‑95 (MQ=255)
                                                                                                                             gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                          >  2:929776/1‑125 (MQ=255)
                                                                                                                             gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                          <  1:929776/125‑1 (MQ=255)
                                                                                                                             gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAAc                 >  2:173256/1‑150 (MQ=255)
                                                                                                                               tCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGa               <  1:92782/150‑1 (MQ=255)
                                                                                                                                ccTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGAc              >  1:883188/1‑150 (MQ=255)
                                                                                                                                 cTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGAccc                              >  1:1707715/1‑133 (MQ=255)
                                                                                                                                 cTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGAccc                              <  2:1707715/133‑1 (MQ=255)
                                                                                                                                 cTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGAcc                               <  2:116585/132‑1 (MQ=255)
                                                                                                                                   gCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                          >  2:1118120/1‑119 (MQ=255)
                                                                                                                                   gCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                          <  1:1118120/119‑1 (MQ=255)
                                                                                                                                     aaCTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGc         <  2:717500/150‑1 (MQ=255)
                                                                                                                                     aaCTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGc         >  2:115645/1‑150 (MQ=255)
                                                                                                                                     aaCTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGc         >  1:1353607/1‑150 (MQ=255)
                                                                                                                                       cTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCa       >  1:1658803/1‑150 (MQ=255)
                                                                                                                                           gAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACt                                                                                       <  2:714297/66‑1 (MQ=255)
                                                                                                                                           gAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACt                                                                                       >  1:714297/1‑66 (MQ=255)
                                                                                                                                            aGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTc  <  1:805281/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GCAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTC  >  NC_002947/1031773‑1032062

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: