Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,330,907 C→A 5.6% S365S (TCG→TCT nspC ← carboxynorspermidine/carboxyspermidine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,330,9070CA 5.6% 236.1 / 3.6 72S365S (TCG→TCTnspCcarboxynorspermidine/carboxyspermidine decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base A (2/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

ATGTTCCGTGACGCGATGGCGATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAA  >  NC_002947/3330759‑3331050
                                                                                                                                                    |                                                                                                                                               
aTGTTCCGTGACGCGATGGCGATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTAc                                                                                                                                                 <  3:271317/149‑1 (MQ=255)
 tGTTCCGTGACGCGATGGCGATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTAc                                                                                                                                                 <  8:267588/148‑1 (MQ=255)
   ttCCGTGACGCGATGGCGATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGa                                                                                                                                               <  6:210969/148‑1 (MQ=255)
      cGTGAAGCGATGGCTATACGATCCTGTTCGTCTTTATGCTTAGCGCACTTGTGGGCAACCACCTTGGCGACACATCCTGCAACAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTTAAAGCCGGTATTACGACAgg                                                                                                                                           <  4:159022/149‑1 (MQ=255)
       gTGACGCGATGGCGATACGACCCAGATCGTCCTTATGCTGCGCGCACTTGTGGGCCACCACCTTGTCGACACCTCCTGCACCAATGATAAGAACGATCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTATTACGACAGGc                                                                                                                                          >  6:289721/1‑149 (MQ=255)
        tGACGCGATGGCGATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCt                                                                                                                                         >  3:45064/1‑149 (MQ=255)
           cgcgATGGCGATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCTGTCGTAAGACCGGCtttc                                                                                                                                      <  8:111812/149‑4 (MQ=255)
              gATGGCGATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAAc                                                                                                                                   <  1:22506/149‑1 (MQ=255)
              gATGGCGATACGACCCAGCTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAAc                                                                                                                                   <  6:348946/149‑1 (MQ=255)
                  gCGATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTa                                                                                                                                <  4:314830/148‑1 (MQ=255)
                  gCGATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTa                                                                                                                                <  5:68500/148‑1 (MQ=255)
                    gATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCTGTCGTAAGACCGGCtttc                                                                                                                                      >  4:340682/1‑137 (MQ=255)
                    gATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCTGTCGTAAGACCGGCtttc                                                                                                                                      <  3:340682/140‑4 (MQ=255)
                    gATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTc                                                                                                                             >  1:183691/1‑149 (MQ=255)
                         gACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGTGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCg                                                                                                                         >  5:357349/1‑148 (MQ=255)
                            ccAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAAc                                                                                                                     >  8:85554/1‑149 (MQ=255)
                                ttCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGaa                                                                                                                 >  1:217046/1‑149 (MQ=255)
                                  cgtcgtTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACt                                                                                                               >  6:8557/1‑149 (MQ=255)
                                  cgtcgtTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACt                                                                                                               >  2:111284/1‑149 (MQ=255)
                                    tcgtTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTcg                                                                                                             >  8:73717/1‑149 (MQ=255)
                                                 cgcACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGa                                                                                                >  2:130760/1‑149 (MQ=255)
                                                         tgGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCg                                                                                        >  7:48828/1‑149 (MQ=255)
                                                           ggCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCTGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCTCGTACGACTTCGACGCTGCCGTc                                                                                      <  8:22914/149‑1 (MQ=255)
                                                           ggCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTc                                                                                      >  6:327498/1‑149 (MQ=255)
                                                            gCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCg                                                                                     <  2:183691/149‑1 (MQ=255)
                                                             ccaccaCCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTc                                                                                      >  7:286210/1‑147 (MQ=255)
                                                              caccacCTTGGCGACACCTCCTGCACTAATGATAAGAACGTTCTTCTTCAATTTACCATTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAg                                                                                   >  1:284872/1‑149 (MQ=255)
                                                               accaccTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTc                                                                                      >  4:1514/1‑145 (MQ=255)
                                                                 caccTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAgctg                                                                                >  8:21861/1‑149 (MQ=255)
                                                                    cTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCtt                                                                             >  1:247488/1‑149 (MQ=255)
                                                                     ttGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAg                                                                                   >  3:166889/1‑142 (MQ=255)
                                                                     ttGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAg                                                                                   <  4:166889/142‑1 (MQ=255)
                                                                       ggCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTg                                                                            <  1:130760/147‑1 (MQ=255)
                                                                                tcctGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATgg                                                                  >  3:366827/1‑148 (MQ=255)
                                                                                  ctGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGaa                                                                <  7:73717/148‑1 (MQ=255)
                                                                                  ctGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGaa                                                                <  2:332120/148‑1 (MQ=255)
                                                                                  ctGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAg                                                               <  6:357349/149‑1 (MQ=255)
                                                                                      aCCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCAt                                                           <  5:8557/149‑1 (MQ=255)
                                                                                        cAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGAACACGATGGAAGGCATCt                                                         >  7:125865/1‑149 (MQ=255)
                                                                                         aaTGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCtt                                                        <  5:327498/149‑1 (MQ=255)
                                                                                                   aCGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGtt                                               >  1:15067/1‑148 (MQ=255)
                                                                                                      ttcttcttcAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAAc                                           >  4:322174/1‑149 (MQ=255)
                                                                                                           cttcAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCa                                         <  8:286210/146‑1 (MQ=255)
                                                                                                           cttcAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGtt                                      >  1:102361/1‑149 (MQ=255)
                                                                                                              cAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGtttt                                    <  1:293316/148‑1 (MQ=255)
                                                                                                                aTTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCt                                 <  3:1514/149‑1 (MQ=255)
                                                                                                                 tttACCGTTTGCTCCTTACTTGATAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGTCGACTTCGACGCTGGCGTCGAGCTGCTTGACGACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCtt                                <  5:289721/149‑1 (MQ=255)
                                                                                                                     cccTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGAcc                            <  2:217046/149‑1 (MQ=255)
                                                                                                                       cTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCAt                          <  2:102361/149‑1 (MQ=255)
                                                                                                                            cTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGaa                                                                >  4:150149/1‑106 (MQ=255)
                                                                                                                            cTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGaa                                                                <  3:150149/106‑1 (MQ=255)
                                                                                                                             tCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGgtgt                     >  7:52952/1‑148 (MQ=255)
                                                                                                                             tCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGgtgt                     >  5:30126/1‑148 (MQ=255)
                                                                                                                             tCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGgtgt                     >  1:181682/1‑148 (MQ=255)
                                                                                                                               cTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAAc                  <  4:45064/149‑1 (MQ=255)
                                                                                                                               cTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAAc                  >  5:101037/1‑149 (MQ=255)
                                                                                                                                 tACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAAccc                >  6:197795/1‑149 (MQ=255)
                                                                                                                                  aCTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATg                         >  8:96604/1‑139 (MQ=255)
                                                                                                                                  aCTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATg                         <  7:96604/139‑1 (MQ=255)
                                                                                                                                  aCTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACccg               >  6:11627/1‑149 (MQ=255)
                                                                                                                                   cTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACccgc              >  5:338030/1‑149 (MQ=255)
                                                                                                                                    ttGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGa                              >  5:312415/1‑132 (MQ=255)
                                                                                                                                    ttGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGa                              <  6:312415/132‑1 (MQ=255)
                                                                                                                                    ttGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACccgcc             >  6:65017/1‑149 (MQ=255)
                                                                                                                                    ttGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACccgcc             <  7:85554/149‑1 (MQ=255)
                                                                                                                                      gAAAGCCGGTCTTACGTCAGGCTGGATACGTAGTCTTCGAAACCGAACTCGCGGTCGAGTTAGACGCTGCCGTCGAGCTGGTTGACGACGATGGTAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTgt                              <  5:311014/130‑2 (MQ=255)
                                                                                                                                      gAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACATCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCttca                              >  6:311014/1‑128 (MQ=255)
                                                                                                                                       ccaGCCTGTCGTAAGACCGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCgtg          >  5:47075/3‑149 (MQ=255)
                                                                                                                                       aaaGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGttttt                                   >  2:153451/1‑124 (MQ=255)
                                                                                                                                       aaaGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGttttt                                   <  1:153451/124‑1 (MQ=255)
                                                                                                                                            cGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGTCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGAt     <  5:65017/149‑1 (MQ=255)
                                                                                                                                            cGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGAt     <  6:101037/149‑1 (MQ=255)
                                                                                                                                             ggTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGAt     >  2:306457/1‑148 (MQ=255)
                                                                                                                                               tCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGaa  <  5:197795/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ATGTTCCGTGACGCGATGGCGATACGACCCAGTTCGTCGTTATGCTGCGCGCACTTGTGGGCCACCACCTTGGCGACACCTCCTGCACCAATGATAAGAACGTTCTTCTTCAATTTACCCTTTGCTCCTTACTTGAAAGCCGGTCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAA  >  NC_002947/3330759‑3331050

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: