Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 519,337 C→G 7.0% I242M (ATC→ATG)  argC → N‑acetyl‑gamma‑glutamyl‑phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947519,3370CG 7.0% 227.2 / 6.4 71I242M (ATC→ATG) argCN‑acetyl‑gamma‑glutamyl‑phosphate reductase
Reads supporting (aligned to +/- strand):  ref base C (36/30);  new base G (2/3);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

TTCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCTGC  >  NC_002947/519191‑519484
                                                                                                                                                  |                                                                                                                                                   
ttCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCa                                                                                                                                                   >  4:355135/1‑149 (MQ=255)
  cTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGGTCCGTGGCATCCAt                                                                                                                                                  <  4:31829/148‑1 (MQ=255)
  cTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCAt                                                                                                                                                  <  7:126779/148‑1 (MQ=255)
    gCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGc                                                                                                                                                >  4:231243/1‑148 (MQ=255)
     cGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATg                                                                                                                                                 <  2:347555/146‑1 (MQ=255)
     cGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCa                                                                                                                                              >  4:165379/1‑149 (MQ=255)
      gAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATGGATGCCAc                                                                                                                                             <  4:209824/149‑1 (MQ=255)
       aGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCAcc                                                                                                                                            <  3:355135/149‑1 (MQ=255)
           cGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGt                                                                                                                                        <  8:9036/149‑1 (MQ=255)
             gCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATGGATGCCAcggatcat                                                                                                                                      <  8:254551/149‑8 (MQ=255)
                aaaGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTAcgc                                                                                                                                    >  1:267473/1‑148 (MQ=255)
                     aTGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGt                                                                                                                              >  5:353410/1‑149 (MQ=255)
                      tGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATTGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTg                                                                                                                             >  1:259048/1‑149 (MQ=255)
                        aaaGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTgg                                                                                                                            <  3:165379/148‑1 (MQ=255)
                         aaGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGt                                                                                                                           <  2:267473/148‑1 (MQ=255)
                            ccTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGa                                                                                                                        >  3:311242/1‑148 (MQ=255)
                             cTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATa                                                                                                                      >  5:271590/1‑149 (MQ=255)
                             cTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATa                                                                                                                      >  8:280954/1‑149 (MQ=255)
                              tACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATAc                                                                                                                     >  7:159452/1‑149 (MQ=255)
                                cgcgGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATAc                                                                                                                     <  3:231243/147‑1 (MQ=255)
                                 gcgGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTc                                                                                                                  >  4:239095/1‑149 (MQ=255)
                                        aGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGAc                                                                                                           <  8:81157/149‑1 (MQ=255)
                                        aGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGAc                                                                                                           <  8:204512/149‑1 (MQ=255)
                                            ccATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGc                                                                                                       <  3:327427/149‑1 (MQ=255)
                                              aTCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCAGATGATCCGTGTCATCCATGCCACACTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAg                                                                                                     >  3:352119/1‑149 (MQ=255)
                                               tCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAgg                                                                                                    <  6:271590/149‑1 (MQ=255)
                                                    aCTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTg                                                                                               <  3:88416/149‑1 (MQ=255)
                                                             agaTCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCCGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGaaa                                                                                      <  6:270881/149‑1 (MQ=255)
                                                               aTCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATAcc                                                                                                                    >  6:343857/1‑117 (MQ=255)
                                                               aTCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATAcc                                                                                                                    <  5:343857/117‑1 (MQ=255)
                                                                tCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGc                                                                                   >  8:217019/1‑149 (MQ=255)
                                                                    ccAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACg                                                                               >  3:4003/1‑149 (MQ=255)
                                                                        ggCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGa                                                                           <  3:198752/149‑1 (MQ=255)
                                                                        ggCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGa                                                                           <  5:245574/149‑1 (MQ=255)
                                                                          ccTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGAt                                                                          >  6:255006/1‑148 (MQ=255)
                                                                           cTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATg                                                                         >  8:272757/1‑148 (MQ=255)
                                                                              cgcCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGaa                                                                       <  8:159452/147‑1 (MQ=255)
                                                                                    gCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTtcgt                                                               <  4:4003/149‑1 (MQ=255)
                                                                                      cgcgTGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCAGTGGCATCCATGACACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGACGATGAACCGTtcgtc                                                              >  6:181665/1‑148 (MQ=255)
                                                                                      cgcgGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTtcgtc                                                              >  1:196565/1‑148 (MQ=255)
                                                                                        cgGGCAAGGACATTGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGAc                                                           <  2:259048/149‑1 (MQ=255)
                                                                                        cgGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGAc                                                           <  7:280954/149‑1 (MQ=255)
                                                                                           gCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGt                                                         >  3:22741/1‑148 (MQ=255)
                                                                                            cAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGGCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGa                                                       <  4:51526/149‑1 (MQ=255)
                                                                                                gACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGcc                                                   <  3:239095/149‑1 (MQ=255)
                                                                                                gACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGc                                                    >  4:37399/1‑148 (MQ=255)
                                                                                                   aTCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCcag                                                 >  8:319887/1‑148 (MQ=255)
                                                                                                     cGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCcagcag                                              >  8:67800/1‑149 (MQ=255)
                                                                                                       gCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggc                                            <  2:242196/149‑1 (MQ=255)
                                                                                                        ccTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTCTAAGCGAACGTGGTCGATACCTCTGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcagg                                             >  5:261856/1‑147 (MQ=255)
                                                                                                         cTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggcag                                          >  5:166772/1‑149 (MQ=255)
                                                                                                         cTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggcag                                          <  7:272757/149‑1 (MQ=255)
                                                                                                         cTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggcag                                          <  7:319887/149‑1 (MQ=255)
                                                                                                              tttCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATc                                     >  5:143092/1‑149 (MQ=255)
                                                                                                                  gTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCgg                                  >  8:239012/1‑148 (MQ=255)
                                                                                                                   tGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCGGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGa                                 >  4:123099/1‑148 (MQ=255)
                                                                                                                   tGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTCCGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGa                                 >  8:155443/1‑148 (MQ=255)
                                                                                                                       gCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAAccc                            <  6:353410/149‑1 (MQ=255)
                                                                                                                              tctgcGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGGCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAgcgt                     <  6:166772/145‑1 (MQ=255)
                                                                                                                               cGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAgcgtg                    >  3:292295/1‑149 (MQ=255)
                                                                                                                                gCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAgcgtgc                   >  6:302811/1‑149 (MQ=255)
                                                                                                                                  cGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAgcgtgcgt                 >  7:114298/1‑149 (MQ=255)
                                                                                                                                  cGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAgcgtgcgt                 >  3:51690/1‑149 (MQ=255)
                                                                                                                                  cGATGATCCGTGGCATCCACGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAgcgtgcgt                 >  3:101989/1‑149 (MQ=255)
                                                                                                                                       atCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGc            <  7:67800/149‑1 (MQ=255)
                                                                                                                                        tCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGc            >  4:270021/1‑148 (MQ=255)
                                                                                                                                        agggTGGCATGGATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGc            >  4:319410/4‑148 (MQ=255)
                                                                                                                                          ggTGGCATGGATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCt                                                                                                >  7:221191/2‑62 (MQ=255)
                                                                                                                                          ggTGGCATGGATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCt                                                                                                <  8:221191/61‑1 (MQ=255)
                                                                                                                                               cATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCt    <  4:22741/149‑1 (MQ=255)
                                                                                                                                                 tCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCTGc  >  6:14185/1‑149 (MQ=255)
                                                                                                                                                  ccATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggca                                           >  4:249806/1‑107 (MQ=255)
                                                                                                                                                  ccATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggca                                           <  3:249806/107‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCTGC  >  NC_002947/519191‑519484

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