Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,026,542 G→T 8.4% V472L (GTG→TTG) ‡ 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,5420GT 8.4% 197.4 / 7.2 65V472L (GTG→TTG) ‡PP_5276phospholipase D family protein
Reads supporting (aligned to +/- strand):  ref base G (27/32);  new base T (3/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

AAGACCTTCATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAGCGGCGATT  >  NC_002947/6026398‑6026690
                                                                                                                                                |                                                                                                                                                    
aaGACCTTCATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGaa                                                                                                                                                  <  5:320773/149‑1 (MQ=255)
  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)
                                 cgcTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGa                                                                                                                                                   >  8:259099/1‑115 (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                                          >  5:62918/1‑143 (MQ=255)
                                                                                                              ccAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCa                                          <  6:62918/143‑1 (MQ=255)
                                                                                                                     ggTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTgg                             >  7:160135/1‑149 (MQ=255)
                                                                                                                          ggCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTg                                                                                                                               <  3:190643/46‑1 (MQ=255)
                                                                                                                          ggCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTg                                                                                                                               >  4:190643/1‑46 (MQ=255)
                                                                                                                           gcgcCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGga                       <  5:11572/149‑1 (MQ=255)
                                                                                                                             gcCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGgag                      >  4:300439/1‑148 (MQ=255)
                                                                                                                                  cgcTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGCGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATg                >  7:216481/1‑149 (MQ=255)
                                                                                                                                     ttaccTATCAGTTTCACCTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGga                       <  4:275272/135‑1 (MQ=255)
                                                                                                                                     ttaccTATCAGTTTCACCTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGga                       >  3:275272/5‑139 (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                  <  2:300628/138‑1 (MQ=255)
                                                                                                                                           aTCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGagaaga                  >  1:300628/1‑138 (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)
                                                                                                                                                |                                                                                                                                                    
AAGACCTTCATCGGCTCGTTCAACTTCGACCCGCGCTCAGTATTGTGGAACACCGAGGTCGGGGTACTGGTGGACAGCCCGGAGCTGGCCGAATACACCCGCGAACTGGCCCAGCAGGGTATGGCGCCTGCGCTGAGCTATCAGGTGAAACTGATAGGTAACAAGCTGGTGTGGGCGACCGAGGACAATGGCCAACGGCACGTGCTGACGTCTGAGCCGGGCGGGATCTGGCGACGCTTCAATGCCTGGATCAGCAAGGCGGTTGGCTTGGAGAAGATGCTATAGCGGCGATT  >  NC_002947/6026398‑6026690

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