Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,492,677 C→T 6.4% G193E (GGG→GAG)  frc ← formyl‑CoA transferase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,492,6770CT 6.4% 201.6 / 2.5 64G193E (GGG→GAG) frcformyl‑CoA transferase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (35/24);  new base T (2/2);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

CGCGGGGAACTGCATCACCAAATGTACCATTCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAGTGGATGCCGCGCCACCCG  >  NC_000913/2492531‑2492818
                                                                                                                                                  |                                                                                                                                             
cgcgGGGAACTGCATCACCAAATGTACCATTCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCAGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGccc                                                                                                                                              <  8:287572/148‑1 (MQ=255)
 gcgGGGAACTGCATCACCAAATGTACCATTCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCa                                                                                                                                             >  7:144418/1‑148 (MQ=255)
           gCATCACCAAATGTACCATTCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGCCGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCAACCTGCATTGACATGGTGACTCGTTGCTCACGCCCCGttt                                                                                                                                   >  2:231326/1‑148 (MQ=255)
             aTCACCAAATGTACCATTCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGttttt                                                                                                                                 >  4:208690/1‑148 (MQ=255)
              tCACCAAATGTACCATTCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGtttttt                                                                                                                                >  6:300615/1‑148 (MQ=255)
               cACCAAATGTACCATTCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTc                                                                                                                               >  5:8164/1‑148 (MQ=255)
                  cAAATGTACCATTCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTcgc                                                                                                                             >  4:73739/1‑147 (MQ=255)
                      tGTACCATTCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATgc                                                                                                                        <  3:292526/148‑1 (MQ=255)
                              tCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAag                                                                                                                 >  6:200785/1‑147 (MQ=255)
                              tCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAag                                                                                                                 >  2:67015/1‑147 (MQ=255)
                               cGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAagc                                                                                                                <  3:73739/147‑1 (MQ=255)
                                  aTACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAgcagca                                                                                                            <  8:144418/148‑1 (MQ=255)
                                     cTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGt                                                                                                         <  3:286409/148‑1 (MQ=255)
                                         gggTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAAc                                                                                                     <  3:208690/148‑1 (MQ=255)
                                          ggTATTCGTCCAGATAAACCAATTTATCGAGACGCTGCAGGTCTCGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACTTGGAGACTCGATGCGCACGCCCCGCTATTTCACGATGCTGCAAAGCATCAAGTAATcc                                                                                                    >  6:163046/1‑148 (MQ=255)
                                           gTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAAcc                                                                                                    >  2:384854/1‑147 (MQ=255)
                                            tATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGa                                                                                                  <  1:67015/148‑1 (MQ=255)
                                                  tcCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGc                                                                                             >  5:367685/1‑147 (MQ=255)
                                                  tcCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTAACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGAATCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGc                                                                                             >  7:118644/1‑147 (MQ=255)
                                                    cAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCaa                                                                                           <  2:284999/147‑1 (MQ=255)
                                                      gATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATg                                                                                        >  7:93039/1‑148 (MQ=255)
                                                        tAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCa                                                                                      <  5:200785/148‑1 (MQ=255)
                                                            ccAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTcc                                                                                  >  3:368021/1‑148 (MQ=255)
                                                            ccAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTcc                                                                                  >  2:363790/1‑148 (MQ=255)
                                                             aaaTTTATCGAGACGCTGCTCGTCACGTAATTTCACGCGGCAAAGCTTCAATACGGCATCCTGCATTGACATCGTCACTCGTTGCCCACGCCCCGTTTTTTCGCGATTCAGCAAAGAAGCAAGTAAACCGATCAGCAAATGCATTCCg                                                                                 >  7:355421/2‑148 (MQ=255)
                                                                tttATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGgtg                                                                              <  5:300615/148‑1 (MQ=255)
                                                                tttATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGgtg                                                                              <  6:8164/148‑1 (MQ=255)
                                                                        gaCGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGt                                                                       >  3:187774/1‑147 (MQ=255)
                                                                            ctgctgGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCAcc                                                                   >  6:20950/1‑147 (MQ=255)
                                                                               ctgGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCaa                                                                <  4:187774/147‑1 (MQ=255)
                                                                               ctgGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAAc                                                               >  3:324517/1‑148 (MQ=255)
                                                                                      cGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAAc                                                                                                     <  2:82758/103‑1 (MQ=255)
                                                                                      cGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAAc                                                                                                     >  1:82758/1‑103 (MQ=255)
                                                                                      cGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTTTTGCTGTCACACAACGCTGCAg                                                        >  3:66284/1‑148 (MQ=255)
                                                                                            ttCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCTCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAgcg                                                      <  7:57238/144‑1 (MQ=255)
                                                                                            ttCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCTCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAgcg                                                      >  8:57238/1‑144 (MQ=255)
                                                                                                      aaaGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCGCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAATCCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGcggcgg                                        <  8:355421/148‑1 (MQ=255)
                                                                                                           ttCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCAt                                   <  4:324517/148‑1 (MQ=255)
                                                                                                                tACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCa                               >  3:20223/1‑147 (MQ=255)
                                                                                                                      aTCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAAccc                        <  1:384854/148‑1 (MQ=255)
                                                                                                                       tcatGCATTGACATTGTGACTAGTTGACCACGCCCAGTTTTTTCGGGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGAACATCCCAAAAACCCg                       <  4:66284/145‑1 (MQ=255)
                                                                                                                        ccTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATg                                                                                                                         >  4:134174/1‑49 (MQ=255)
                                                                                                                        ccTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATg                                                                                                                         <  3:134174/49‑1 (MQ=255)
                                                                                                                        ccTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAACACCCg                       >  4:57269/1‑147 (MQ=255)
                                                                                                                        ccTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCg                       >  2:151074/1‑147 (MQ=255)
                                                                                                                           gCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAGt                   >  6:219452/1‑148 (MQ=255)
                                                                                                                            cATTGACATGGTGACTCGTTGCTCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAg                    <  1:231326/146‑1 (MQ=255)
                                                                                                                            cATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCaa                                                                                           >  2:325919/1‑75 (MQ=255)
                                                                                                                            cATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCaa                                                                                           <  1:325919/75‑1 (MQ=255)
                                                                                                                             aTTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAgcggc                                          >  7:352504/1‑123 (MQ=255)
                                                                                                                             aTTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAgcggc                                          <  8:352504/123‑1 (MQ=255)
                                                                                                                              ttGACATGGTGAATCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAGTgg                 <  3:241342/147‑1 (MQ=255)
                                                                                                                               tGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACAATCCCAAAAACCCGGAGAGGa                >  4:241342/1‑147 (MQ=255)
                                                                                                                                gACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCaaaa                            <  4:272866/134‑1 (MQ=255)
                                                                                                                                gACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCaaaa                            >  3:272866/1‑134 (MQ=255)
                                                                                                                                 aCATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAGTGGATg              <  5:219452/147‑1 (MQ=255)
                                                                                                                                     ggTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAgcgc                                                     <  6:104819/104‑1 (MQ=255)
                                                                                                                                     ggTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAgcgc                                                     >  5:104819/1‑104 (MQ=255)
                                                                                                                                     ggTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGAAAAGCAGCTAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAGTGGATGCcgcg         >  1:8124/1‑148 (MQ=255)
                                                                                                                                       tGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAGTGGATGCCgcgc        >  6:367633/1‑147 (MQ=255)
                                                                                                                                          cTCGTTGCCCACGCCCCGTTTTTTCTCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACAAGCGGCGGACCATACCAAAAACCCGTAGTGGATGCCGCGCCAcc    >  8:7965/1‑148 (MQ=255)
                                                                                                                                          cTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAGTGGATGCCGCGCCAcc    >  4:69690/1‑148 (MQ=255)
                                                                                                                                          cTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATACCAAAAACCCGTAGTGGATGCCGCGCCAcc    >  1:239104/1‑148 (MQ=255)
                                                                                                                                            cGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAGTGGATGCCGCGCCACCCg  >  8:11898/1‑148 (MQ=255)
                                                                                                                                             gTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCAGCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAGTGGATGCCGCGCCACCCg  <  4:368021/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CGCGGGGAACTGCATCACCAAATGTACCATTCGGATACTGCGGGTATTCTTCCAGATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCCTGCATTGACATGGTGACTCGTTGCCCACGCCCCGTTTTTTCGCGATGCAGCAAAGCAGCAAGTAAACCGATCAGCAAATGCATTCCGGTGTTGCTGTCACCCAACGCTGCAGCGCTTACCAGCGGCGGACCATCCCAAAAACCCGTAGTGGATGCCGCGCCACCCG  >  NC_000913/2492531‑2492818

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