Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 585,109 T→A 8.0% Q49H (CAA→CAT)  ansA ← type 1 L‑asparaginase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947585,1090TA 8.0% 259.1 / 9.2 87Q49H (CAA→CAT) ansAtype 1 L‑asparaginase
Reads supporting (aligned to +/- strand):  ref base T (39/41);  new base A (3/4);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCCATGCTGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCA  >  NC_002947/584963‑585254
                                                                                                                                                   |                                                                                                                                                 
gCCATGCTGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCAtt                                                                                                                                                  <  7:337232/148‑1 (MQ=255)
gCCATGCTGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTg                                                                                                                                                 <  1:327671/149‑1 (MQ=255)
       tGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTgg                                                                                                                                          >  6:25301/1‑149 (MQ=255)
            gTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTg                                                                                                                                                 >  1:125074/1‑137 (MQ=255)
            gTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTg                                                                                                                                                 <  2:125074/137‑1 (MQ=255)
             tACCGTGTAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCAATGCAccccaggtg                                                                                                                                      <  4:88776/147‑8 (MQ=255)
             tACCGTGTAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCAATGCAccccaggtg                                                                                                                                      >  3:88776/1‑140 (MQ=255)
               ccGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCAtt                                                                                                                                  >  1:106869/1‑149 (MQ=255)
                cGTGCAGCACCAACACGCCATCGTGGCGCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTg                                                                                                                                 <  8:284917/149‑1 (MQ=255)
                cGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTg                                                                                                                                 >  4:203549/1‑149 (MQ=255)
                cGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCAATGCACCCCAGGTGCAATg                                                                                                                                 <  1:343345/149‑1 (MQ=255)
                 gTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTgg                                                                                                                                <  6:171172/149‑1 (MQ=255)
                   gcagcaCCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGcc                                                                                                                              <  3:210284/149‑1 (MQ=255)
                      gcaCCAACACGCCATCGTGGCCCGCGCCATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCTTATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCCCCTGGGGTGCATTGGCCagc                                                                                                                           >  6:150676/1‑147 (MQ=255)
                         ccAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTgg                                                                                                                        >  1:173258/1‑149 (MQ=255)
                          cAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTggg                                                                                                                       <  2:203539/149‑1 (MQ=255)
                              acGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                   <  4:202889/149‑1 (MQ=255)
                                 ccATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAgtg                                                                                                                >  1:229861/1‑149 (MQ=255)
                                   aTCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAgtgt                                                                                                               >  2:172225/1‑148 (MQ=255)
                                    tCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTc                                                                                                             <  5:25301/149‑1 (MQ=255)
                                     cGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCAATGCACCCCAGGTGCAATgga                                                                                                                               >  1:218620/1‑129 (MQ=255)
                                     cGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCAATGCACCCCAGGTGCAATgga                                                                                                                               <  2:218620/130‑2 (MQ=255)
                                         gCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCTCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                   <  5:331744/138‑1 (MQ=255)
                                         gCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                   >  6:331744/1‑138 (MQ=255)
                                          cccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTCGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCAt                                                                                                       >  8:165940/1‑149 (MQ=255)
                                           ccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                   <  4:230013/136‑1 (MQ=255)
                                           ccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                   >  3:230013/1‑136 (MQ=255)
                                           ccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATc                                                                                                      >  6:388921/1‑149 (MQ=255)
                                              cgACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGt                                                                                                   >  5:152207/1‑149 (MQ=255)
                                               gACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTAGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGt                                                                                                   <  1:299089/149‑1 (MQ=255)
                                                    ccACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTc                                                                                              >  2:321882/1‑148 (MQ=255)
                                                      aCCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaa                                                                                            >  5:12288/1‑148 (MQ=255)
                                                       ccgccgCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaaaa                                                                                          >  1:234766/1‑149 (MQ=255)
                                                       ccgccgCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaaaa                                                                                          <  6:62341/149‑1 (MQ=255)
                                                       ccgccgCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaaaa                                                                                          >  2:372082/1‑149 (MQ=255)
                                                           cgcgACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAaccac                                                                                      >  7:118146/1‑149 (MQ=255)
                                                           cgcgACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAaccac                                                                                      <  8:118146/149‑1 (MQ=255)
                                                           cgcgACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCAAATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAaccac                                                                                      >  3:280514/1‑149 (MQ=255)
                                                             cgacgaTGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCAccg                                                                                    >  4:271551/1‑149 (MQ=255)
                                                                   tGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGt                                                                              <  1:328321/149‑1 (MQ=255)
                                                                       atcgcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCa                                                                          <  2:229861/149‑1 (MQ=255)
                                                                       atcgcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCa                                                                          <  1:172225/149‑1 (MQ=255)
                                                                        tcgcccatagcCAGCCAGTTGTGCTGTTGCATATTGGCGATGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTTCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAg                                                                         <  4:101680/149‑1 (MQ=255)
                                                                          gcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAgg                                                                        <  4:109791/148‑1 (MQ=255)
                                                                            gcctcgcCAGCCAGGTGTGCTGTTGCATCTTTGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAg                                                                         <  4:96344/142‑1 (MQ=255)
                                                                               tcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTcc                                                                  >  2:183785/1‑149 (MQ=255)
                                                                                 gccagccaGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCgg                                                                <  3:271551/149‑1 (MQ=255)
                                                                                 gccagccaGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCgg                                                                <  2:173258/149‑1 (MQ=255)
                                                                                 gccagccaGTCGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCgg                                                                <  7:165940/149‑1 (MQ=255)
                                                                                      ccagTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGttt                                                           <  3:203549/149‑1 (MQ=255)
                                                                                        agTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCg                                                         <  6:12288/149‑1 (MQ=255)
                                                                                         gTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGa                                                        <  1:372082/149‑1 (MQ=255)
                                                                                          ttGTGCTGTTGCAAATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAg                                                       <  4:280514/149‑1 (MQ=255)
                                                                                            gtgCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGc                                                      <  6:152207/148‑1 (MQ=255)
                                                                                             tgCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCAt                                                    <  5:388921/149‑1 (MQ=255)
                                                                                             tgCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCAt                                                    >  3:184147/1‑149 (MQ=255)
                                                                                                    gCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCAACTTGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAACACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCAGGCCGATg                                             >  5:264149/1‑149 (MQ=255)
                                                                                                     cATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTg                                                                                                  >  5:128263/1‑95 (MQ=255)
                                                                                                     cATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTg                                                                                                  <  6:128263/95‑1 (MQ=255)
                                                                                                           ggCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCAc                                      >  4:166494/1‑149 (MQ=255)
                                                                                                              gcTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCgg                                   <  2:234766/149‑1 (MQ=255)
                                                                                                              gcTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCgg                                   >  3:235441/1‑149 (MQ=255)
                                                                                                                  tCGAGCAGCGGGTTCAGCTGGT‑TCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGta                                 <  6:264149/147‑1 (MQ=255)
                                                                                                                  tCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATg                                             <  7:13533/135‑1 (MQ=255)
                                                                                                                  tCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATg                                             >  8:13533/1‑135 (MQ=255)
                                                                                                                  tCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGta                                 <  2:106869/147‑1 (MQ=255)
                                                                                                                  ctcgacagcGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTTTTGACGCCTCCGTTCTTCAAAACCACGGGCCGGTGCCAGGGCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGtata                               >  1:253365/6‑149 (MQ=255)
                                                                                                                    gAGCAGCGGGTTCAGCTCCT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTAACTCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGtata                               >  6:67270/1‑147 (MQ=255)
                                                                                                                       cagcGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGc                                                      <  8:168327/121‑1 (MQ=255)
                                                                                                                       cagcGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGc                                                      >  7:168327/1‑121 (MQ=255)
                                                                                                                            ggTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACgaga                     <  1:183785/149‑1 (MQ=255)
                                                                                                                              ttCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACgag                      >  3:350352/1‑146 (MQ=255)
                                                                                                                               tCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACgag                      >  3:108624/1‑145 (MQ=255)
                                                                                                                                 aGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGt                                  <  5:386189/131‑1 (MQ=255)
                                                                                                                                 aGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGt                                  >  6:386189/1‑131 (MQ=255)
                                                                                                                                 aGCTCGT‑GCAGTGTCCATTGCACCATGGGTGCATTGGCCATCAGGGCAAAGTGTTCACGCATACGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATTCCGATGGTGCCACCGGTATAGAGAACGAGAAGAt                 >  7:51944/1‑148 (MQ=255)
                                                                                                                                  gCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGc                                                                       <  8:308714/93‑1 (MQ=255)
                                                                                                                                  gCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGc                                                                       >  7:308714/1‑93 (MQ=255)
                                                                                                                                  gCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCTATGGTGCCACCGGTATagag                            <  2:229512/136‑1 (MQ=255)
                                                                                                                                  gCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATagag                            >  1:229512/1‑136 (MQ=255)
                                                                                                                                       t‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCGACCGGTATAGAGAACGAGAAGATTCTTTAc          <  1:30746/149‑1 (MQ=255)
                                                                                                                                             tgtCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgcgc     >  7:104502/1‑149 (MQ=255)
                                                                                                                                               tCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgcgcgc   <  4:184147/149‑1 (MQ=255)
                                                                                                                                                ccATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCa  >  7:88086/1‑149 (MQ=255)
                                                                                                                                                ccAATGCACCCCAGGTGCAATGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGagaaga                  <  6:150011/133‑1 (MQ=255)
                                                                                                                                                ccAATGCACCCCAGGTGCAATGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGagaaga                  >  5:150011/1‑133 (MQ=255)
                                                                                                                                                  aTTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgc       >  4:324706/1‑142 (MQ=255)
                                                                                                                                                  aTTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgc       <  3:324706/142‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCATGCTGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCA  >  NC_002947/584963‑585254

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