Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,258,520:1 +G 8.1% coding (184/1005 nt) asd → aspartate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,258,5201.G 8.1% 287.0 / 19.5 74coding (184/1005 nt)asdaspartate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (31/37);  new base G (4/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCC  >  NC_002947/2258374‑2258666
                                                                                                                                                    |                                                                                                                                                   
ccGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTc                                                                                                                                                   <  1:251734/149‑1 (MQ=255)
ccGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTc                                                                                                                                                   <  8:16393/149‑1 (MQ=255)
 cGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCg                                                                                                                                                  <  5:198193/149‑1 (MQ=255)
   gCAGCGTCGGTGAAGCCTTTGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGAc                                                                                                                                                <  2:150064/149‑1 (MQ=255)
      gCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTc                                                                                                                                             >  5:17764/1‑149 (MQ=255)
      gCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTc                                                                                                                                             <  2:143521/149‑1 (MQ=255)
                  ccTTGGTGCAGATCCTCGAAGAACTGGCATTCGCCGGTTGCCACTCTGCACCTGATGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTc                                                                                                                                    >  1:201593/1‑147 (MQ=255)
                  ccTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAg                                                                                                                                 >  5:172385/1‑149 (MQ=255)
                      ggTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTgg                                                                                                                             >  7:243385/1‑149 (MQ=255)
                      ggTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTg                                                                                                                              <  1:203868/148‑1 (MQ=255)
                       gTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACATGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCTCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGc                                                                                                                            >  8:83523/1‑149 (MQ=255)
                        tGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGcc                                                                                                                           <  2:80485/149‑1 (MQ=255)
                        tGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGc                                                                                                                            <  6:172385/148‑1 (MQ=255)
                                      gaaCTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCgg                                                                                                             >  8:45050/1‑149 (MQ=255)
                                          tGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCct                                                                                                                          >  4:50060/1‑132 (MQ=255)
                                          tGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCct                                                                                                                          <  3:50060/132‑1 (MQ=255)
                                             cATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCgccgcc                                                                                                      <  8:23895/149‑1 (MQ=255)
                                             cATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCgccgcc                                                                                                      <  4:295086/149‑1 (MQ=255)
                                               ttCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGt                                                                                                    >  8:102812/1‑149 (MQ=255)
                                                          ccACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGc                                                                                          >  7:259103/1‑149 (MQ=255)
                                                           cACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCt                                                                                         >  7:182718/1‑149 (MQ=255)
                                                            aCTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCt                                                                                         >  1:116935/1‑148 (MQ=255)
                                                            aCTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCt                                                                                         <  2:116935/148‑1 (MQ=255)
                                                             ctctGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCttt                                                                                       >  5:320006/1‑149 (MQ=255)
                                                                tGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGa                                                                                                                                                 >  6:219403/1‑88 (MQ=255)
                                                                tGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGa                                                                                                                                                 <  5:219403/88‑1 (MQ=255)
                                                                tGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGcc                                                                                    >  5:234361/1‑149 (MQ=255)
                                                                     ctgctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAAGCTGTCGACTTCg                                                                                                                                            <  3:206978/88‑1 (MQ=39)
                                                                     ctgctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAAGCTGTCGACTTCg                                                                                                                                            >  4:206978/1‑88 (MQ=39)
                                                                       gctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAgg                                                                             >  4:165406/1‑149 (MQ=255)
                                                                        ctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGgc                                                                            <  7:83523/149‑1 (MQ=255)
                                                                           gCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCt                                                                         <  6:183731/149‑1 (MQ=255)
                                                                             cAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGc                                                                       >  6:9936/1‑149 (MQ=255)
                                                                             cAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGc                                                                       <  8:243385/149‑1 (MQ=255)
                                                                             cAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGc                                                                       <  6:17764/149‑1 (MQ=255)
                                                                               gCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAg                                                                     >  5:171529/1‑149 (MQ=255)
                                                                                cATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAgg                                                                    <  6:320006/149‑1 (MQ=255)
                                                                                cATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAgg                                                                    >  2:292773/1‑149 (MQ=255)
                                                                                  tGGAATAGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCg                                                                 <  3:240352/149‑1 (MQ=255)
                                                                                   ggAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCGAGGCGCTGCAGGCgg                                                                >  5:273703/1‑149 (MQ=255)
                                                                                    gAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCggg                                                               >  2:84999/1‑149 (MQ=255)
                                                                                               cagcagCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCt                                                             <  4:245387/140‑1 (MQ=255)
                                                                                               cagcagCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCt                                                             >  3:245387/1‑140 (MQ=255)
                                                                                                  cagcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCa                                                          >  6:157668/1‑140 (MQ=255)
                                                                                                  cagcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCa                                                          <  5:157668/140‑1 (MQ=255)
                                                                                                      gTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCt                                             >  3:102811/1‑149 (MQ=255)
                                                                                                        gATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGt                                           >  3:308914/1‑149 (MQ=255)
                                                                                                            ttCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCcg                                        <  5:9936/148‑1 (MQ=255)
                                                                                                                cTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCggcggc                                   >  2:103019/1‑149 (MQ=255)
                                                                                                                   gCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTg                                <  1:292773/149‑1 (MQ=255)
                                                                                                                     aaGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGga                              <  7:45050/149‑1 (MQ=255)
                                                                                                                       gAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGgacg                            >  8:268873/1‑149 (MQ=255)
                                                                                                                          aCTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGacgacg                         >  3:252240/1‑149 (MQ=255)
                                                                                                                              aaGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCGAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccct                     <  6:273703/149‑1 (MQ=255)
                                                                                                                              aaGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCa                                                          >  2:314273/1‑112 (MQ=255)
                                                                                                                              aaGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCa                                                          <  1:314273/112‑1 (MQ=255)
                                                                                                                                ggTGCGCGAAGTCGACAGCTTCGCCTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctgg                   >  2:179188/1‑149 (MQ=255)
                                                                                                                                 ctgggcgAAGTCGAAGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggc                  >  6:218466/5‑149 (MQ=255)
                                                                                                                                  tgggcgAAGTCGAAGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                               <  6:318327/116‑1 (MQ=38)
                                                                                                                                  tgggcgAAGTCGAAGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggcc                 >  4:306256/4‑149 (MQ=255)
                                                                                                                                  tgggcgAAGTCGAAGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACCGTCATCGAc                                               >  5:318327/4‑119 (MQ=38)
                                                                                                                                     gcgAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccctg              <  6:234361/149‑1 (MQ=255)
                                                                                                                                      cgAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTg                                            <  6:248187/118‑1 (MQ=255)
                                                                                                                                      cgAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTg                                            >  5:248187/1‑118 (MQ=255)
                                                                                                                                      cgAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGACGGCCTGGACGACGccctggccctgg             <  4:102811/149‑1 (MQ=255)
                                                                                                                                       gAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGccctggccctgg             <  1:84999/149‑1 (MQ=255)
                                                                                                                                           tCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCg        <  3:165406/149‑1 (MQ=255)
                                                                                                                                           tCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCg        <  2:201593/149‑1 (MQ=255)
                                                                                                                                             gACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCgg       <  4:308914/148‑1 (MQ=255)
                                                                                                                                              aCAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa     <  4:103944/149‑1 (MQ=255)
                                                                                                                                               cAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAg    <  7:268873/149‑1 (MQ=255)
                                                                                                                                               cAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAg    <  4:252240/149‑1 (MQ=255)
                                                                                                                                                 gCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGcc  <  6:171529/149‑1 (MQ=255)
                                                                                                                                                  cTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGcc  >  1:93377/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCC  >  NC_002947/2258374‑2258666

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: