Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,258,516 Δ1 bp 8.3% coding (180/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,5160C. 8.3% 227.4 / 14.1 72coding (180/1005 nt)asdaspartate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (31/35);  new base . (4/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

CCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAG  >  NC_002947/2258374‑2258664
                                                                                                                                               |                                                                                                                                                      
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                                                                                                                        <  3:50060/132‑1 (MQ=255)
                                          tGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCct                                                                                                                        >  4:50060/1‑132 (MQ=255)
                                             cATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCgccgcc                                                                                                    <  4:295086/149‑1 (MQ=255)
                                             cATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCgccgcc                                                                                                    <  8:23895/149‑1 (MQ=255)
                                               ttCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGt                                                                                                  >  8:102812/1‑149 (MQ=255)
                                                         gCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAg                                                                                                                                                      <  4:306238/89‑1 (MQ=255)
                                                         gCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAg                                                                                                                                                      >  3:306238/1‑89 (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                                                                                                                                          >  4:206978/1‑88 (MQ=39)
                                                                     ctgctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAAGCTGTCGACTTCg                                                                                                                                          <  3:206978/88‑1 (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:17764/149‑1 (MQ=255)
                                                                             cAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGc                                                                     <  8:243385/149‑1 (MQ=255)
                                                                             cAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGc                                                                     >  6:9936/1‑149 (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                                                           >  3:245387/1‑140 (MQ=255)
                                                                                               cagcagCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCt                                                           <  4:245387/140‑1 (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)
                                                                                                                                               |                                                                                                                                                      
CCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAG  >  NC_002947/2258374‑2258664

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: