Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,914,212 T→C 5.7% intergenic (+48/+105) 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,2120TC 5.7% 213.7 / 2.2 70intergenic (+48/+105)spuI/argAglutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base C (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

AAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAG  >  NC_002947/5914063‑5914356
                                                                                                                                                     |                                                                                                                                                
aaGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGgtgt                                                                                                                                                  <  1:1723316/150‑1 (MQ=255)
  gTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCG                                                                                                                                                >  1:341059‑M1/1‑150 (MQ=255)
       tCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTC                                                                                                                                           >  1:1746457‑M1/1‑150 (MQ=255)
       tCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTC                                                                                                                                           >  2:210632‑M1/1‑150 (MQ=255)
              cTTCGCCGCGTGAAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGG                                                                                                                                    <  1:1020086‑M1/150‑1 (MQ=255)
                    cgcgTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTGCAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGC                                                                                                                              <  2:1514091‑M1/150‑1 (MQ=255)
                    cgcgTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGC                                                                                                                              <  1:1048459‑M1/150‑1 (MQ=255)
                    cgcgTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGC                                                                                                                              <  1:1865973‑M1/150‑1 (MQ=255)
                     gcgTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCC                                                                                                                             <  2:341059‑M1/150‑1 (MQ=255)
                        tGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGC                                                                                                                          <  2:1332381‑M1/150‑1 (MQ=255)
                         gCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTT                                                                                                                                         <  2:1984562‑M1/134‑1 (MQ=255)
                         gCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTT                                                                                                                                         >  1:1984562‑M1/1‑134 (MQ=255)
                          cAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC                                                                                                                        >  1:1857355‑M1/1‑150 (MQ=255)
                          cAAGGAAAGCAAGCTCGAAGAGTTCGAGCACTCGATCTCCGCCCCCGAGTACAACTGGTCCCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC                                                                                                                        <  1:502284‑M1/150‑1 (MQ=255)
                           aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCC                                                                                                                       >  2:409278‑M1/1‑150 (MQ=255)
                              gAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGc                                                                                                                    >  2:643836/1‑150 (MQ=255)
                                 aGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                                                                 <  2:2188248/150‑1 (MQ=255)
                                 aGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                                                                 <  1:1478268/150‑1 (MQ=255)
                                  gCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                                            <  1:894922‑M1/138‑1 (MQ=255)
                                  gCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                                            >  2:894922‑M1/1‑138 (MQ=255)
                                                   gAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGaa                                                                                               <  1:1947027/150‑1 (MQ=255)
                                                     gCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGc                                                                                             <  1:643836/150‑1 (MQ=255)
                                                       aCTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCt                                                                                           >  1:999847/1‑150 (MQ=255)
                                                         tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                                         <  2:1746457/150‑1 (MQ=255)
                                                         tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                                         <  2:1361491/150‑1 (MQ=255)
                                                         tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                                         >  2:1274299/1‑150 (MQ=255)
                                                           gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                                                       <  1:409278/150‑1 (MQ=255)
                                                           gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                                                       <  2:1857355/150‑1 (MQ=255)
                                                            aTCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGTTAATTTCGAAGGCTGAAGCCTGTGAc                                                                                      <  1:428638/150‑1 (MQ=255)
                                                               tcCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattac                                                                                   >  2:2292032‑M1/1‑115 (MQ=255)
                                                                 cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAcc                                                                                 >  1:2096973/1‑150 (MQ=255)
                                                                 cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAcc                                                                                 >  2:1448028/1‑150 (MQ=255)
                                                                      tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                                                      >  1:1091104/1‑124 (MQ=255)
                                                                      tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                                                      <  2:1091104/124‑1 (MQ=255)
                                                                      tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                                                            >  1:1274299‑M2/1‑148 (MQ=255)
                                                                          gTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAgg                                                                                                    <  2:1392924/122‑1 (MQ=255)
                                                                          gTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAgg                                                                                                    >  1:1392924/1‑122 (MQ=255)
                                                                             cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcggg                                                                     <  2:2238315‑M2/150‑10 (MQ=255)
                                                                              aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGTGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggc                                                                    >  1:67893‑M2/1‑140 (MQ=255)
                                                                              aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCgcagggtagttcacaggcttcagccttcgaaattaccctgcgg                                                                           >  1:1409745‑M1/1‑100 (MQ=255)
                                                                              aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCgcagggtagttcacaggcttcagccttcgaaattaccctgcgg                                                                           <  2:1409745‑M1/143‑44 (MQ=255)
                                                                              aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggc                                                                    <  2:41425‑M1/150‑51 (MQ=255)
                                                                              aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggc                                                                    <  1:2292032‑M1/150‑51 (MQ=255)
                                                                                cTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgt                                                                  >  2:999847‑M2/1‑138 (MQ=255)
                                                                                      ccTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                                                            >  2:1979363‑M2/1‑132 (MQ=255)
                                                                                      ccTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                                                            <  1:1979363‑M2/134‑3 (MQ=255)
                                                                                       cTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                                         >  1:1224484/1‑120 (MQ=255)
                                                                                       cTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                                         <  2:1224484/120‑1 (MQ=255)
                                                                                          cATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcgaag                                                        <  1:1365381‑M2/150‑23 (MQ=255)
                                                                                           aTACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGCAGGCTAAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcgaaga                                                       >  2:1054816‑M2/1‑127 (MQ=255)
                                                                                           aTACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcgaaga                                                       >  2:2096973‑M2/1‑127 (MQ=255)
                                                                                                  tcaaacaAAAACGCCGCCCTTGAGAGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcaca                                                                                                        <  2:1052441‑M1/92‑15 (MQ=255)
                                                                                                  tcaaacaAAAACGCCGCCCTTGAGAGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcaca                                                                                                        >  1:1052441‑M1/3‑80 (MQ=255)
                                                                                                       cGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAgg                                           <  2:1365381/150‑1 (MQ=255)
                                                                                                       cGAAAACGCCCCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAgg                                           <  1:2238315/150‑1 (MQ=255)
                                                                                                        gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggcgtgcccgcgaagagg                                                     <  1:167823‑M1/139‑66 (MQ=255)
                                                                                                        gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggcgtgcccgcgaagagg                                                     >  2:167823‑M1/1‑74 (MQ=255)
                                                                                                           aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTa                                       >  1:1509754/1‑150 (MQ=255)
                                                                                                           aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTa                                       >  1:1278678/1‑150 (MQ=255)
                                                                                                              gccgccCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcgaaga                                                       <  2:1070901‑M2/131‑24 (MQ=255)
                                                                                                              gccgccCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcgaaga                                                       >  1:1070901‑M2/1‑108 (MQ=255)
                                                                                                                     cTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgc                                                                >  1:1421555‑M2/1‑101 (MQ=255)
                                                                                                                     cTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgc                                                                <  2:1421555‑M2/115‑15 (MQ=255)
                                                                                                                            gcggcgTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcg                                                                       <  2:1889363‑M1/101‑48 (MQ=255)
                                                                                                                            gcggcgTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcg                                                                       >  1:1889363‑M1/1‑54 (MQ=255)
                                                                                                                                  tttttGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACTCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCg                >  2:1047234/1‑150 (MQ=255)
                                                                                                                                  tttttGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGc                                                                                             <  1:544823/73‑1 (MQ=255)
                                                                                                                                  tttttGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGc                                                                                             >  2:544823/1‑73 (MQ=255)
                                                                                                                                              aCCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgg                                                                            >  1:1727236‑M1/1‑36 (MQ=255)
                                                                                                                                              aCCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgg                                                                            <  2:1727236‑M1/78‑43 (MQ=255)
                                                                                                                                                cGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAg  >  2:1206259/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAG  >  NC_002947/5914063‑5914356

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