Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,914,128 A→G 7.5% intergenic (+20/+133) 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,1280AG 7.5% 188.6 / 4.8 69intergenic (+20/+133)spuI/argAglutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base G (3/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

AGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTC  >  NC_002947/5913981‑5914275
                                                                                                                                                   |                                                                                                                                                   
agaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                    <  3:69864/149‑1 (MQ=255)
    ccTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGggcg                                                                                                                                                <  2:251422/149‑1 (MQ=255)
     cTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcg                                                                                                                                                <  7:62478/148‑1 (MQ=255)
        aaCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcggc                                                                                                                                              <  3:86320/147‑1 (MQ=255)
            aGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCACCCCTCAAGGGCGGCGttttt                                                                                                                                        <  3:90930/149‑1 (MQ=255)
                cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTCGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                    <  8:296510/149‑1 (MQ=255)
                 aTCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                    >  8:268600/1‑148 (MQ=255)
                  tCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                    >  5:273254/1‑147 (MQ=255)
                     aTCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGAc                                                                                                                               >  2:288386/1‑149 (MQ=255)
                      tCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGAcc                                                                                                                              >  8:5308/1‑149 (MQ=255)
                        cGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCgg                                                                                                                            >  1:55601/1‑149 (MQ=255)
                        cGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCg                                                                                                                             <  6:273254/148‑1 (MQ=255)
                         gAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGACCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGgt                                                                                                                           <  7:134672/149‑1 (MQ=255)
                            gTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTC                                                                                                                        >  3:359752‑M1/1‑149 (MQ=255)
                                aTCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCC                                                                                                                    >  6:50153‑M1/1‑149 (MQ=255)
                                   gACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCT                                                                                                                 <  6:236321‑M1/149‑1 (MQ=255)
                                         ttCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGG                                                                                                           >  1:86384‑M1/1‑149 (MQ=255)
                                           cgtcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGG                                                                                                           <  5:166205‑M1/147‑1 (MQ=255)
                                           cgtcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGG                                                                                                           <  2:86384‑M1/147‑1 (MQ=255)
                                            gtcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCA                                                                                                        >  8:324647‑M1/1‑149 (MQ=255)
                                                 gTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGCCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                   >  3:113134‑M1/1‑149 (MQ=255)
                                                     aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCTTGAGGGGCGGCGtt                                                                                                                                           <  6:32641/105‑1 (MQ=255)
                                                     aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCTTGAGGGGCGGCGtt                                                                                                                                           >  5:32641/1‑105 (MQ=255)
                                                          aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCa                                                                                          <  5:50153/149‑1 (MQ=255)
                                                          aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCa                                                                                          <  4:359752/149‑1 (MQ=255)
                                                          aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCa                                                                                          <  7:123294/149‑1 (MQ=255)
                                                          aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCa                                                                                          <  1:288386/149‑1 (MQ=255)
                                                                      agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCACGGGCACGCCCGCTCCCGCAGGGTAATTTCg                                                                               >  4:163199/1‑148 (MQ=255)
                                                                          ttCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                   >  7:66649‑M1/1‑124 (MQ=255)
                                                                          ttCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                   <  8:66649‑M1/124‑1 (MQ=255)
                                                                             gAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCCCCCCCGCGCCCGCCGGGGAAATTCGAAGGcgg                                                                        >  6:268268/1‑146 (MQ=255)
                                                                              aGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGaa                                                                      <  8:314524/149‑1 (MQ=255)
                                                                                 aCTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagcctt                                                                   >  1:199375‑M1/1‑123 (MQ=255)
                                                                                     gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCg                                                                               <  7:122549/133‑1 (MQ=255)
                                                                                         tcCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACt                                                            >  3:18138/1‑148 (MQ=255)
                                                                                          ccGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTa                                                          >  5:284693/1‑149 (MQ=255)
                                                                                           cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCCCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgtg                                                               <  7:5308/143‑1 (MQ=255)
                                                                                             aCCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCACGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCt                                                       <  3:163199/149‑1 (MQ=255)
                                                                                               cTCGAGTACAACTGGTACCTGCATCCCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                                     <  7:324647/149‑1 (MQ=255)
                                                                                                  gAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTGCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACGTACCTGTgg                                                   <  6:84885‑M2/148‑3 (MQ=255)
                                                                                                  gAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTg                                                    >  4:18138‑M2/1‑146 (MQ=255)
                                                                                                    gTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggag                                                <  8:17623‑M2/149‑6 (MQ=255)
                                                                                                       cAACTGGTACCTGCATACGGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGAGCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgg                                             <  5:201018‑M1/149‑49 (MQ=255)
                                                                                                       cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgg                                             >  7:17623‑M2/1‑141 (MQ=255)
                                                                                                         aCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggc                                           <  2:199375‑M1/149‑51 (MQ=255)
                                                                                                          cTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCACGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcg                                          <  5:97487‑M2/149‑12 (MQ=255)
                                                                                                            ggTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                >  1:305662/1‑125 (MQ=255)
                                                                                                            ggTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                <  2:305662/125‑1 (MQ=255)
                                                                                                             gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtg                                        >  3:88134‑M2/1‑135 (MQ=255)
                                                                                                              tACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggccgccc                                      <  3:153668‑M2/149‑16 (MQ=255)
                                                                                                                 cTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCATGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgc                                   >  7:268600‑M2/1‑131 (MQ=255)
                                                                                                                 cTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggc                                           >  8:359409‑M2/1‑131 (MQ=255)
                                                                                                                 cTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggc                                           <  7:359409‑M2/141‑11 (MQ=255)
                                                                                                                      tACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgccc                                     >  3:279171‑M2/1‑126 (MQ=255)
                                                                                                                      tACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgccc                                     <  4:279171‑M2/142‑17 (MQ=255)
                                                                                                                           tgtAAAGGAAAACGCCGCCCCTCACGGGGGGGTTTTTTGTTTGTGACCGGTGTCGCCCTATTCGCGGCCACGCCGGCTCCCGATGGTTAATTTCGACGGCTGATGCGTGTGACTTCCCTGTGGGAGGGGGCGTGCCCGCGAAGAGggcg                         <  2:16232/149‑1 (MQ=255)
                                                                                                                               aaCGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAggg                                                                                       >  3:350064/1‑83 (MQ=255)
                                                                                                                               aaCGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAggg                                                                                       <  4:350064/83‑1 (MQ=255)
                                                                                                                                  gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgg                                             <  3:312036‑M1/122‑49 (MQ=255)
                                                                                                                                  gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgg                                             >  4:312036‑M1/1‑74 (MQ=255)
                                                                                                                                  gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgg                                                   <  4:153286‑M1/116‑43 (MQ=255)
                                                                                                                                  gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgg                                                   >  3:153286‑M1/1‑74 (MQ=255)
                                                                                                                                    aaaCGCCGCCCTTGAGGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTTATTTCGAAGGCTGAAGCCTGTGTCTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAgggt                >  4:323601/1‑147 (MQ=38)
                                                                                                                                    aaaCGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGGTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGCGGGCGGCCCAgggg                >  7:63991/1‑147 (MQ=38)
                                                                                                                                     cccgccgcCCCGCCCGGGCGGCGTTTTTGTTGGAGCCCGGGGGCGCCCTCTGCGCGGGCCCGCCCGCTCCCGCCGGGTAATTTCGACGGCTGCCCCCTGTGACTTCCCTGTGGGCGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCtt               <  5:268268/147‑1 (MQ=255)
                                                                                                                                          cgccCCTCATGGGGGGTTTTTTTTTTTGAGAGCGGTGTCCCCGTGTTCGGGGCCACGCCGGCTCCGGCAGGTTATTTTCTATGGTTTATGCCTGTGAGTTCCCTGTGGTTGGGGGCTTGCCGGCGAAGAGGGCGGCCCAGGCTTACAgg          <  2:55601/149‑1 (MQ=255)
                                                                                                                                              ccTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCt      >  2:75771/1‑149 (MQ=255)
                                                                                                                                              ccTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCACGGTAATTTAGAAGGCTGAAGCCTGTGACTTAACAGTGGGCGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCATACAGGGGCt      >  2:304770/1‑149 (MQ=39)
                                                                                                                                              ccTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTAGCGGGCCCCCCCCCGCCCGCCGGGGAATTTTGAAGGCGTGACCCGGTGGGGTACCGGTGGGGGGGGGCGTGGCCCCGGAGAGGGCGGCCCCGGGGTAAAgggggg      >  4:341417/1‑147 (MQ=255)
                                                                                                                                                  aaGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCCCGCCCGCAGGGGAATTTAGAAGGCTGAAGCCTGTGACGTACCTGTGGGAGCGGGCGTTCCCGCGAAGAGGGCGGCCCAGGGTTTAAGGGGCttgtg  >  1:129296/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTC  >  NC_002947/5913981‑5914275

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