Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,428 A→C 6.3% intergenic (‑330/+268) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,153,4280AC 6.3% 228.8 / 7.2 80intergenic (‑330/+268)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (29/46);  new base C (3/2);  total (32/48)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

CACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGG  >  NC_000913/4153281‑4153576
                                                                                                                                                   |                                                                                                                                                    
cACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGa                                                                                                                                                      <  1:388181/148‑1 (MQ=255)
 acgacgCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCTCTGGAGTTTGGCAACAGTGCCTGATGAGGCGCGAgc                                                                                                                                                    >  7:178963/1‑149 (MQ=255)
  cgacgCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcg                                                                                                                                                   >  4:214478/1‑149 (MQ=255)
   gacgCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                                                                                                                                                  <  4:192730/149‑1 (MQ=255)
   gacgCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                                                                                                                                                  >  2:19029/1‑149 (MQ=255)
   gacgCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                                                                                                                                                  >  4:227690/1‑149 (MQ=255)
       cTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCCCATACGACAATATTGCGGCACACGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                              >  1:179597/1‑149 (MQ=25)
        tttAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                                                             >  2:175261/1‑149 (MQ=255)
        tttAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                                                             >  1:59038/1‑149 (MQ=255)
           aaaTAAGCGTAACTTATGGAAATTTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                                           <  4:24251/148‑1 (MQ=255)
             aTAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                                        >  7:66584/1‑149 (MQ=255)
                aGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                                                                     <  7:247875/149‑1 (MQ=255)
                aGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                                                                     <  3:126750/149‑1 (MQ=255)
                   gTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAg                                                                                                                                  <  2:76957/149‑1 (MQ=255)
                    tAACTTATGGAAATGTTAAAAAATCTCCCCAAGTAACACCACAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGt                                                                                                                                 <  6:346694/149‑1 (MQ=255)
                    tAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGt                                                                                                                                 <  3:104769/149‑1 (MQ=255)
                       cTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                                                              <  3:403/149‑1 (MQ=255)
                       cTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                                                              <  1:19029/149‑1 (MQ=255)
                       cTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                                                              <  4:332041/149‑1 (MQ=255)
                          aTGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGc                                                                                                                            <  1:175261/148‑1 (MQ=255)
                           tGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCa                                                                                                                           <  8:44451/148‑1 (MQ=255)
                             gAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATc                                                                                                                         <  3:344640/148‑1 (MQ=255)
                                aTGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGttt                                                                                                                     <  5:99886/149‑1 (MQ=255)
                                aTGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGttt                                                                                                                     <  4:256856/149‑1 (MQ=255)
                                   ttAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGGTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGa                                                                                                                   <  3:264107/148‑1 (MQ=255)
                                       aaaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCAc                                                                                                              <  1:279142/149‑1 (MQ=255)
                                        aaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACt                                                                                                             <  2:186483/149‑1 (MQ=255)
                                          aTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgt                                                                                                           >  6:351674/1‑149 (MQ=255)
                                           tCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtc                                                                                                          >  1:270745/1‑149 (MQ=255)
                                              ccccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtaggt                                                                                                       >  1:273486/1‑149 (MQ=255)
                                               cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggtc                                                                                                      <  7:28712/149‑1 (MQ=255)
                                               cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTAGGCAACAGTGACGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggtc                                                                                                      >  1:275925/1‑149 (MQ=255)
                                                  aaGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGa                                                                                                   <  5:351674/149‑1 (MQ=255)
                                                   aGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGAt                                                                                                  <  3:214478/149‑1 (MQ=255)
                                                    gTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATa                                                                                                 <  8:119660/149‑1 (MQ=255)
                                                     tAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGAt                                                                                                  >  3:180795/1‑147 (MQ=255)
                                                       acacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCAGCGTTTGAGTCACTGTCGGTCGGATAAg                                                                                               <  1:165506/148‑1 (MQ=255)
                                                         acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTAAAGTTGGGCATCCTTTGCGTCACTCTCCGTCGGATAAGaag                                                                                            >  7:337150/1‑147 (MQ=14)
                                                         acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTAGTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATg                                                                                            >  7:350206/1‑149 (MQ=255)
                                                           cAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                                                                                                                                                 <  3:32320/94‑1 (MQ=40)
                                                           cAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                                                                                                                                                 >  4:32320/1‑94 (MQ=40)
                                                           cAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCGTGGTTTGAGTCACTGTCGGTCGGATAAGATgcg                                                                                          >  7:177902/1‑149 (MQ=255)
                                                              aGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCaa                                                                                       <  3:312738/149‑1 (MQ=39)
                                                                 tgtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGt                                                                                     >  3:134124/1‑148 (MQ=38)
                                                                   tAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCAc                                                                                                              >  4:310353/1‑121 (MQ=34)
                                                                   tAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCAc                                                                                                              <  3:310353/121‑1 (MQ=34)
                                                                    aGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATc                                                                                  <  7:321805/148‑1 (MQ=38)
                                                                        cGGATAAGATGCTCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATc                                                                             <  1:102771/149‑1 (MQ=37)
                                                                            tAAGATGCGCAAGTATCGAATCCGACTTTATTGCGGCACTGGAGTTGTCCAACATTGCCGGATGCTTCGCGAGCGCCTTATCCGGACTACAGTTGGGCATCTTTTGAGTCTCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGAc                                                                         <  2:275925/149‑1 (MQ=255)
                                                                             aaGCTGCGCAAGTCTCGCATCCGACATTATTGCGGCACGGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTGTGAGTCACTGTCGGGCGGATAAGATGCGCAAGTATCGCAt                                                                              <  7:366533/143‑1 (MQ=11)
                                                                             aaGATGCGCAAGTATCGCATCCGACATTAATGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAGCCGGCATACAGTTGGGCATCGTGGGAGTCGCTGTCGGTCGGATAAGATGCGCAAGTATCGCag                                                                              >  8:366533/1‑142 (MQ=11)
                                                                               gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACa                                                                        >  7:33425/1‑147 (MQ=37)
                                                                               gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACa                                                                        <  8:33425/147‑1 (MQ=37)
                                                                                   cgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACattatt                                                                   <  6:153593/148‑1 (MQ=37)
                                                                                      aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAg                                                                                               >  8:323630/1‑117 (MQ=35)
                                                                                      aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAg                                                                                               <  7:323630/117‑1 (MQ=35)
                                                                                           tagaaTCCGACATTATTGCGGCACTGGAGTTTGTCAACAGTGCCGTATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTg                                                          <  1:235772/145‑1 (MQ=255)
                                                                                           tCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                                           >  7:335274/1‑116 (MQ=35)
                                                                                           ccgatctcGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTCGCCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTg                                                          <  5:288876/142‑1 (MQ=255)
                                                                                             gCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGGTGCGGCGCGAGCGGCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                        <  7:14754/149‑1 (MQ=255)
                                                                                                   gACATTATTGCGGCACTGGAGATTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATAAGGCGTACAGTTGGGCATCGGATGAGAAAATGTCGTTAGGATAGGATGCGCGGGGGTCGCCTGCGACAGTGTTGGGGCGCTGGAgatcgg                                                  >  6:288876/1‑144 (MQ=2)
                                                                                                         attGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGCTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACa                                             <  4:180795/148‑1 (MQ=37)
                                                                                                             cGGCACTGGAGTTTGGCAACTGTGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGGCGGATAAGATGCGCAAGTATCGCATCCGACTTTATTGCGGCACTGGAGTTTGGCAACAGTGcc                                        <  2:385157/149‑1 (MQ=255)
                                                                                                                   tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCAt                                                                              >  7:241528/1‑105 (MQ=17)
                                                                                                                     gAGTTTGGCAACAGTTCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCTCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcg                                 <  8:157176/148‑1 (MQ=255)
                                                                                                                      aGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcggc                               <  4:134124/149‑1 (MQ=39)
                                                                                                                           ggCTACAGTGCAGGATGCGGCTCGAGCTCCTTATCATTCATACATTTGGGAATCTTTTGATTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGa                           <  4:109672/148‑1 (MQ=14)
                                                                                                                            gCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgc                         <  8:110623/149‑1 (MQ=39)
                                                                                                                              ctgtaggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                       <  4:327867/143‑1 (MQ=14)
                                                                                                                              aaCAGTGGAGTATGCGGCTCGATCGCCTTATCAGGCCTACAGTTGGGCATAGTTTGAGTAACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                       <  8:238498/149‑1 (MQ=14)
                                                                                                                              aaCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGt                                                      >  1:351282/1‑118 (MQ=39)
                                                                                                                               aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                        >  1:326107/1‑115 (MQ=39)
                                                                                                                               aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                      >  7:301419/1‑149 (MQ=39)
                                                                                                                               aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                      <  8:66584/149‑1 (MQ=39)
                                                                                                                               aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGAACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGTGcc                      >  8:83259/1‑149 (MQ=255)
                                                                                                                                    gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                       >  2:372727/1‑127 (MQ=14)
                                                                                                                                    gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                       <  1:372727/127‑1 (MQ=14)
                                                                                                                                    gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACAGTTGGGCATCGTGTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATc                 >  3:262465/1‑149 (MQ=14)
                                                                                                                                    gCCGGATAAGGCGCTCGCGCCGCATCCGGCCGACAGTTGAGCATCGTGTGAGTCCCTGTCGGGCGGATAAGATGCGCAAGTATCGCATCCGACCTTATTGCGGCCCTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATc                 >  5:165547/1‑149 (MQ=14)
                                                                                                                                      atgaagcggccCGAGCTCCTTATCCGGCCTACAGTTGCGCATCCACTGATTCACTGTAGGTCGGATAAGATTCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATTCGGCGCGAGCGCCTTATcc                <  7:182885/143‑1 (MQ=11)
                                                                                                                                                  gAGCGCCTTATCCGGCCTACAGTTGCGCAGCGTGTGAGTCCCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg   <  8:177902/149‑1 (MQ=255)
                                                                                                                                                   aGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTggg  >  1:10416/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGG  >  NC_000913/4153281‑4153576

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