Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,090,047 A→G 10.3% L297P (CTG→CCG)  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,0470AG10.3% 177.9 / 10.8 68L297P (CTG→CCG) astA‑IIarginine N‑succinyltransferase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (3/4);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

GTCTGTTGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCTC  >  NC_002947/5089902‑5090195
                                                                                                                                                 |                                                                                                                                                    
gTCTGTTGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGc                                                                                                                                                   <  2:357560/149‑1 (MQ=255)
     ttGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTg                                                                                                                                              >  2:92017/1‑149 (MQ=255)
     ttGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGACTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGtt                                                                                                                                               <  1:67274/148‑1 (MQ=255)
        ccAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACGGCCTGTTGca                                                                                                                                            <  5:138373/148‑1 (MQ=255)
        ccAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACGGCCTGTTGca                                                                                                                                            <  3:104028/148‑1 (MQ=255)
          aGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACGGCCTGTTGcagga                                                                                                                                         <  5:205083/149‑4 (MQ=255)
              cTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAgg                                                                                                                                     >  8:316794/1‑149 (MQ=255)
                 aaaCCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAgg                                                                                                                                                    >  6:181794/1‑131 (MQ=255)
                 aaaCCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAgg                                                                                                                                                    <  5:181794/131‑1 (MQ=255)
                   aCCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACgg                                                                                                                                <  5:94616/149‑1 (MQ=255)
                   aCCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACgg                                                                                                                                <  7:56614/149‑1 (MQ=255)
                     cgcgACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACGGCCTGTTGcag                                                                                                                                           >  8:307515/1‑135 (MQ=255)
                     cgcgACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACGGCCTGTTGcag                                                                                                                                           <  7:307515/136‑2 (MQ=255)
                      gcgACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACCCCAGGTACGGgcg                                                                                                                             >  3:229032/1‑149 (MQ=255)
                      gcgACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGgcg                                                                                                                             >  6:7453/1‑149 (MQ=255)
                              gcgcACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACAGCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCtt                                                                                                                     >  1:333289/1‑149 (MQ=255)
                              gcgcACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAg                                                                                                                                      <  5:297832/132‑1 (MQ=255)
                              gcgcACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAg                                                                                                                                      >  6:297832/1‑132 (MQ=255)
                                gcACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACgg                                                                                                                                >  5:329740/1‑136 (MQ=255)
                                gcACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACgg                                                                                                                                <  6:329740/136‑1 (MQ=255)
                                    gCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCg                                                                                                               >  1:94142/1‑149 (MQ=255)
                                      tGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCggg                                                                                                             <  2:94142/149‑1 (MQ=255)
                                        gcgcCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCggg                                                                                                             <  6:42030/147‑1 (MQ=255)
                                             cTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCtt                                                                                                      <  8:223864/149‑1 (MQ=255)
                                                aCCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGTCTGACCGGATTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGa                                                                                                   >  2:369113/1‑149 (MQ=255)
                                                  ccACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATc                                                                                                 >  6:130036/1‑149 (MQ=255)
                                                  ccACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATc                                                                                                 <  6:267522/149‑1 (MQ=255)
                                                  ccACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATc                                                                                                 >  5:166974/1‑149 (MQ=255)
                                                     cGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTc                                                                                              <  3:98433/149‑1 (MQ=255)
                                                       cccAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCggg                                                                                                             >  6:360636/1‑132 (MQ=255)
                                                       cccAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCggg                                                                                                             <  5:360636/132‑1 (MQ=255)
                                                        ccAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCAcc                                                                                           >  5:197533/1‑149 (MQ=255)
                                                        ccAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCAcc                                                                                           <  5:221557/149‑1 (MQ=255)
                                                               gTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGcacca                                                                                    >  8:74367/1‑149 (MQ=255)
                                                                     tgctgcCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATc                                                                                                 <  1:393600/130‑1 (MQ=255)
                                                                     tgctgcCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATc                                                                                                 >  2:393600/1‑130 (MQ=255)
                                                                          ccACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTgg                                                                         <  1:92017/149‑1 (MQ=255)
                                                                          ccACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTgg                                                                         >  7:13124/1‑149 (MQ=255)
                                                                           cACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTggg                                                                        <  6:197533/149‑1 (MQ=255)
                                                                            aCGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTggg                                                                        >  4:7132/1‑148 (MQ=255)
                                                                            aCGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTggg                                                                        <  3:7132/148‑1 (MQ=255)
                                                                                tCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTCCAACAGGCCGTTGCACACAAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGCCTCTCGGCGATc                                                                   >  2:334225/1‑149 (MQ=255)
                                                                                   ggCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAg                                                                <  5:224418/149‑1 (MQ=255)
                                                                                    gCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCgg                                                                                                              >  4:280219/1‑102 (MQ=255)
                                                                                    gCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCgg                                                                                                              <  3:280219/102‑1 (MQ=255)
                                                                                    gCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGc                                                               >  3:26132/1‑149 (MQ=255)
                                                                                      tGACCGGCTTGCCGGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCg                                                              >  7:384381/1‑148 (MQ=255)
                                                                                         ccGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATg                                                          >  6:406205/1‑149 (MQ=255)
                                                                                                   aGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTgc                                                >  8:385139/1‑149 (MQ=255)
                                                                                                     gCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTgcgc                                              >  1:105434/1‑149 (MQ=255)
                                                                                                       aCCCAGTCCAGGTCCAGCAACCCCGCGCGGTAGTCATCCAACAGCCCGTTCCACACCAGTTACGGGCGCCCGCCTTTGCCCGGGGCATCTACGATCTTCACCGGCCCCCCCCGCCAATGGGCGATCAAGCGGATGCCACAGGTgcgcg                                             >  6:38547/1‑148 (MQ=255)
                                                                                                         ccAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCAt                                          <  2:105434/149‑1 (MQ=255)
                                                                                                            gtccaggtccagCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTg                                                                                                                                              >  8:89566/1‑46 (MQ=255)
                                                                                                            gtccaggtccagCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTg                                                                                                                                              <  7:89566/46‑1 (MQ=255)
                                                                                                                aggtccagCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGt                                   <  7:316794/149‑1 (MQ=255)
                                                                                                                     cagCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCggg                               <  6:166974/148‑1 (MQ=255)
                                                                                                                      agCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGgcc                             <  7:74367/149‑1 (MQ=255)
                                                                                                                       gCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGgccg                            >  8:225759/1‑149 (MQ=255)
                                                                                                                               gcgcgGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGaa                    <  4:26132/149‑1 (MQ=255)
                                                                                                                                   ggTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCTGCTCTGGGCGATCGAGCGGATGCCAGAGCTGCGCGCATCCAGGGTCGGGCCGCCGTCGAAGATg                >  5:78551/1‑149 (MQ=255)
                                                                                                                                      tggtgtGCAACGGCCTGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGCTGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTc              >  8:298711/6‑148 (MQ=255)
                                                                                                                                       gTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGa            <  5:7453/149‑1 (MQ=255)
                                                                                                                                          gtGCAACGGCCTGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGt         >  3:53973/2‑149 (MQ=255)
                                                                                                                                          cTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGACCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATg          >  3:242612/1‑148 (MQ=255)
                                                                                                                                             cAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCATAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAg       <  1:369113/148‑1 (MQ=255)
                                                                                                                                              aaCAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGa                  <  2:69906/136‑1 (MQ=255)
                                                                                                                                              aaCAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGa                  >  1:69906/1‑136 (MQ=255)
                                                                                                                                              aaCAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTg     >  2:20403/1‑149 (MQ=255)
                                                                                                                                                 aGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCTc  >  6:218690/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTCTGTTGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCAGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCTC  >  NC_002947/5089902‑5090195

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