Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,213,981 G→T 6.1% P48Q (CCG→CAG)  yjaB ← putative N‑acetyltransferase YjaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,213,9810GT 6.1% 273.5 / 5.1 82P48Q (CCG→CAG) yjaBputative N‑acetyltransferase YjaB
Reads supporting (aligned to +/- strand):  ref base G (50/27);  new base T (2/3);  total (52/30)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

AGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAACCATTTTCTCT  >  NC_000913/4213835‑4214130
                                                                                                                                                  |                                                                                                                                                     
agagCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCTTCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCa                                                                                                                                                    <  2:1353179/150‑1 (MQ=255)
   gcgcATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACCTCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACCGCGCTTCCGGCagga                                                                                                                                                 >  1:866018/1‑150 (MQ=255)
   gcgcATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCagga                                                                                                                                                 <  1:823306/150‑1 (MQ=255)
      cATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagg                                                                                                                                              >  1:1046174/1‑150 (MQ=255)
        tGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAggaaggaa                                                                                                                                            >  2:1127506/1‑150 (MQ=255)
         gCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCTGCcggaagcgcc                                                                                                                                           <  1:177293/150‑5 (MQ=255)
         gCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaac                                                                                                                                           <  1:170832/150‑1 (MQ=255)
          cTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCGACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacg                                                                                                                                          <  2:533885/150‑1 (MQ=255)
           tCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCGTGCcggaagcgccgt                                                                                                                                         <  2:1240866/150‑7 (MQ=255)
           tCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacga                                                                                                                                         >  2:1228615/1‑150 (MQ=255)
           tCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacga                                                                                                                                         >  2:74938/1‑150 (MQ=255)
                          ccTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGc                                                                                                                          >  2:666785/1‑150 (MQ=255)
                            tACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTc                                                                                                                        >  2:1555571/1‑150 (MQ=255)
                              cgccgcAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCgg                                                                                                                      <  1:1416525/150‑1 (MQ=255)
                                ccgcAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTc                                                                                                                    >  1:319399/1‑150 (MQ=255)
                                 cgcAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTcc                                                                                                                   >  2:166475/1‑150 (MQ=255)
                                    aGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGAt                                                                                                                <  2:1477030/150‑1 (MQ=255)
                                        gcgcACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATActc                                                                                                            >  2:944950/1‑150 (MQ=255)
                                            acaTCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTctgct                                                                                                        <  1:1625402/150‑1 (MQ=255)
                                            acaTCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTctgct                                                                                                        <  2:1026915/150‑1 (MQ=255)
                                               tCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGAt                                                                                                     <  2:932015/150‑1 (MQ=255)
                                                cAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAgg                                                                                                                                   >  2:88510/1‑119 (MQ=255)
                                                cAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAgg                                                                                                                                   <  1:88510/119‑1 (MQ=255)
                                                         tAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGt                                                                                           <  1:578568/150‑1 (MQ=255)
                                                         tAAACAGCGCATCCATATGCTGCACACTTAGCAACATAAATCCAAACGGCTGATCCCGCTCATTAACCGCGACCCACAACCCCGCTTCCGCCAGGAAGCAACGAACAAGGTCCTCCAGCTCGGTCCGATACTCTCCTGATACAAAAACGt                                                                                           >  2:666336/1‑150 (MQ=255)
                                                              aGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTg                                                                                      >  2:991197/1‑150 (MQ=255)
                                                               gcgcATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTgg                                                                                     >  1:155129/1‑150 (MQ=255)
                                                                cgcATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCAGGCAGGAAGGAACGCACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGc                                                                                    >  2:1616111/1‑150 (MQ=255)
                                                                 gcATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCa                                                                                   >  2:1185563/1‑150 (MQ=255)
                                                                    tCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCg                                                                                >  1:223209/1‑150 (MQ=255)
                                                                    tCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCg                                                                                <  1:944950/150‑1 (MQ=255)
                                                                      cATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCTAGCTCGGTCCGATACTCTGCTGATTGAAAATCGTGAGTGGCATCGAc                                                                              >  2:414303/1‑150 (MQ=255)
                                                                       atatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACa                                                                             <  1:186777/150‑1 (MQ=255)
                                                                           gctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAAc                                                                         <  1:1267686/150‑1 (MQ=255)
                                                                            ctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACg                                                                        <  1:1555571/150‑1 (MQ=255)
                                                                             tgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACgc                                                                       >  1:1223242/1‑149 (MQ=255)
                                                                             tgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGa                                                                       >  1:1161191/1‑150 (MQ=255)
                                                                             tgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGa                                                                       >  2:1017925/1‑150 (MQ=255)
                                                                              gcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGac                                                                      >  2:138136/1‑150 (MQ=255)
                                                                                 cACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACAcc                                                                   >  2:336018/1‑150 (MQ=255)
                                                                                     tAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAAt                                                               <  1:1565116/150‑1 (MQ=255)
                                                                                        cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGc                                                            <  1:1266275/150‑1 (MQ=255)
                                                                                         aaCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCa                                                           >  2:216985/1‑150 (MQ=255)
                                                                                          aCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCaa                                                          >  1:395689/1‑150 (MQ=255)
                                                                                          aCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATAGTGAGTGGCATCGACAGAACGACACCAAATCGCaa                                                          >  1:1081368/1‑150 (MQ=255)
                                                                                           cATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAAc                                                         >  1:928584/1‑150 (MQ=255)
                                                                                            aTAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACg                                                        >  2:815422/1‑150 (MQ=255)
                                                                                               aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACGACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCACCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGACCGACACCAAATCGCAACGAGt                                                     >  2:817165/1‑150 (MQ=255)
                                                                                               aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGt                                                     >  1:1570204/1‑150 (MQ=255)
                                                                                                   cAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTtcct                                                 >  1:1084511/1‑150 (MQ=255)
                                                                                                   cAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTtcct                                                 >  2:705626/1‑150 (MQ=255)
                                                                                                          cTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTc                                          <  2:395689/150‑1 (MQ=255)
                                                                                                           tGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTccc                                              <  2:486381/145‑1 (MQ=255)
                                                                                                           tGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTccc                                              >  1:486381/1‑145 (MQ=255)
                                                                                                            gATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTTATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCt                                        <  2:1046174/150‑1 (MQ=255)
                                                                                                             aTCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTc                                       <  2:1599619/150‑1 (MQ=255)
                                                                                                              tCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCa                                      >  2:130176/1‑150 (MQ=255)
                                                                                                                     cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTg                               >  2:900550/1‑150 (MQ=255)
                                                                                                                     cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTg                               <  1:166475/150‑1 (MQ=255)
                                                                                                                        tAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAgc                            >  2:1109269/1‑150 (MQ=255)
                                                                                                                        tAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAgc                            >  2:8361/1‑150 (MQ=255)
                                                                                                                        tAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAgc                            >  1:616067/1‑150 (MQ=255)
                                                                                                                        tAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAgc                            >  1:256973/1‑150 (MQ=255)
                                                                                                                         aaCCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAgcg                           >  2:1351023/1‑150 (MQ=255)
                                                                                                                         aaCCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAgcg                           <  1:991197/150‑1 (MQ=255)
                                                                                                                               gACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAAt                     <  2:1444252/150‑1 (MQ=255)
                                                                                                                                aCCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATa                    <  1:666336/150‑1 (MQ=255)
                                                                                                                                  ccACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACt                  <  2:1570204/150‑1 (MQ=255)
                                                                                                                                  ccACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACt                  <  1:1228615/150‑1 (MQ=255)
                                                                                                                                   cacaACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGCACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACta                 >  1:574908/1‑150 (MQ=255)
                                                                                                                                   cacaACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATAGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACta                 >  1:1377371/1‑150 (MQ=255)
                                                                                                                                     caACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCACCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTGCCTCCCCTTCCTCATGACGTCAGCGGCGAATACaaat               >  2:901373/1‑150 (MQ=255)
                                                                                                                                         gttccTTCCTGCCGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACtaata              <  2:1264196/143‑1 (MQ=255)
                                                                                                                                         gttccTTCCTGCCGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACtaata              >  1:1264196/5‑147 (MQ=255)
                                                                                                                                         gttccTTCCTGCCGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAAcc           >  2:67396/5‑150 (MQ=255)
                                                                                                                                             cTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAACCAttt       >  1:1539672/1‑150 (MQ=255)
                                                                                                                                             cTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAACCAttt       >  2:1569298/1‑150 (MQ=255)
                                                                                                                                               tCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAACCATTTtc     >  1:1621358/1‑150 (MQ=255)
                                                                                                                                                ccGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAACCATTTtct    >  1:1548827/1‑150 (MQ=255)
                                                                                                                                                 cGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAACCATTTtctc   >  1:1574588/1‑150 (MQ=255)
                                                                                                                                                 cGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAACCATTTtctc   >  1:1575799/1‑150 (MQ=255)
                                                                                                                                                 cGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAACCATTTtctc   >  1:228754/1‑150 (MQ=255)
                                                                                                                                                  ggCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAACCATTTtctct  <  1:900550/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTCCTCCCCTTCCTCATGCCGTGAGCGGCGAATACTAATAACCATTTTCTCT  >  NC_000913/4213835‑4214130

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