Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,340,372 A→G 10.3% intergenic (+83/+45) ubiG → / ← yfaL bifunctional 3‑demethylubiquinone‑9 3‑methyltransferase/ 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,340,3720AG10.3% 264.1 / 13.6 85intergenic (+83/+45)ubiG/yfaLbifunctional 3‑demethylubiquinone‑9 3‑methyltransferase/ 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin
Reads supporting (aligned to +/- strand):  ref base A (40/36);  new base G (4/5);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATACTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTG  >  NC_000913/2340226‑2340516
                                                                                                                                                  |                                                                                                                                                
tAATACTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAg                                                                                                                                                >  2:10749/1‑149 (MQ=255)
tAATACTTTTAAACTCGGCCCCGGCGTGAATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAg                                                                                                                                                <  3:841483/149‑1 (MQ=255)
 aaTACTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAgg                                                                                                                                               <  2:416519/149‑1 (MQ=255)
    aCTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGt                                                                                                                                            <  4:489131/149‑1 (MQ=255)
     cTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGtt                                                                                                                                           <  3:640999/149‑1 (MQ=255)
      ttttAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGtaa                                                                                                                                          >  1:482860/1‑147 (MQ=255)
      ttttAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGttt                                                                                                                                          >  2:1103057/1‑149 (MQ=255)
      ttttAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGACTTATCAGGCCAACAAATGCTCCCCGTAGGCCGGATAAGGCGttt                                                                                                                                          >  1:1453949/1‑149 (MQ=255)
          aaaCTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATaa                                                                                                                                                 <  1:763603/138‑1 (MQ=255)
          aaaCTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATaa                                                                                                                                                 >  2:763603/1‑138 (MQ=255)
           aaCTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                     >  2:305205/1‑149 (MQ=255)
               cGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                 <  3:64953/149‑5 (MQ=38)
               cGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                 <  2:126555/149‑5 (MQ=38)
                ggCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                >  4:1850817/1‑149 (MQ=255)
                  ccccGTCGTGGATGTGAACTTTATGCTGCACACGCAGAATAAGTCAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTTGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                              <  2:1453949/149‑1 (MQ=255)
                   cccGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                             <  3:949244/149‑1 (MQ=38)
                    ccGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                            <  1:10749/149‑1 (MQ=255)
                              tgtgAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                 <  2:1433269/134‑5 (MQ=38)
                              tgtgAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                 >  1:1433269/1‑130 (MQ=38)
                                    cTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGa                                                                                                            <  1:1103057/149‑1 (MQ=255)
                                     tataTGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                      >  4:577350/1‑122 (MQ=255)
                                     tataTGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                      <  3:577350/122‑1 (MQ=255)
                                       tatGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                <  1:1407101/126‑6 (MQ=38)
                                       tatGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                >  2:1407101/1‑121 (MQ=38)
                                              cacaCGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAc                                                                                                  <  3:1850817/149‑1 (MQ=255)
                                              cacaCGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGGGATTAc                                                                                                  >  1:631179/1‑149 (MQ=255)
                                                cacGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATaa                                                                                                                                                 >  4:497609/1‑100 (MQ=255)
                                                cacGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATaa                                                                                                                                                 <  3:497609/100‑1 (MQ=255)
                                                 acGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAAcgcct                                                                                                                                    >  4:535912/1‑111 (MQ=37)
                                                  cGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACt                                                                                                               >  4:218512/1‑132 (MQ=255)
                                                  cGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACt                                                                                                               <  3:218512/132‑1 (MQ=255)
                                                   gCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt                                                                                             >  4:515986/1‑149 (MQ=255)
                                                       aaTAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCAcc                                                                                         <  2:631179/149‑1 (MQ=255)
                                                             tGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAt                                                                                               <  1:619164/137‑1 (MQ=255)
                                                             tGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAt                                                                                               >  2:619164/1‑137 (MQ=255)
                                                                     aTGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACaaaaa                                                                           >  3:47227/1‑149 (MQ=255)
                                                                       gTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTa                                                                                                              <  4:1862601/112‑1 (MQ=255)
                                                                       gTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTa                                                                                                              >  3:1862601/1‑112 (MQ=255)
                                                                       gTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAATcc                                                                         <  3:515986/149‑1 (MQ=255)
                                                                         ttGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATc                                                                                   >  3:1707093/1‑137 (MQ=255)
                                                                         ttGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATc                                                                                   <  4:1707093/137‑1 (MQ=255)
                                                                                  ccAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTgg                                                              <  1:305205/149‑1 (MQ=255)
                                                                                         cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt                                                                                             >  4:1666006/1‑111 (MQ=255)
                                                                                         cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt                                                                                             <  3:1666006/111‑1 (MQ=255)
                                                                                                 aCGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTc                                               >  1:1456991/1‑149 (MQ=255)
                                                                                                 aCGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTc                                               >  1:1457811/1‑149 (MQ=255)
                                                                                                   gCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTccc                                             >  4:378456/1‑149 (MQ=255)
                                                                                                        ccGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAAccc                                                                        <  1:167840/117‑1 (MQ=255)
                                                                                                        ccGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAAccc                                                                        >  2:167840/1‑117 (MQ=255)
                                                                                                        ccGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCt                                        >  1:1812827/1‑149 (MQ=255)
                                                                                                            gTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAAtcatca                                                  <  1:879270/135‑1 (MQ=255)
                                                                                                            gTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTCCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAAtcatca                                                  >  2:879270/1‑135 (MQ=255)
                                                                                                             tcgcaTCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCTTTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccac                                   <  2:1664582/145‑1 (MQ=255)
                                                                                                             tCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccac                                   <  2:1456991/149‑1 (MQ=255)
                                                                                                             tCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccac                                   >  1:85495/1‑149 (MQ=255)
                                                                                                             tCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccac                                   <  3:378456/149‑1 (MQ=255)
                                                                                                                  tCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccacgtcac                              >  1:1850257/1‑149 (MQ=255)
                                                                                                                  tCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccacgccac                              >  4:879288/1‑149 (MQ=255)
                                                                                                                    aGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAg                                                                                                         >  4:686703/1‑72 (MQ=255)
                                                                                                                    aGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAg                                                                                                         <  1:973611/72‑1 (MQ=255)
                                                                                                                    aGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAg                                                                                                         >  2:973611/1‑72 (MQ=255)
                                                                                                                    aGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAg                                                                                                         <  3:686703/72‑1 (MQ=255)
                                                                                                                    aGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCcacgccacgc                            >  4:1285310/1‑149 (MQ=255)
                                                                                                                    aGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCcacgccacgc                            >  1:1074132/1‑149 (MQ=255)
                                                                                                                     ggCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCt                           >  2:1435371/1‑149 (MQ=255)
                                                                                                                       ccTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACa                                                                                                          >  2:1861140/1‑68 (MQ=255)
                                                                                                                       ccTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACa                                                                                                          <  1:1861140/68‑1 (MQ=255)
                                                                                                                         tACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAt                                                                                               <  3:1695451/77‑1 (MQ=255)
                                                                                                                         tACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAt                                                                                               >  4:1695451/1‑77 (MQ=255)
                                                                                                                            aaaTGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGc                                                                     <  1:1810109/100‑1 (MQ=255)
                                                                                                                            aaaTGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGc                                                                     >  2:1810109/1‑100 (MQ=255)
                                                                                                                            aaaTGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccac                                   <  4:547812/134‑1 (MQ=255)
                                                                                                                            aaaTGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccac                                   >  3:547812/1‑134 (MQ=255)
                                                                                                                               tGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAgg                 <  2:85495/149‑1 (MQ=255)
                                                                                                                                gCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCt                                                                <  3:1703190/101‑1 (MQ=255)
                                                                                                                                gCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCt                                                                >  4:1703190/1‑101 (MQ=255)
                                                                                                                                  tCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAatc                                                      >  2:1360976/1‑109 (MQ=255)
                                                                                                                                  tCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAatc                                                      <  1:1360976/109‑1 (MQ=255)
                                                                                                                                  tCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCgcgc                    >  4:399077/1‑143 (MQ=255)
                                                                                                                                  tCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCgcgc                    <  3:399077/143‑1 (MQ=255)
                                                                                                                                  gaaccGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGac              >  2:102005/4‑149 (MQ=255)
                                                                                                                                        tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGTCACGCTGCCGCGCAGGGACACTCGc        <  2:1850257/149‑1 (MQ=255)
                                                                                                                                            ccGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGcc                                                                    >  1:288285/1‑85 (MQ=255)
                                                                                                                                            ccGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGcc                                                                    <  2:288285/85‑1 (MQ=255)
                                                                                                                                            ccGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgac    <  2:1457811/149‑1 (MQ=255)
                                                                                                                                            ccGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgac    >  3:614714/1‑149 (MQ=255)
                                                                                                                                              ggATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgactg  >  3:5980/1‑149 (MQ=255)
                                                                                                                                                aTAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgac    >  2:1005080/1‑145 (MQ=255)
                                                                                                                                                aTAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgac    <  1:1005080/145‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TAATACTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTG  >  NC_000913/2340226‑2340516

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