Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,702,835 G→A 5.5% L12L (CTG→TTG)  dppF ← dipeptide/heme ABC transporter ATPas

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,702,8350GA 5.5% 245.4 / 3.5 73L12L (CTG→TTG) dppFdipeptide/heme ABC transporter ATPas
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base A (2/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

TACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGC  >  NC_000913/3702688‑3702983
                                                                                                                                                   |                                                                                                                                                    
tACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGc                                                                                                                                                     >  4:213618/1‑149 (MQ=255)
tACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGc                                                                                                                                                     <  2:192060/149‑1 (MQ=255)
 aCCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCg                                                                                                                                                    <  1:519554/149‑1 (MQ=255)
    gCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                                                                 <  1:164514/149‑1 (MQ=255)
    gCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                                                                 <  3:316571/149‑1 (MQ=255)
       gtcAGATTCGCCCACTACTGCCAGCGTTTTTCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgtt                                                                                                                                              <  3:506226/147‑1 (MQ=255)
        ccAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCg                                                                                                                                                    >  1:156231/1‑142 (MQ=255)
        ccAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCg                                                                                                                                                    <  2:156231/142‑1 (MQ=255)
         cAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGc                                                                                                                                            <  3:65627/149‑1 (MQ=255)
            aTTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                                                                         <  3:328657/149‑1 (MQ=255)
             ttCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAggg                                                                                                                                        <  2:188484/149‑1 (MQ=255)
              tCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGt                                                                                                                                       <  1:186854/149‑1 (MQ=255)
              tCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGt                                                                                                                                       <  4:163227/149‑1 (MQ=255)
                 cccACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGc                                                                                                                                    <  3:117493/149‑1 (MQ=255)
                     ctactGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCctct                                                                                                                                >  1:66123/1‑149 (MQ=255)
                      tactGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCtt                                                                                                                               >  2:609202/1‑149 (MQ=255)
                         tGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGca                                                                                                                            >  1:117841/1‑148 (MQ=255)
                            cAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                                                                         <  2:66123/149‑1 (MQ=255)
                                 ttttGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATa                                                                                                                    <  2:191467/149‑1 (MQ=255)
                                 ttttGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATa                                                                                                                    >  2:389503/1‑149 (MQ=255)
                                   ttGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgt                                                                                                                  <  2:138691/149‑1 (MQ=255)
                                      ccACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTc                                                                                                               <  2:285027/149‑1 (MQ=255)
                                       cACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCg                                                                                                              >  3:236864/1‑149 (MQ=255)
                                       cACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCg                                                                                                              >  3:429323/1‑149 (MQ=255)
                                       cACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCg                                                                                                              >  4:569041/1‑149 (MQ=255)
                                             cAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTcc                                                                                                        <  2:462486/149‑1 (MQ=255)
                                                 gTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGc                                                                                                    >  4:107680/1‑149 (MQ=255)
                                                   tAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcat                                                                                                  <  2:347656/149‑1 (MQ=255)
                                                    aaacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatc                                                                                                 >  2:390358/1‑149 (MQ=255)
                                                     aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatca                                                                                                <  3:152828/149‑1 (MQ=255)
                                                           acgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGt                                                                                          <  1:390358/149‑1 (MQ=255)
                                                           acgCCATCCAGCGCTTTAACCAAACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCTGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGt                                                                                          <  2:578362/149‑1 (MQ=255)
                                                               cATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGt                                                                                      >  1:3016/1‑149 (MQ=255)
                                                                 tCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                                                                 <  1:523719/88‑1 (MQ=255)
                                                                 tCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                                                                 >  2:523719/1‑88 (MQ=255)
                                                                    aGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATgg                                                                                 >  2:135264/1‑149 (MQ=255)
                                                                    aGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATgg                                                                                 >  1:452474/1‑149 (MQ=255)
                                                                        cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCAtt                                                                             <  2:485192/149‑1 (MQ=255)
                                                                        cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCAtt                                                                             <  1:389503/149‑1 (MQ=255)
                                                                         tttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCAttt                                                                            <  3:480863/149‑1 (MQ=255)
                                                                         tttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCAttt                                                                            <  1:212449/149‑1 (MQ=255)
                                                                          ttAACCAGACGTTCCGGTGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                                                                           <  3:498386/149‑1 (MQ=255)
                                                                          ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                                                                           <  2:437079/149‑1 (MQ=255)
                                                                            aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTgg                                                                                        >  3:616633/1‑134 (MQ=255)
                                                                            aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTgg                                                                                        <  4:616633/134‑1 (MQ=255)
                                                                            aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGa                                                                         >  4:196160/1‑149 (MQ=255)
                                                                               cAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgtt                                                                                                                                              >  2:454026/1‑77 (MQ=255)
                                                                               cAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgtt                                                                                                                                              <  1:454026/77‑1 (MQ=255)
                                                                                    gTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGcc                                                                 >  4:365910/1‑149 (MQ=255)
                                                                                             cgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGa                                                        <  1:318710/149‑1 (MQ=255)
                                                                                               aaCATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGc                                                      <  1:403146/149‑1 (MQ=255)
                                                                                               aaCATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGc                                                      >  1:234150/1‑149 (MQ=255)
                                                                                                     ccTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTc                                                <  3:272616/149‑1 (MQ=255)
                                                                                                        ttcttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAgc                                             <  1:620660/149‑1 (MQ=255)
                                                                                                                cGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGttc                                     <  1:609202/149‑1 (MQ=255)
                                                                                                                 ggATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGttct                                    >  1:140761/1‑149 (MQ=255)
                                                                                                                  gATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGttctt                                   >  1:447465/1‑149 (MQ=255)
                                                                                                                   aTAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGttcttc                                  >  4:94445/1‑149 (MQ=255)
                                                                                                                    tAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCa                                 >  4:50194/1‑149 (MQ=255)
                                                                                                                         ttttttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGcgcg                            >  1:735/1‑149 (MQ=255)
                                                                                                                          tttttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGa                           <  3:107680/149‑1 (MQ=255)
                                                                                                                                   tCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCa                  >  2:301990/1‑149 (MQ=255)
                                                                                                                                       tAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTgg              >  3:438929/1‑149 (MQ=255)
                                                                                                                                        acccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGc             >  4:304953/3‑149 (MQ=255)
                                                                                                                                          ccTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCAt           >  3:398490/1‑149 (MQ=255)
                                                                                                                                           cTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTgg                                                                                        <  2:186151/71‑1 (MQ=38)
                                                                                                                                           cTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTgg                                                                                        >  1:186151/1‑71 (MQ=38)
                                                                                                                                            tGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAg         >  1:508010/1‑149 (MQ=255)
                                                                                                                                            tGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAg         >  4:261748/1‑149 (MQ=255)
                                                                                                                                            tGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAg         >  4:403960/1‑149 (MQ=255)
                                                                                                                                               aaCAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAgggg      >  2:81237/1‑149 (MQ=255)
                                                                                                                                                aCAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGc     >  2:191588/1‑149 (MQ=255)
                                                                                                                                                aCAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGc     >  3:502990/1‑149 (MQ=255)
                                                                                                                                                   gCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAgc  >  1:523405/1‑149 (MQ=255)
                                                                                                                                                   gCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAgc  >  1:301221/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGC  >  NC_000913/3702688‑3702983

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: