Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,026,547 A→C 7.6% K473N (AAA→AAC) ‡ PP_5276 → phospholipase D family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,026,5470AC 7.6% 202.5 / 7.9 66K473N (AAA→AAC) ‡PP_5276phospholipase D family protein
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base C (3/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GACCTTCATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAGCGGCGATTCCGGG  >  NC_002947/6026400‑6026695
                                                                                                                                                   |                                                                                                                                                    
gACCTTCATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAAc                                                                                                                                                     >  7:303552/1‑149 (MQ=255)
    ttCATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGa                                                                                                                                                  >  8:156368/1‑148 (MQ=255)
      cATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGa                                                                                                                                                  <  7:156368/146‑1 (MQ=255)
      cATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGAt                                                                                                                                                 <  8:272153/147‑1 (MQ=255)
      cATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATa                                                                                                                                                <  2:183532/148‑1 (MQ=255)
      cATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATa                                                                                                                                                <  3:101978/148‑1 (MQ=255)
      cATAGGCTCGTTGAACTTCGACCCGCGCGCAGTATTGTGGATCACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATa                                                                                                                                                <  8:321498/148‑1 (MQ=255)
       aTCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAg                                                                                                                                               <  3:97781/148‑1 (MQ=255)
       aTCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGCGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAgg                                                                                                                                              <  1:291254/149‑1 (MQ=255)
          atcTCGTTCTACTTCGACCCGGGCTGAGTATTGTGGAACACCGAGGTCGGTGTACTGTTGGACAGCCCGGAGCTGGCCGATTACACCCGCGAACTGGCCCAGTAGGGTATTGCGCCTGCTCTTAGCTATCAGGTGAAACTGATAGGTaa                                                                                                                                           <  6:64336/147‑1 (MQ=255)
           gCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTaa                                                                                                                                           >  4:95463/1‑148 (MQ=255)
           gCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAAc                                                                                                                                          <  8:158253/149‑1 (MQ=255)
            cTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGACGCCTGCGCTCAGCTATCAGATGAAACTGATAGGTAAc                                                                                                                                          >  8:147423/1‑148 (MQ=255)
             tCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAAc                                                                                                                                          >  6:11572/1‑147 (MQ=255)
               gctCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGc                                                                                                                                      <  6:337580/147‑1 (MQ=255)
               gTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGc                                                                                                                                      >  7:231770/1‑149 (MQ=255)
               gTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGACGAATACACCCGCGAACTGGCCCAGTAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACagat                                                                                                                                      >  5:64336/1‑146 (MQ=255)
                  cAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACACGCGAACTGGCCCAGCATGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCt                                                                                                                                     <  2:256254/147‑1 (MQ=255)
                        cGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTgg                                                                                                                              <  3:95463/148‑1 (MQ=255)
                          aCCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGc                                                                                                                            >  8:274794/1‑148 (MQ=255)
                          aCCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGc                                                                                                                            >  3:304135/1‑148 (MQ=255)
                           cccGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGa                                                                                                                          <  6:351126/149‑1 (MQ=255)
                               cgcTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGa                                                                                                                      <  3:174899/149‑1 (MQ=255)
                                gcTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAg                                                                                                                     <  7:274794/149‑1 (MQ=255)
                                    aGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGa                                                                                                                   <  4:304135/147‑1 (MQ=255)
                                             gAACACCGAGGTCGGGGTACTGGTGGACAGCCTGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTg                                                                                                                                    <  6:17988/121‑1 (MQ=255)
                                             gAACACCGAGGTCGGGGTACTGGTGGACAGCCAGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTg                                                                                                                                    >  5:17988/1‑121 (MQ=255)
                                                       ggaggtgTACTGGTGGACAGCCCTGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGc                                                                                               <  8:303552/142‑1 (MQ=255)
                                                           gggTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAgg                                                                                                                    <  1:266216/123‑1 (MQ=255)
                                                           gggTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAgg                                                                                                                    >  2:266216/1‑123 (MQ=255)
                                                            ggTACTGGTGGACAGCCCTGAGCTGTCCGAATACACCCGCGATCTGGCCCATCAGGGTATGGTTCATGCGCTGAGCTATCAGGTGAAACTGATAGTTAACAAGCTGGTGTTGGAGACCGAGGACAATGGCCAAAGGCAAGTGCTGACGt                                                                                         <  8:231770/149‑1 (MQ=255)
                                                                   gtggACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGc                                                                                                            <  5:346183/123‑1 (MQ=255)
                                                                   gtggACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGc                                                                                                            >  6:346183/1‑123 (MQ=255)
                                                                   gtggACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGcc                                                                                  >  8:319095/1‑149 (MQ=255)
                                                                                     gCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACg                                                                >  5:150271/1‑149 (MQ=255)
                                                                                        gAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCtt                                                             <  7:147423/149‑1 (MQ=255)
                                                                                         aaTACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTc                                                            <  1:38239/149‑1 (MQ=255)
                                                                                                 cgcgAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCcgggc                                                                              >  7:119195/1‑123 (MQ=255)
                                                                                                 cgcgAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCcgggc                                                                              <  8:119195/123‑1 (MQ=255)
                                                                                                     aaCTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGAt                                                 >  6:16331/1‑148 (MQ=255)
                                                                                                        tGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCcgggcggg                                                                           >  6:340076/1‑119 (MQ=255)
                                                                                                        tGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCcgggcggg                                                                           <  5:340076/119‑1 (MQ=255)
                                                                                                          gCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGc                                             >  7:267481/1‑147 (MQ=255)
                                                                                                            ccAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATc                                                >  3:46987/1‑142 (MQ=255)
                                                                                                            ccAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATc                                                <  4:46987/142‑1 (MQ=255)
                                                                                                            ccAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCa                                               <  6:62918/143‑1 (MQ=255)
                                                                                                            ccAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCa                                               >  5:62918/1‑143 (MQ=255)
                                                                                                                   ggTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTgg                                  >  7:160135/1‑149 (MQ=255)
                                                                                                                        ggCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTg                                                                                                                                    >  4:190643/1‑46 (MQ=255)
                                                                                                                        ggCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTg                                                                                                                                    <  3:190643/46‑1 (MQ=255)
                                                                                                                         gcgcCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGga                            <  5:11572/149‑1 (MQ=255)
                                                                                                                           gcCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGgag                           >  4:300439/1‑148 (MQ=255)
                                                                                                                                cgcTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGCGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATg                     >  7:216481/1‑149 (MQ=255)
                                                                                                                                   ttaccTATCAGTTTCACCTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGga                            >  3:275272/5‑139 (MQ=255)
                                                                                                                                   ttaccTATCAGTTTCACCTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGga                            <  4:275272/135‑1 (MQ=255)
                                                                                                                                      ccTATCAGTTTCACCTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGga                            <  3:28134/135‑1 (MQ=255)
                                                                                                                                      ccTATCAGTTTCACCTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGga                            >  4:28134/2‑136 (MQ=255)
                                                                                                                                      ccTATCAGTTTCACCTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAg               >  5:244424/2‑149 (MQ=255)
                                                                                                                                       cTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGGCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAg               <  6:150271/148‑1 (MQ=255)
                                                                                                                                         aTCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGagaaga                       >  1:300628/1‑138 (MQ=255)
                                                                                                                                         aTCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGagaaga                       <  2:300628/138‑1 (MQ=255)
                                                                                                                                          tCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAgcggc           >  7:89855/1‑149 (MQ=255)
                                                                                                                                          tCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAgcggc           >  7:101708/1‑149 (MQ=255)
                                                                                                                                              gTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTc                                                            <  8:204341/96‑1 (MQ=255)
                                                                                                                                              gTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTc                                                            >  7:204341/1‑96 (MQ=255)
                                                                                                                                              gTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAGCGGCGAtt       >  1:5799/1‑149 (MQ=255)
                                                                                                                                                   aCTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAGCGGCGATTCCggg  <  8:160135/149‑1 (MQ=255)
                                                                                                                                                   aCTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAGCGGCGATTCCggg  >  7:8016/1‑149 (MQ=255)
                                                                                                                                                   aCTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAGCGGCGATTCCggg  <  2:5799/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACCTTCATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAGCGGCGATTCCGGG  >  NC_002947/6026400‑6026695

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: ‑