Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,914,180 C→T 10.9% intergenic (+16/+137) spuI → / ← argA glutamylpolyamine synthetase/amino acid N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,914,1800CT10.9% 249.5 / 17.1 92intergenic (+16/+137)spuI/argAglutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (41/41);  new base T (4/6);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGC  >  NC_002947/5914035‑5914307
                                                                                                                                                 |                                                                                                                               
cGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCaa                                                                                                                             >  1:2446666/1‑150 (MQ=255)
  agaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTGTAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAgg                                                                                                                           <  1:1983985/150‑1 (MQ=255)
   gaTCCTGAACAAGTACATCGATCCGAATTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAggg                                                                                                                          >  1:712216/1‑150 (MQ=255)
          aaCAAGTACATCGATCCGAAGTACATCGACCTCTTCGTCGCGTGCAAGGAAAGCGAGCTCGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCTTGAGGGGCGGCGtt                                                                                                                   <  2:1277641/150‑1 (MQ=39)
          aaCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGCCCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGtt                                                                                                                   <  1:1519355/150‑1 (MQ=255)
           aCAAGTACATCGATCCGAATTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttt                                                                                                                  <  2:712216/150‑1 (MQ=255)
           aCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttt                                                                                                                  <  1:380145/150‑1 (MQ=255)
             aaGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGTCGCCCTTGAGGGGCGGCGttttc                                                                                                                <  2:1181393/150‑2 (MQ=38)
              aGTACATCGATCCGCAATACATCGCCCTCCCCGTCGCGTGCAAGGCCAGCGAGCTGGAAGAGATCGAGCACTCGACCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttg                                                                                                               <  2:1084938/150‑1 (MQ=255)
                 aCATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                            <  2:764636/150‑1 (MQ=255)
                   aTCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTga                                                                                                          >  2:2123275/1‑150 (MQ=255)
                   aTCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTga                                                                                                          >  2:1740546/1‑150 (MQ=255)
                        tctgAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCg                                                                                                     <  1:2281655/147‑1 (MQ=255)
                        tCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCg                                                                                                     >  1:1420289/1‑150 (MQ=255)
                        tCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCg                                                                                                     >  2:932430/1‑150 (MQ=255)
                               tACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGC                                                                                              <  1:2109539‑M1/150‑1 (MQ=255)
                                aCATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCC                                                                                             <  2:1574291‑M1/150‑1 (MQ=255)
                                  aTCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGgt                                                                                                   <  1:2200958/142‑1 (MQ=255)
                                  aTCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGgt                                                                                                   >  2:2200958/1‑142 (MQ=255)
                                   tCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcg                                                                                                                        <  1:738725/120‑1 (MQ=255)
                                   tCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcg                                                                                                                        >  2:738725/1‑120 (MQ=255)
                                                    tGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGC                                                                         <  2:1907401‑M1/150‑1 (MQ=255)
                                                    tGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGC                                                                         >  2:668721‑M1/1‑150 (MQ=255)
                                                    tGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGGGCGGCGtt                                                                                                                   >  2:1009665/1‑108 (MQ=38)
                                                    tGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGGGCGGCGtt                                                                                                                   <  1:1009665/108‑1 (MQ=38)
                                                       aaGGGAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCC                                                                      >  1:636596‑M1/1‑150 (MQ=255)
                                                       aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGCgaa                                                                      <  1:301834‑M1/150‑4 (MQ=255)
                                                        aGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATAACGTGTCCTCGGTATCGCCGCCCCTCAAGGGGGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCC                                                                     <  1:2123275‑M1/150‑1 (MQ=255)
                                                          gAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCaa                                                                                                                             <  2:588645/92‑1 (MQ=255)
                                                          gAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCaa                                                                                                                             >  1:588645/1‑92 (MQ=255)
                                                          gAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGc                                                                   >  2:930097/1‑150 (MQ=255)
                                                               cGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGt                                                              <  2:1420289/150‑1 (MQ=255)
                                                                    tGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                           <  2:1073478‑M1/132‑1 (MQ=255)
                                                                    tGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                           >  1:1073478‑M1/1‑132 (MQ=255)
                                                                     ggAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcac                                                        >  2:312289‑M1/1‑137 (MQ=255)
                                                                      gAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcaca                                                       <  2:1001593‑M1/150‑15 (MQ=255)
                                                                      gAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTCTGATGTTAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcaca                                                       >  1:24968‑M1/1‑136 (MQ=255)
                                                                        agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGc                                                                   <  2:1290190/136‑1 (MQ=255)
                                                                        agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGc                                                                   >  1:1290190/1‑136 (MQ=255)
                                                                            ttCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCt                                                 <  2:1563681/150‑1 (MQ=255)
                                                                                 gCACTCGATTTCCGACCTCGAGTACAACTGGTACCTGCATACCGGGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGc                                            <  2:669861/150‑1 (MQ=255)
                                                                                 gCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGc                                            <  2:636596/150‑1 (MQ=255)
                                                                                  cACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagcctt                                           >  2:1014011‑M1/1‑124 (MQ=255)
                                                                                     tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcga                                        >  2:548422‑M1/1‑121 (MQ=255)
                                                                                       gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                      >  1:1411251/1‑150 (MQ=255)
                                                                                          ctcCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACtt                                   >  2:198920/1‑150 (MQ=255)
                                                                                              gACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCt                               <  2:2446666/150‑1 (MQ=255)
                                                                                               aCCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTgag                                                                                                         >  2:1284740/1‑75 (MQ=255)
                                                                                               aCCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTgag                                                                                                         <  1:1284740/75‑1 (MQ=255)
                                                                                                ccTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgc                             <  1:312289‑M1/150‑41 (MQ=255)
                                                                                                ccTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgc                             <  1:1014011‑M1/150‑41 (MQ=255)
                                                                                                 cTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTg                            >  1:932430‑M2/1‑149 (MQ=255)
                                                                                                  tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCTTGAGGGGCGGCGttttc                                                                                                                >  2:564172/1‑64 (MQ=21)
                                                                                                  tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCTTGAGGGGCGGCGttttc                                                                                                                <  1:564172/65‑2 (MQ=21)
                                                                                                  tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                     >  1:2180289/1‑124 (MQ=255)
                                                                                                  tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                     <  2:2180289/124‑1 (MQ=255)
                                                                                                  tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCt                                                 >  2:722046/1‑128 (MQ=255)
                                                                                                  tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCt                                                 <  1:722046/128‑1 (MQ=255)
                                                                                                  tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                           >  1:1799836‑M2/1‑148 (MQ=255)
                                                                                                   cGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTggg                          >  2:1444593‑M2/1‑147 (MQ=255)
                                                                                                    gAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCTCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTggga                         <  2:1225607‑M2/150‑5 (MQ=255)
                                                                                                    gAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTggga                         >  1:668721‑M2/1‑146 (MQ=255)
                                                                                                      gTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagc                       >  2:2247332‑M2/1‑144 (MQ=255)
                                                                                                       tACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGTCACTCCAGTGCTGCTTTTTTGTTTGAGACCGGTTACGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcg                      <  1:1740546‑M1/150‑48 (MQ=255)
                                                                                                       tACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcggga                         >  2:842775‑M1/1‑103 (MQ=255)
                                                                                                       tACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcggga                         <  1:842775‑M1/147‑45 (MQ=255)
                                                                                                       tACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcg                      >  2:1293909‑M2/1‑143 (MQ=255)
                                                                                                        aCAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgg                     <  1:548422‑M1/150‑49 (MQ=255)
                                                                                                        aCAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgg                     >  1:930097‑M2/1‑142 (MQ=255)
                                                                                                          aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                             >  2:2066661/1‑140 (MQ=255)
                                                                                                          aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                             <  1:2066661/140‑1 (MQ=255)
                                                                                                          aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgg                     <  2:2386332‑M1/148‑49 (MQ=255)
                                                                                                          aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgg                     >  1:2386332‑M1/1‑100 (MQ=255)
                                                                                                           cacggtacCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGG                                                                                   <  1:1063766‑M1/82‑1 (MQ=255)
                                                                                                           cacggtacCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGG                                                                                   >  2:1063766‑M1/4‑85 (MQ=255)
                                                                                                           aCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcg                  <  2:1770803‑M2/150‑12 (MQ=255)
                                                                                                            cTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgt                 >  1:198920‑M2/1‑138 (MQ=255)
                                                                                                            cTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCATCTTCGCGGGCACGCACGATCCGGCGTGATGATTTCGTATTCTGAAGCCTGTGACTTACCTGTgggagcgggcgt                 >  1:1225607‑M2/1‑138 (MQ=255)
                                                                                                                tACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgccc             >  1:1770803‑M2/1‑134 (MQ=255)
                                                                                                                   cTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgg                     >  1:69135‑M2/1‑131 (MQ=255)
                                                                                                                   cTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgg                     <  2:69135‑M2/139‑9 (MQ=255)
                                                                                                                           cGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcgaagagggc  <  2:2385594‑M2/150‑28 (MQ=255)
                                                                                                                           cGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTACCTGTgggagcgggcgt                 <  2:538301‑M2/134‑13 (MQ=255)
                                                                                                                           cGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTACCTGTgggagcgggcgt                 >  1:538301‑M2/1‑122 (MQ=255)
                                                                                                                                aaaCAAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                      <  1:1886766/109‑1 (MQ=255)
                                                                                                                                aaaCAAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                      >  2:1886766/1‑109 (MQ=255)
                                                                                                                                 aaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcggg                    <  1:2058150‑M2/126‑10 (MQ=255)
                                                                                                                                 aaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcggg                    >  2:2058150‑M2/1‑117 (MQ=255)
                                                                                                                                 aaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagc                       <  1:936940‑M2/123‑7 (MQ=255)
                                                                                                                                 aaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagc                       >  2:936940‑M2/1‑117 (MQ=255)
                                                                                                                                    gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                     >  1:1917488/1‑90 (MQ=255)
                                                                                                                                    gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                     <  2:1917488/90‑1 (MQ=255)
                                                                                                                                                 cTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                     <  2:2450963/77‑1 (MQ=255)
                                                                                                                                                 cTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                     >  1:2450963/1‑77 (MQ=255)
                                                                                                                                                 |                                                                                                                               
CGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGC  >  NC_002947/5914035‑5914307

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