Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,090,935 Δ1 bp 5.2% coding (50/1020 nt) astA‑II ← arginine N‑succinyltransferase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,090,9350C. 5.2% 367.9 / 13.9 116coding (50/1020 nt)astA‑IIarginine N‑succinyltransferase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (59/51);  new base . (3/3);  total (62/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 119 total aligned reads displayed.

GTCTTCGAGCACGAAAAAGTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGCCGACCCGGATCA  >  NC_002947/5090790‑5091083
                                                                                                                                                  |                                                                                                                                                    
   ttCGAGCACGAAAAAGTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCAAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCAc                                                                                                                                               <  1:2202814/150‑1 (MQ=255)
   ttCGAGCACGAAAAAGTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGACGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCAc                                                                                                                                               >  2:12749/1‑150 (MQ=255)
      gAGCACGAAAAAGTAACTCTCTTCAACGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCAcgtc                                                                                                                                            <  1:2776338/150‑1 (MQ=255)
         cACGAAAAAGTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCCCGTCGtt                                                                                                                                         >  2:748084/1‑150 (MQ=255)
         cACGAAAAAGTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGtt                                                                                                                                         <  1:1324248/150‑1 (MQ=255)
            gAAAAAGTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCTCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAg                                                                                                                                      >  2:2166821/1‑150 (MQ=255)
            gAAAAAGTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAg                                                                                                                                      <  1:2272008/150‑1 (MQ=255)
            gAAAAAGTAACTCTCTCCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGtt                                                                                                                                         <  2:73851/147‑1 (MQ=255)
            gAAAAAGTAACTCTCTCCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGtt                                                                                                                                         >  1:73851/1‑147 (MQ=255)
              aaaaGTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGAt                                                                                                                                    <  2:2157571/150‑1 (MQ=255)
                aaGTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCa                                                                                                                                  <  1:1320559/150‑1 (MQ=255)
                  gTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCTCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGc                                                                                                                                <  1:2166821/150‑1 (MQ=255)
                   tAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGcc                                                                                                                               <  2:2504544/150‑1 (MQ=255)
                           ttCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTgcgc                                                                                                                       >  1:1479937/1‑150 (MQ=255)
                           ttCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTgcgc                                                                                                                       >  2:2558232/1‑150 (MQ=255)
                            tCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGcgcg                                                                                                                      <  1:1944242/150‑1 (MQ=255)
                             cACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTCCAGCCATGCGCTGCAGGTCGTTCAGATCAGCCATTTGCGCgg                                                                                                                     >  1:814709/1‑150 (MQ=255)
                                  ttGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGgcgc                                                                                                                >  2:1167002/1‑150 (MQ=255)
                                      aGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCa                                                                                                            >  1:2433205/1‑150 (MQ=255)
                                        cTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCAcc                                                                                                          <  2:735661/150‑1 (MQ=255)
                                          gACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAg                                                                                                        >  1:1244752/1‑150 (MQ=255)
                                            ccTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCa                                                                                                      >  2:820086/1‑150 (MQ=255)
                                             cTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCAt                                                                                                     >  1:1880482/1‑150 (MQ=255)
                                               cggcggCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCatgctgg                                                                                                   <  1:2558232/150‑7 (MQ=255)
                                                 gcggcgAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGgtg                                                                                                 <  1:1103041/150‑1 (MQ=255)
                                                      gAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATTCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACt                                                                                            >  1:935585/1‑150 (MQ=255)
                                                      gAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACt                                                                                            >  1:850244/1‑150 (MQ=255)
                                                          gATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCtt                                                                                        >  1:862430/1‑150 (MQ=255)
                                                           aTGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCttt                                                                                       <  2:140187/150‑1 (MQ=255)
                                                           aTGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCttt                                                                                       <  2:862430/150‑1 (MQ=255)
                                                            tgcctcAGAGGCGGCAATCTTGTCGCCCAGACGGGCAGTGGCGGCCGGGAGCGAGGGGACACCAATGGGGCTGGCTGCAGCCATGCGCTGCACGTTGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTg                                                                                      <  2:814709/147‑1 (MQ=255)
                                                            tGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTg                                                                                      <  2:167557/150‑1 (MQ=255)
                                                             gTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTgg                                                                                     >  2:752002/1‑150 (MQ=255)
                                                              tctcAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCttt                                                                                       >  2:1543897/1‑147 (MQ=255)
                                                              tctcAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCttt                                                                                       <  1:1543897/147‑1 (MQ=255)
                                                                  agagGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGAt                                                                                                                                    <  1:1387673/98‑1 (MQ=255)
                                                                  agagGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGAt                                                                                                                                    >  2:1387673/1‑98 (MQ=255)
                                                                   gagGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACg                                                                               <  2:1479937/150‑1 (MQ=255)
                                                                    agGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACgg                                                                              <  1:1167002/150‑1 (MQ=255)
                                                                      gcggcAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGc                                                                            <  1:820086/150‑1 (MQ=255)
                                                                           aaTCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGGCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCtt                                                                                        >  1:371928/1‑133 (MQ=255)
                                                                             tCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCACTGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCAc                                                                     >  2:430212/1‑150 (MQ=255)
                                                                               ttGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATgcg                                                                                                                                                     >  1:1156168/1‑68 (MQ=255)
                                                                               ttGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATgcg                                                                                                                                                     <  2:1156168/68‑1 (MQ=255)
                                                                                tGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTc                                                                  <  2:850244/150‑1 (MQ=255)
                                                                                      ccAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACg                                                            >  1:383810/1‑150 (MQ=255)
                                                                                      ccAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACg                                                            >  2:1793146/1‑150 (MQ=255)
                                                                                            gAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGgtg                                                                                                 <  2:185922/107‑1 (MQ=255)
                                                                                            gAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGaaaaa                                                      >  1:18340/1‑150 (MQ=255)
                                                                                              gCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACggg                                                                             >  1:778968/1‑125 (MQ=255)
                                                                                               cAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGcgcg                                                                                                                      >  2:763840/1‑83 (MQ=255)
                                                                                               cAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGcgcg                                                                                                                      <  1:763840/83‑1 (MQ=255)
                                                                                               cAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAAc                                                     <  1:2312675/148‑1 (MQ=255)
                                                                                               cAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAAc                                                     >  2:2312675/1‑148 (MQ=255)
                                                                                                aGTGTCGTTCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGgt                                                                                                  <  1:2045008/102‑1 (MQ=255)
                                                                                                  tgtCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCg                                                                 >  1:4401/1‑133 (MQ=255)
                                                                                                  tgtCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCg                                                                 <  2:4401/133‑1 (MQ=255)
                                                                                                      gtcCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTc                                                                                           <  2:2780358/103‑1 (MQ=255)
                                                                                                      gtcCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTc                                                                                           >  1:2780358/1‑103 (MQ=255)
                                                                                                       tcCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGGGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGcaggcacc                                           >  1:2359717/1‑150 (MQ=255)
                                                                                                         cGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGcaggcacccc                                         >  2:2996168/1‑149 (MQ=255)
                                                                                                         cGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGcaggcaccca                                         <  1:752002/150‑1 (MQ=255)
                                                                                                         cGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGcaggcaccca                                         >  2:2570682/1‑150 (MQ=255)
                                                                                                          ggCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGgtg                                                                                                 >  1:2016433/1‑93 (MQ=255)
                                                                                                          ggCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGgtg                                                                                                 <  2:2016433/93‑1 (MQ=255)
                                                                                                           gCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGcaggcacccagg                                       >  2:538174/1‑150 (MQ=255)
                                                                                                             aGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGcaggcacccaggca                                     >  2:1084236/1‑150 (MQ=255)
                                                                                                              gCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGcaggcacccaggcac                                    <  2:2433205/150‑1 (MQ=255)
                                                                                                                 aGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTg                                 <  1:1294821/150‑1 (MQ=255)
                                                                                                                     gACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACt                             >  1:2192991/1‑150 (MQ=255)
                                                                                                                      acacCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTg                            >  1:1510512/1‑150 (MQ=255)
                                                                                                                        acCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCAt                                                                                                     >  2:2744254/1‑75 (MQ=255)
                                                                                                                        acCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCAt                                                                                                     <  1:2744254/75‑1 (MQ=255)
                                                                                                                        acCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGa                          >  1:388233/1‑150 (MQ=255)
                                                                                                                         tgatctgaacgacgTGCAGCGCATGGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACgg                                                           <  2:334924‑M1/102‑1 (MQ=255)
                                                                                                                         tgatctgaacgacgTGCAGCGCATGGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACgg                                                           >  1:334924‑M1/15‑116 (MQ=255)
                                                                                                                          cAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAAt                        >  1:2524078/1‑150 (MQ=255)
                                                                                                                           aaTGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGcgcg                                                                                                                      >  2:1379575/1‑55 (MQ=255)
                                                                                                                           aaTGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGcgcg                                                                                                                      <  1:1379575/55‑1 (MQ=255)
                                                                                                                           aaTGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACggc                                                  <  2:1121261/123‑1 (MQ=255)
                                                                                                                           aaTGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACggc                                                  >  1:1121261/1‑123 (MQ=255)
                                                                                                                           aaTGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATc                       <  2:2168917/150‑1 (MQ=255)
                                                                                                                            aTGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGgc                                                                                                                  >  1:1648645/1‑58 (MQ=255)
                                                                                                                              ggggCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGGCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCggg                    <  1:2570682/150‑1 (MQ=255)
                                                                                                                               gggCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATgg                                                                                                   >  2:2995683/1‑70 (MQ=255)
                                                                                                                               gggCTGTCGTCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATgg                                                                                                   <  1:2995683/70‑1 (MQ=255)
                                                                                                                               gaacgacgTGCAGCGCATGGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCAt                                                                                                     >  2:2008253‑M1/9‑68 (MQ=255)
                                                                                                                               gaacgacgTGCAGCGCATGGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCAt                                                                                                     <  1:2008253‑M1/60‑1 (MQ=255)
                                                                                                                                  cTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATc                <  1:1084236/150‑1 (MQ=255)
                                                                                                                                    gTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGc              >  1:766521/1‑150 (MQ=255)
                                                                                                                                     tCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGcc             >  1:719258/1‑150 (MQ=255)
                                                                                                                                     tCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGcc             >  2:838178/1‑150 (MQ=255)
                                                                                                                                     tCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCAATCTTGGCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGcc             >  1:2806912/1‑150 (MQ=255)
                                                                                                                                     cgTGCAGCGCATGGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACAAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCgg                     >  2:2670737‑M1/3‑142 (MQ=255)
                                                                                                                                     cgTGCAGCGCATGGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACAAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCgg                     <  1:2670737‑M1/140‑1 (MQ=255)
                                                                                                                                       TGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGCCGa           <  2:506166‑M1/150‑1 (MQ=255)
                                                                                                                                       TGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGCCGa           >  2:2524078‑M1/1‑150 (MQ=255)
                                                                                                                                              cATGCGCTGAACGTCGTTCAGATCAGTCATTTGCGCGGGGTTCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGCCGACCCGGa     <  2:2192991/150‑1 (MQ=255)
                                                                                                                                                tGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGCCGACCCGGATc   <  1:2601470/150‑1 (MQ=255)
                                                                                                                                                 gcgcTGCACGTCGTTCAGATCAGCCATTTGGGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGCCGACCCGGATCa  <  2:388233/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTCTTCGAGCACGAAAAAGTAACTCTCTTCACCGTTGAAGCTGACCTCGGCGGCGAAGGATGTCTCAGAGGCGGCAATCTTGTCGCCCAGACGAGCAGTGTCGTCCGGCAGCGAGGTGACACCAATGGGGCTGTCTGCAGCCATGCGCTGCACGTCGTTCAGATCAGCCATTTGCGCGGGGCGCATCACCAGCATGGTGTCACTCCTTTGGAATACGGGCCATTCACGTCGGCACGGAAAAACGGCAGGCACCCAGGCACCTGCACTGGAATCGGGTATCGCCGACCCGGATCA  >  NC_002947/5090790‑5091083

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: