Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,281,586 A→T 7.6% intergenic (‑163/+555) PP_3751 ← / ← PP_3753 cyanate transporter/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,281,5860AT 7.6% 192.1 / 5.5 66intergenic (‑163/+555)PP_3751/PP_3753cyanate transporter/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (26/35);  new base T (3/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

TTGTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCTTGAACACTCGCTGCTGGAGC  >  NC_002947/4281439‑4281732
                                                                                                                                                   |                                                                                                                                                  
ttGTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAg                                                                                                                                                   >  1:186957/1‑149 (MQ=255)
 tGTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAg                                                                                                                                                   >  2:17419/1‑148 (MQ=255)
  gTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAg                                                                                                                                                   <  5:319128/147‑1 (MQ=255)
  gTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAg                                                                                                                                                   >  5:134721/1‑147 (MQ=255)
   tAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTa                                                                                                                                                 <  2:186957/148‑1 (MQ=255)
   tAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTCGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTa                                                                                                                                                 <  3:32769/148‑1 (MQ=255)
    aaCAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTCCAGTAgc                                                                                                                                               >  7:191344/1‑149 (MQ=255)
     aCAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAgcg                                                                                                                                              >  5:352597/1‑149 (MQ=255)
     aCAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATCGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAgc                                                                                                                                               <  5:109435/148‑1 (MQ=255)
       aGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAgcgc                                                                                                                                             >  6:129670/1‑148 (MQ=255)
       aGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAgcgc                                                                                                                                             >  7:126314/1‑148 (MQ=255)
       aGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGcc                                                                                                                                            <  5:247392/149‑1 (MQ=255)
          aaCTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACTGTAGCGccgct                                                                                                                                         <  8:317415/149‑4 (MQ=255)
          aaCTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACTGTAGCGccgc                                                                                                                                          <  7:136089/148‑3 (MQ=255)
          aaCTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACTGTAGCGccgc                                                                                                                                          >  8:136089/1‑146 (MQ=255)
           aCTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCatat                                                                                                                                        <  3:192402/149‑1 (MQ=255)
              ggCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCa                                                                                                                                     >  5:268365/1‑149 (MQ=255)
               gCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCAc                                                                                                                                    >  1:21648/1‑149 (MQ=255)
                 cATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCAccc                                                                                                                                  <  2:358234/149‑1 (MQ=255)
                  aTTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCAccc                                                                                                                                  >  6:282098/1‑148 (MQ=255)
                      tCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCcaca                                                                                                                              <  7:56978/148‑1 (MQ=255)
                        ccGGCTCGGCTGGCCTGGCAGAAATGAAGAGTTTGGATGCAGGTCAGTGGCCCCCTGAATGGACGATTGGCCTCCCGCGTATGCGCAAAATTGGCTAGCGGGGGGGCCAATAAGCGGCGCGACAGTAGCGCCCTATCccccagagcgg                                                                                                                            >  7:390230/1‑141 (MQ=255)
                           gcgcggctGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTg                                                                                                                        <  8:126314/146‑1 (MQ=255)
                           cctcgcctGGCCTGGCAGAACTGAACACAGTGGATCCAGGCCAGGGGCCACCTGAATGGACGATCGGCCCCCCGCTTATGCGCAAAATTGGCCAGCGGGGCGGCCAATAAGCGGCGCTACAGCAGCGCCATATCCACCCCACACGCTTg                                                                                                                        <  8:191344/148‑1 (MQ=255)
                             tcggctGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGga                                                                                                                      <  4:90006/149‑1 (MQ=255)
                               ggctGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGga                                                                                                                      <  5:260431/147‑1 (MQ=255)
                                         gcAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATg                                                                                                            <  2:21648/147‑1 (MQ=255)
                                         gcAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGca                                                                                                          <  2:45948/149‑1 (MQ=255)
                                         gcAGAAATGAACAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGACTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGc                                                                                                           <  8:366531/148‑1 (MQ=255)
                                          cAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGcac                                                                                                         <  1:435491/149‑1 (MQ=255)
                                                   agagATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATc                                                                                                <  6:352597/149‑1 (MQ=255)
                                                   agagATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATc                                                                                                <  1:281500/149‑1 (MQ=255)
                                                      gaTTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGaa                                                                                             <  7:81477/149‑1 (MQ=255)
                                                      gaTTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGaa                                                                                             >  5:103387/1‑149 (MQ=255)
                                                      gaTTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGaa                                                                                             >  3:309691/1‑149 (MQ=255)
                                                       aTTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAAt                                                                                            >  6:196316/1‑149 (MQ=255)
                                                           gATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGc                                                                                        >  1:367719/1‑149 (MQ=255)
                                                            aTGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGc                                                                                        <  5:196316/148‑1 (MQ=255)
                                                               cAGGTCACTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCaa                                                                                     <  6:268365/148‑1 (MQ=255)
                                                                aGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAac                                                                                   <  3:229136/149‑1 (MQ=255)
                                                                            aCCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAg                                                                       >  7:363241/1‑149 (MQ=255)
                                                                              cTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGcc                                                                     >  3:61692/1‑149 (MQ=255)
                                                                                   tGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATc                                                                 >  2:23840/1‑148 (MQ=255)
                                                                                     gACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAg                                                              >  1:198960/1‑149 (MQ=255)
                                                                                       cGAGTGGCCTCCCTCTTTTGAGCAAAATTGGCTATCGGGGTGTCCAATATGCGGCGCTACAGTAGCGCCATATCCACCAGACTCGCTTGGAGACTTACATGCACAACCGTATCGAATGGGCCAAACACGCGGCGTAAGCCTATATAGcc                                                            <  2:340570/149‑1 (MQ=255)
                                                                                        gATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACCCTTCGAGACTGACATGCACAACCGTATCGAATCCGCCAAACACGCGCCGTAAGCCTATCAATCCa                                                           >  8:97610/1‑149 (MQ=255)
                                                                                          ttGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATg                                                         >  5:309325/1‑149 (MQ=255)
                                                                                                    gCTTCTGCGCAAAATTGGCTCTCTGGGCGGGCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTgg                                               <  8:324290/149‑1 (MQ=255)
                                                                                                     cTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGa                                              <  4:61692/149‑1 (MQ=255)
                                                                                                        aTGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGc                                            <  5:165216/148‑1 (MQ=255)
                                                                                                           cgcAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGgc                                        <  2:198960/149‑1 (MQ=255)
                                                                                                           cgcAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGgc                                        >  5:288917/1‑149 (MQ=255)
                                                                                                           cgcAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGgc                                        >  7:259219/1‑149 (MQ=255)
                                                                                                               aaaTTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTg                                    <  6:134721/149‑1 (MQ=255)
                                                                                                                aaTTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACAGGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTgg                                   <  7:97610/149‑1 (MQ=255)
                                                                                                                aaTTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTgg                                   <  5:282098/149‑1 (MQ=255)
                                                                                                                        aTCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTTTCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGc                           >  5:334133/1‑149 (MQ=255)
                                                                                                                              gTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCtt                     <  2:367719/149‑1 (MQ=255)
                                                                                                                                 gCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCTTGaa                  <  6:103387/149‑1 (MQ=255)
                                                                                                                                    aaTAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCTTGAacac               >  6:176647/1‑149 (MQ=255)
                                                                                                                                      gatatggcgCTACTGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCTTGAACACTc             >  3:298007/6‑149 (MQ=255)
                                                                                                                                             cgcTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCTTGAACACTCGctgctg      <  8:363241/149‑1 (MQ=255)
                                                                                                                                              gcTACTGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCTTGAACACTCGctgctg      >  3:200771/1‑148 (MQ=255)
                                                                                                                                               cTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTCTCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCTTGAACACTCGCTGCTGGa    <  8:189030/149‑1 (MQ=255)
                                                                                                                                                 aCAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCTTGAACACTCGCTGCTGGAGc  >  7:27238/1‑149 (MQ=255)
                                                                                                                                                  cAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCTTGAACACTCGCTGCTGGAg   <  8:259219/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGCCACCTGAATGGACGATTGGCCTCCCGCTTATGCGCAAAATTGGCTATCGGGGTGGCCAATAAGCGGCGCTACAGTAGCGCCATATCCACCCCACACGCTTGGAGACTGACATGCACAACCGTATCGAATGGGCCAAACACGCGCCGTAAGCCTATCAAGCCATGGTCGGCCTGGAGCAGGCGCTGGCCAACAGCGGCCTTGAACACTCGCTGCTGGAGC  >  NC_002947/4281439‑4281732

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: