Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,627,339 C→A 8.8% E55* (GAA→TAA) ‡ PP_4944 ← carbamoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,627,3390CA 8.8% 207.7 / 8.7 68E55* (GAA→TAA) ‡PP_4944carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (33/29);  new base A (4/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

CGATTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCGGACAGGCC  >  NC_002947/5627191‑5627483
                                                                                                                                                    |                                                                                                                                                
cgaTTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTc                                                                                                                                                  <  7:280835/149‑1 (MQ=255)
  aTTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCa                                                                                                                                                <  8:185538/149‑1 (MQ=255)
   ttGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAg                                                                                                                                               <  2:106359/149‑1 (MQ=255)
      cccATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGc                                                                                                                                             <  2:350547/148‑1 (MQ=255)
        cATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAg                                                                                                                                           >  6:98183/1‑148 (MQ=255)
        cATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGGAACAGGCAgg                                                                                                                                          <  5:202327/149‑2 (MQ=255)
         aTTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCTGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGGAACAGGCAggc                                                                                                                                         <  3:245233/149‑3 (MQ=255)
                tGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGcc                                                                                                                                  >  4:250771/1‑149 (MQ=255)
                tGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGcc                                                                                                                                  >  8:115522/1‑149 (MQ=255)
                   cATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCgcc                                                                                                                               >  6:257510/1‑147 (MQ=255)
                    aTCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACt                                                                                                                              <  4:196032/149‑1 (MQ=255)
                      cGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCg                                                                                                                            >  7:149731/1‑149 (MQ=255)
                         gTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCa                                                                                                                                                <  3:334886/126‑1 (MQ=255)
                         gTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCa                                                                                                                                                >  4:334886/1‑126 (MQ=255)
                          tGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAgg                                                                                                                        <  7:115522/149‑1 (MQ=255)
                            cACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAggg                                                                                                                       >  4:147530/1‑148 (MQ=255)
                             aCGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCa                                                                                                                     >  7:70815/1‑149 (MQ=255)
                             aCGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCa                                                                                                                     >  6:274373/1‑149 (MQ=255)
                               gATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATg                                                                                                                   <  3:250771/149‑1 (MQ=255)
                               gATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATg                                                                                                                   >  5:202929/1‑149 (MQ=255)
                                        cGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGttcttc                                                                                                          <  7:23595/149‑1 (MQ=255)
                                          tACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCTTAGGGCATGCGGTTCTTCGc                                                                                                        <  2:114836/149‑1 (MQ=255)
                                               aTAGCGCTTGGAATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTAACACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCt                                                                                                   >  8:164225/1‑149 (MQ=255)
                                                tAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCtt                                                                                                  <  5:98183/149‑1 (MQ=255)
                                                  gcgcTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGt                                                                                                >  4:35726/1‑149 (MQ=255)
                                                     cTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTcac                                                                                             >  5:134286/1‑149 (MQ=255)
                                                     cTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTcac                                                                                             <  6:364005/149‑1 (MQ=255)
                                                              gTGCCACCGCGCCTTGCCGAACAGGCTGATTGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAgc                                                                                    >  6:198752/1‑149 (MQ=255)
                                                              gTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAg                                                                                     <  5:57405/148‑1 (MQ=255)
                                                                    ccGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCtctt                                                                              <  5:153493/149‑1 (MQ=255)
                                                                          cTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGc                                                                                                                                   >  1:219305/1‑90 (MQ=255)
                                                                          cTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGc                                                                                                                                   <  2:219305/90‑1 (MQ=255)
                                                                           ttGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCGCGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTcggc                                                                       <  5:318947/149‑1 (MQ=255)
                                                                               cGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGATTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTcggcggcg                                                                   <  1:110623/149‑1 (MQ=255)
                                                                                gAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTcggcggcgg                                                                  >  7:214129/1‑149 (MQ=255)
                                                                                gAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTcggcggcgg                                                                  >  2:234221/1‑149 (MQ=255)
                                                                                      gCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGttct                                                                                                            >  6:348394/1‑101 (MQ=255)
                                                                                      gCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCa                                                            <  3:147530/149‑1 (MQ=255)
                                                                                            cGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTg                                                       <  8:70815/148‑1 (MQ=255)
                                                                                                gCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGcc                                                     >  2:168035/1‑146 (MQ=255)
                                                                                                gCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGcc                                                     <  1:168035/146‑1 (MQ=255)
                                                                                                 cGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTc                                                  >  4:248028/1‑148 (MQ=255)
                                                                                                 cGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTc                                                  >  8:322318/1‑148 (MQ=255)
                                                                                                                ccACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCtctt                                                                              >  7:137711/1‑105 (MQ=255)
                                                                                                                ccACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCcgcg                                  >  1:333557/1‑149 (MQ=255)
                                                                                                                ccACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCCCATGCTTGTCACGCACGAAGCGCtctt                                                                              <  8:137711/105‑1 (MQ=255)
                                                                                                                 cACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCgg                                 <  5:198752/149‑1 (MQ=255)
                                                                                                                     tcgacgtcgGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAgg                             <  8:214129/149‑1 (MQ=255)
                                                                                                                     tcgacgtcgGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAgg                             <  7:322318/149‑1 (MQ=255)
                                                                                                                     tcgacgtcgGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGGAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAgg                             <  6:202929/149‑1 (MQ=255)
                                                                                                                          gtcgGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGTGTCAt                        <  4:369259/149‑1 (MQ=255)
                                                                                                                          gtcgGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCAt                        <  7:164225/149‑1 (MQ=255)
                                                                                                                               aCGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCt                    >  1:2019/1‑148 (MQ=255)
                                                                                                                                cGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTaa                  <  3:35726/149‑1 (MQ=255)
                                                                                                                                 ggCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCATGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAg                 >  1:14283/1‑149 (MQ=255)
                                                                                                                                 ggCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAg                 >  4:213975/1‑149 (MQ=255)
                                                                                                                                 ggCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAg                 <  6:134286/149‑1 (MQ=255)
                                                                                                                                  gCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAgg                <  2:2019/149‑1 (MQ=255)
                                                                                                                                    ttGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAACGTACAGGGCCGCGGAAGGGTCATGGCTAAGGgc              >  1:348944/1‑149 (MQ=255)
                                                                                                                                        gttctgcctgGAACAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCg          >  3:235349/4‑149 (MQ=255)
                                                                                                                                        gttctgcctgGAACAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGACg          >  5:363035/4‑149 (MQ=255)
                                                                                                                                         ttctgcctgGAACAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCg          >  3:76530/3‑148 (MQ=255)
                                                                                                                                         ttctgcctgGAACAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCg          >  2:317518/3‑148 (MQ=255)
                                                                                                                                          cctgcctgTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCGGa        >  7:313596/1‑149 (MQ=255)
                                                                                                                                           ctgcctgTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCGGa        >  2:146536/1‑148 (MQ=255)
                                                                                                                                               ctgTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCGGACAGGc   >  3:85353/1‑149 (MQ=255)
                                                                                                                                                 gTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCGGACAGGcc  >  6:130924/1‑148 (MQ=255)
                                                                                                                                                  ttCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCGGACAGGcc  >  8:308618/1‑147 (MQ=255)
                                                                                                                                                  ttCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCGGACAGGcc  >  6:230597/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGATTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGGTCATGGCTAAGGGCGCCGGACAGGCC  >  NC_002947/5627191‑5627483

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: