Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,267,957 T→A 6.6% L143Q (CTG→CAG)  azoR1 → FMN‑dependent NADH‑azoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,267,9570TA 6.6% 176.9 / 2.4 61L143Q (CTG→CAG) azoR1FMN‑dependent NADH‑azoreductase
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base A (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

CCCATGTACAACTTCACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCA  >  NC_002947/3267809‑3268102
                                                                                                                                                    |                                                                                                                                                 
cccATGTACAACTTCACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCt                                                                                                                                                   <  5:245143/149‑1 (MQ=255)
      tACAACTTCACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGc                                                                                                                                             >  5:25158/1‑149 (MQ=255)
       aCAACTTCACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCt                                                                                                                                            >  6:64570/1‑149 (MQ=255)
       aCAACTTCACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCt                                                                                                                                            >  8:251372/1‑149 (MQ=255)
        cAACTTCACCGTGCCGACCCAGCTCATAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTg                                                                                                                                           <  4:71778/149‑1 (MQ=255)
         aaCTTCACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTgg                                                                                                                                          >  1:203238/1‑149 (MQ=255)
         aaCTTCACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTgg                                                                                                                                          >  6:44032/1‑149 (MQ=255)
            ttCACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGc                                                                                                                                         >  8:239206/1‑147 (MQ=255)
              cACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAg                                                                                                                                      <  3:310945/148‑1 (MQ=255)
              cACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGc                                                                                                                                     >  6:89460/1‑149 (MQ=255)
               aCCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGcc                                                                                                                                    <  1:146144/149‑1 (MQ=255)
                     ccGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTgg                                                                                                                              <  8:320616/149‑1 (MQ=255)
                     ccGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGCTGGCCAGCGTGCAGGCCACCCgcagtg                                                                                                                              <  6:27552/149‑6 (MQ=255)
                       gACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGCGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGgcg                                                                                                                            >  5:220235/1‑149 (MQ=255)
                          ccAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCg                                                                                                                         >  7:255911/1‑149 (MQ=255)
                          ccAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCg                                                                                                                         <  2:162854/149‑1 (MQ=255)
                          ccAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCg                                                                                                                         >  6:209335/1‑149 (MQ=255)
                                cAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCAGCGTGCAGGCCac                                                                                                                                      <  6:103969/130‑2 (MQ=255)
                                cAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCAGCGTGCAGGCCac                                                                                                                                      >  5:103969/1‑129 (MQ=255)
                                 aaaGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATg                                                                                                                   >  8:139734/1‑148 (MQ=255)
                                 aaaGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATg                                                                                                                   >  5:35388/1‑148 (MQ=255)
                                  aaGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGaa                                                                                                                 <  1:337784/149‑1 (MQ=255)
                                     ccTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGAc                                                                                                              >  5:346238/1‑149 (MQ=255)
                                         gATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCTAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACtt                                                                                                            <  6:66674/147‑1 (MQ=255)
                                           ttGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACt                                                                                                                                >  2:112980/1‑125 (MQ=255)
                                           ttGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACt                                                                                                                                <  1:112980/125‑1 (MQ=255)
                                           ttGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCt                                                                                                         >  4:72404/1‑148 (MQ=255)
                                              aCCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGaaa                                                                                                     <  6:35388/149‑1 (MQ=255)
                                              aCCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGa                                                                                                       >  1:249156/1‑147 (MQ=255)
                                              aCCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGa                                                                                                       >  4:174416/1‑147 (MQ=255)
                                                 gcgTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAgtg                                                                                                  <  5:89460/149‑1 (MQ=255)
                                                        cGTCGCCGGCAAGACCTACCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGCTGCACGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGCCATGAAGACTACCTGAAAGTGTTCCTgg                                                                                           >  5:46334/1‑149 (MQ=255)
                                                                   aGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCggg                                                                                <  7:251372/149‑1 (MQ=255)
                                                                    gACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGa                                                                               <  3:72404/149‑1 (MQ=255)
                                                                       cTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTa                                                                            >  3:230633/1‑149 (MQ=255)
                                                                         tCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACt                                                                          >  5:26432/1‑149 (MQ=255)
                                                                          ccGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACt                                                                          <  6:26432/148‑1 (MQ=255)
                                                                                cacCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTgg                                                                   <  5:209335/149‑1 (MQ=255)
                                                                                cacCGAAGCCGGCCCGCAAGGCCTGCGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCt                                                                     <  6:220235/147‑1 (MQ=255)
                                                                                   cGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCGGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGAGCCGGGCTTGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGaa                                                                 >  5:363903/1‑148 (MQ=255)
                                                                                   cGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGGGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCATGAAAGTGTTGCTGGGCTTCATCGGGATTACTGACCTGGaaa                                                                >  7:254823/1‑149 (MQ=255)
                                                                                    gAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAAt                                                               >  7:286906/1‑149 (MQ=255)
                                                                                     aaGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATc                                                              <  8:262668/149‑1 (MQ=255)
                                                                                        ccGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTc                                                           >  8:269424/1‑149 (MQ=255)
                                                                                              cGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGcc                                                     <  6:346238/149‑1 (MQ=255)
                                                                                                cAAGCCCTGGGCGGCAACAAGAAAGCGGGGCTGGGTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCa                                                   <  8:255911/149‑1 (MQ=255)
                                                                                                      cTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACggcct                                             <  4:230633/149‑1 (MQ=255)
                                                                                                          gCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACggcctggcc                                         <  5:64570/149‑1 (MQ=255)
                                                                                                               aacaaGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATgg                                    >  8:61409/1‑149 (MQ=255)
                                                                                                                   aGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATTAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCg                                <  7:239206/149‑1 (MQ=255)
                                                                                                                     aaaGTGGTGCTGGTTTCCACTGCGGGGGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCATAGGCCAGGCCGTGGGCCCGGa                              <  2:192430/149‑1 (MQ=255)
                                                                                                                        gtggtgCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGagca                           <  2:203238/149‑1 (MQ=255)
                                                                                                                        gtggtgCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGagca                           <  6:25158/149‑1 (MQ=255)
                                                                                                                            tgCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAgcagcg                        >  7:135141/1‑148 (MQ=255)
                                                                                                                             gCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTCGCCTATGGCCCGGAgcagcg                        >  4:352773/1‑147 (MQ=255)
                                                                                                                                  tttCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGGTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAg                 <  5:44032/149‑1 (MQ=255)
                                                                                                                                       aCTGCTGGTGGCCTGCCCTCTGGCCAGCCGACTGGCGCCGGGCATGACGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTTGAAATCGTCCGTGCCCACTGCCTGGCCTATGGCCCGGAGCAGCGCAGCCACGcgat            <  6:363903/149‑1 (MQ=255)
                                                                                                                                       aCTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGcgat            >  4:33545/1‑149 (MQ=255)
                                                                                                                                         tGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGcgatc           >  2:198225/1‑148 (MQ=255)
                                                                                                                                           cGGGTGGCCTGCACGCTGGCGAGCCTACTGGCGGAGTGCATGAAGCCTTCCTGATAGTGGTCCTGGTCTTCATCGGGTTTTCTGACCTGGAAATCGTCCGGGCCCACGGGCTGGCCTATGGCGCGGAGGTGCGCAGCCAGGcgatcgat        <  6:46334/149‑1 (MQ=255)
                                                                                                                                              ctGGCCAGCGTGCAGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGAtgct     >  2:4750/2‑149 (MQ=255)
                                                                                                                                                   cTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGca  >  1:161729/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCCATGTACAACTTCACCGTGCCGACCCAGCTCAAAGCCTGGATTGACCGCGTAGCCGTCGCCGGCAAGACCTTCCGCTACACCGAAGCCGGCCCGCAAGGCCTGTGCGGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCA  >  NC_002947/3267809‑3268102

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑