Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,074,810 Δ1 bp 8.7% coding (32/1533 nt) yeeR → CP4‑44 prophage; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,074,8100T. 8.7% 330.6 / 23.3 92coding (32/1533 nt)yeeRCP4‑44 prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base T (42/42);  new base . (4/4);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGA  >  NC_000913/2074684‑2074933
                                                                                                                              |                                                                                                                            
cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                        >  3:67046/1‑148 (MQ=255)
 ggTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttctt                                                                                                       >  3:105838/1‑148 (MQ=255)
      tcatGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCt                                                                                                  >  1:869447/1‑148 (MQ=255)
        atGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCGATCAGCAGGAt                                                                                                                  >  4:869599/1‑130 (MQ=255)
        atGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCGATCAGCAGGAt                                                                                                                  <  3:869599/130‑1 (MQ=255)
         tGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctgctg                                                                                                                             >  2:873528/1‑119 (MQ=255)
         tGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctgctg                                                                                                                             <  1:873528/119‑1 (MQ=255)
               ggACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCag                                                                                         <  4:67046/148‑1 (MQ=255)
                gACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGcc                                                                                                   <  1:1049759/137‑1 (MQ=255)
                gACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGcc                                                                                                   >  2:1049759/1‑137 (MQ=255)
                  cccGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGc                                                                                      <  2:869447/148‑1 (MQ=255)
                    cGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                        >  4:993906/1‑128 (MQ=255)
                    cGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                        >  3:560143/1‑128 (MQ=255)
                    cGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                        <  3:993906/128‑1 (MQ=255)
                    cGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                        <  4:560143/128‑1 (MQ=255)
                    cGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCAt                                                                                    >  2:825041/1‑148 (MQ=255)
                      gACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTa                                                                                  >  1:1114458/1‑148 (MQ=255)
                        ccGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGAGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAAt                                                                                <  1:825041/148‑1 (MQ=255)
                             gcAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCtttt                                                                                               <  2:97806/128‑1 (MQ=255)
                             gcAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCtttt                                                                                               >  1:97806/1‑128 (MQ=255)
                                   gATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCa                                                                                     <  2:219981/132‑1 (MQ=255)
                                   gATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCa                                                                                     >  1:219981/1‑132 (MQ=255)
                                   gATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCtt                                                                     <  2:1114458/148‑1 (MQ=255)
                                        tCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGAt                                                                                                                           >  2:17853/1‑90 (MQ=255)
                                        tCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGAt                                                                                                                           <  1:17853/90‑1 (MQ=255)
                                           aCGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCGATCAGCAGGa                                                                                                                   <  4:1021465/94‑1 (MQ=255)
                                           aCGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCGATCAGCAGGa                                                                                                                   >  3:1021465/1‑94 (MQ=255)
                                            cGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTg                                                                  <  1:214760/142‑1 (MQ=255)
                                            cGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTg                                                                  >  2:214760/1‑142 (MQ=255)
                                              ccACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCa                                                          >  1:1127165/1‑148 (MQ=255)
                                                acaTTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAgg                                                        >  3:479445/1‑148 (MQ=255)
                                                       cTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTc                                                                                           <  3:1133869/106‑1 (MQ=255)
                                                       cTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTc                                                                                           >  4:1133869/1‑106 (MQ=255)
                                                           cATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                             >  1:67136/1‑148 (MQ=255)
                                                           cATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTGTGGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATAGTGc                                             <  4:105838/148‑1 (MQ=255)
                                                             ttAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCttt                                                                                                <  1:953292/95‑1 (MQ=255)
                                                             ttAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCttt                                                                                                >  2:953292/1‑95 (MQ=255)
                                                               aaCAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtgtg                                         <  2:1127165/148‑1 (MQ=255)
                                                                 cAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTg                                                                                                               <  4:908263/76‑1 (MQ=255)
                                                                 cAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTg                                                                                                               >  3:908263/1‑76 (MQ=255)
                                                                 cAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTc                                       <  4:479445/148‑1 (MQ=255)
                                                                          gATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                             <  1:1022050/133‑1 (MQ=255)
                                                                          gATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                             >  2:1022050/1‑133 (MQ=255)
                                                                           atatTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                        <  2:397157/73‑1 (MQ=255)
                                                                           atatTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                        >  1:397157/1‑73 (MQ=255)
                                                                              ttCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTgcg                                                                >  3:347997/1‑110 (MQ=255)
                                                                              ttCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTgcg                                                                <  4:347997/110‑1 (MQ=255)
                                                                               tcatcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGAt                                                                                                                           <  1:1124301/51‑1 (MQ=255)
                                                                               tcatcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGAt                                                                                                                           >  2:1124301/1‑51 (MQ=255)
                                                                               tcatcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTc                                                                                           >  3:65813/1‑82 (MQ=255)
                                                                               tcatcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTc                                                                                           <  4:65813/82‑1 (MQ=255)
                                                                                catcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCtt                                                                                                 >  1:413524/1‑75 (MQ=255)
                                                                                catcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCtt                                                                                                 <  2:413524/75‑1 (MQ=255)
                                                                                  tcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCa                                                                                     >  3:942125/1‑85 (MQ=255)
                                                                                  tcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCa                                                                                     <  4:942125/85‑1 (MQ=255)
                                                                                   caCGGAGTGACTATGTTACAGATAGTCCGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTAc                     >  3:73113/1‑148 (MQ=255)
                                                                                     cGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGc                                                                            >  3:934590/1‑91 (MQ=255)
                                                                                     cGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGc                                                                            <  4:934590/91‑1 (MQ=255)
                                                                                     cGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAgg                                                        <  4:222319/111‑1 (MQ=255)
                                                                                     cGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAgg                                                        >  3:222319/1‑111 (MQ=255)
                                                                                     cGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActg                   >  3:256568/1‑148 (MQ=255)
                                                                                        aGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActggct                <  4:73113/148‑1 (MQ=255)
                                                                                           gACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtg                                           >  2:83818/1‑118 (MQ=255)
                                                                                           gACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtg                                           <  1:83818/118‑1 (MQ=255)
                                                                                             cTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTg                                                                  >  1:206509/1‑93 (MQ=255)
                                                                                             cTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTg                                                                  <  2:206509/93‑1 (MQ=255)
                                                                                             cTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCta           >  3:452842/1‑148 (MQ=255)
                                                                                               aTGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCa                                                                                          <  1:100725/67‑1 (MQ=255)
                                                                                               aTGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCa                                                                                          >  2:100725/1‑67 (MQ=255)
                                                                                               aTGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGc                                                                                      >  2:156525/1‑71 (MQ=255)
                                                                                               aTGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGc                                                                                      <  1:156525/71‑1 (MQ=255)
                                                                                               aTGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtg                                           <  3:842951/114‑1 (MQ=255)
                                                                                               aTGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtg                                           >  4:842951/1‑114 (MQ=255)
                                                                                                tGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGTATTAAt                                                                                >  1:565474/1‑76 (MQ=255)
                                                                                                tGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGTATTAAt                                                                                <  2:565474/76‑1 (MQ=255)
                                                                                                 gTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATc       <  4:256568/148‑1 (MQ=255)
                                                                                                   tACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActggctgg              <  4:764933/139‑1 (MQ=255)
                                                                                                   tACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActggctgg              >  3:764933/1‑139 (MQ=255)
                                                                                                   tACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCAc     <  4:452842/148‑1 (MQ=255)
                                                                                                    aCAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCAcc    <  2:67136/148‑1 (MQ=255)
                                                                                                     cAGATAGTCGGCGCGCTGCTCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCag                                                                                         <  3:607907/62‑1 (MQ=255)
                                                                                                     cAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCag                                                                                         >  4:607907/1‑62 (MQ=255)
                                                                                                        aTAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcccgg                           <  3:798837/121‑1 (MQ=255)
                                                                                                        aTAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcccgg                           >  4:798837/1‑121 (MQ=255)
                                                                                                             cGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCt                                                                      >  2:926285/1‑73 (MQ=255)
                                                                                                             cGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCt                                                                      <  1:926285/73‑1 (MQ=255)
                                                                                                             cGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGc                                                      <  1:529745/89‑1 (MQ=255)
                                                                                                             cGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGc                                                      >  2:529745/1‑89 (MQ=255)
                                                                                                                      aaTCCTGCGATCAGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCTCGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCAGCCTTACTGGCTGGCTATATCACCGa  >  1:916031/2‑132 (MQ=255)
                                                                                                                      aaTCCTGCGATCAGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCCGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGa  <  2:916031/131‑1 (MQ=255)
                                                                                                                      aaTCCTGCGATCAGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActggctgg              <  2:237659/119‑1 (MQ=255)
                                                                                                                      aaTCCTGCGATCAGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActggctgg              >  1:237659/2‑120 (MQ=255)
                                                                                                                       aTCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcccggcc                         <  2:650863/108‑1 (MQ=255)
                                                                                                                       aTCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcccggcc                         >  1:650863/1‑108 (MQ=255)
                                                                                                                          ctgctgATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATAtttc                                             <  1:978320/85‑3 (MQ=255)
                                                                                                                          ctgctgATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                             >  2:978320/1‑85 (MQ=255)
                                                                                                                           tgctgATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACtggctggc             >  2:207331/1‑116 (MQ=255)
                                                                                                                           tgctgATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACtggctggc             <  1:207331/116‑1 (MQ=255)
                                                                                                                             ctgATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtg                                           <  1:359164/84‑1 (MQ=255)
                                                                                                                             ctgATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtg                                           >  2:359164/1‑84 (MQ=255)
                                                                                                                              |                                                                                                                            
CGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGA  >  NC_000913/2074684‑2074933

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