Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 582,152 C→T pseudogene (653/653 nt) ybcY ← pseudogene, DLP12 prophage, methyltransferase homology,Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913582,1520CT98.9% 297.7 / NA 93pseudogene (653/653 nt)ybcYpseudogene, DLP12 prophage, methyltransferase homology,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (47/45);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTTCTGATT  >  NC_000913/582034‑582286
                                                                                                                      |                                                                                                                                      
gAACCGTGTTGATACATCAACTGCACGTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   <  2:1055694/139‑1 (MQ=255)
gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   >  1:492668/1‑139 (MQ=255)
gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   <  1:357870/139‑1 (MQ=255)
     gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              <  1:1370425/139‑1 (MQ=255)
     gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              <  2:492668/139‑1 (MQ=255)
        ttGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATAATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATtac                                                                                                           >  2:210193/1‑139 (MQ=255)
            tACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                       <  2:1526278/139‑1 (MQ=255)
            tACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                       <  1:19864/139‑1 (MQ=255)
             aCATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAg                                                                                                      <  2:1444612/139‑1 (MQ=255)
                tCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTAc                                                                                                   >  1:1190481/1‑139 (MQ=255)
                        aCCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                           <  1:210003/139‑1 (MQ=255)
                        aCCTGAGATTGAGTGGCCTACGAACCCTGACAGGGAGTAATAATTGGGACTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTAAGGTTATTTTTTTCCTGTAGTGGAAAAAAATCCTATAGTACTTTGaata                                                                                           >  2:1561179/1‑137 (MQ=255)
                         ccTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                          <  2:525685/139‑1 (MQ=255)
                                  gAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                                 <  1:145293/139‑1 (MQ=255)
                                   aGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATtact                                                                                                          >  1:1495790/1‑113 (MQ=255)
                                   aGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATtact                                                                                                          <  2:1495790/113‑1 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTg                                                                                                                           >  2:1448327/1‑92 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTg                                                                                                                           <  1:1448327/92‑1 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCgg                                                                                                                     >  2:775761/1‑98 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCgg                                                                                                                     <  1:775761/98‑1 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc                                                                            >  2:675685/1‑139 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc                                                                            >  1:278168/1‑139 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   <  1:401001/99‑1 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   >  2:401001/1‑99 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              >  1:751621/1‑104 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              <  2:751621/104‑1 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              <  1:1542832/104‑1 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtct                                                                           >  2:1420897/1‑139 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATAATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATACTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              >  2:1542832/1‑104 (MQ=255)
                                         cTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtctc                                                                          <  1:210193/139‑1 (MQ=255)
                                          tACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCg                                                                         >  1:497506/1‑139 (MQ=255)
                                                 ccTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                                  <  1:1533607/139‑1 (MQ=255)
                                                   tGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                          <  1:1436785/113‑1 (MQ=255)
                                                   tGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                          >  2:1436785/1‑113 (MQ=255)
                                                               aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGa                                                                        >  2:385476/1‑119 (MQ=255)
                                                               aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGa                                                                        <  1:385476/119‑1 (MQ=255)
                                                               aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATAtgt                                                    >  2:1139343/1‑139 (MQ=255)
                                                                     gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAg                                              >  1:1216352/1‑139 (MQ=255)
                                                                       gATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGagag                                                           <  1:1266943/124‑1 (MQ=255)
                                                                       gATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGagag                                                           >  2:1266943/1‑124 (MQ=255)
                                                                        aTTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                                 >  2:797734/1‑69 (MQ=255)
                                                                        aTTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                                 <  1:797734/69‑1 (MQ=255)
                                                                          tATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCcttc                                         >  1:1604654/1‑139 (MQ=255)
                                                                            tGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                       <  1:721346/75‑1 (MQ=255)
                                                                            tGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                       >  2:721346/1‑75 (MQ=255)
                                                                                 cagcaCGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   >  2:167037/1‑58 (MQ=255)
                                                                                 cagcaCGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   <  1:167037/58‑1 (MQ=255)
                                                                                 cagcaCGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              <  1:924125/111‑1 (MQ=255)
                                                                                 cagcaCGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              >  2:924125/1‑111 (MQ=255)
                                                                                  agcaCGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                                 <  2:966522/91‑1 (MQ=255)
                                                                                  agcaCGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                                 >  1:966522/1‑91 (MQ=255)
                                                                                     aCGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCtt                              <  1:776677/139‑1 (MQ=255)
                                                                                      cGTCTTTAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTAGAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCctt                                          <  2:1083477/126‑1 (MQ=255)
                                                                                      cGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                            <  2:573514/60‑1 (MQ=255)
                                                                                      cGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                            >  1:573514/1‑60 (MQ=255)
                                                                                      cGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCctt                                          >  1:1083477/1‑126 (MQ=255)
                                                                                      cGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCttt                             <  2:497506/139‑1 (MQ=255)
                                                                                        tCTTAAGCAAGAACATGCTCCGGTTGATTCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTAGTTTTAACATTCTCGAAATGCTCTGACAGTATATGTGTTAAGCCTTCTTCGGACTCTTTtg                           <  2:1069368/139‑1 (MQ=255)
                                                                                         cTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTtgt                          >  2:1253091/1‑139 (MQ=255)
                                                                                              gCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTg                     >  1:407/1‑139 (MQ=255)
                                                                                                 aGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGaaa                  <  1:1253091/139‑1 (MQ=255)
                                                                                                  gAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTg                                                               >  1:648379/1‑93 (MQ=255)
                                                                                                  gAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTg                                                               <  2:648379/93‑1 (MQ=255)
                                                                                                  gAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGaaaa                 >  1:1263745/1‑139 (MQ=255)
                                                                                                   aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                                      >  1:1239086/1‑85 (MQ=255)
                                                                                                   aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                                      <  2:1239086/85‑1 (MQ=255)
                                                                                                   aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              <  1:496262/93‑1 (MQ=255)
                                                                                                   aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              <  2:715670/93‑1 (MQ=255)
                                                                                                   aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              >  1:715670/1‑93 (MQ=255)
                                                                                                   aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              >  2:496262/1‑93 (MQ=255)
                                                                                                   aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAg                >  1:630037/1‑139 (MQ=255)
                                                                                                   aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAg                >  2:487229/1‑139 (MQ=255)
                                                                                                     cATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGAt              <  2:407/139‑1 (MQ=255)
                                                                                                         ctgcGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCctt                                          >  1:1435148/1‑107 (MQ=255)
                                                                                                         ctgcGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCctt                                          <  2:1435148/107‑1 (MQ=255)
                                                                                                         ctgcGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCt          >  2:1108698/1‑139 (MQ=255)
                                                                                                         ctgcGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCt          >  1:153936/1‑139 (MQ=255)
                                                                                                          tgcGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              <  1:152699/38‑1 (MQ=25)
                                                                                                          tgcGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              >  2:152699/1‑38 (MQ=25)
                                                                                                            cGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATtac                                                                                                           >  1:234760/1‑39 (MQ=25)
                                                                                                            cGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATtac                                                                                                           <  2:234760/39‑1 (MQ=25)
                                                                                                              ggttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                           >  1:1556462/5‑54 (MQ=38)
                                                                                                              ggttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                           <  2:1556462/50‑1 (MQ=39)
                                                                                                              ggttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                                  <  2:462648/75‑1 (MQ=255)
                                                                                                              ggttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                                  >  1:462648/5‑79 (MQ=255)
                                                                                                              ggttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                       <  1:265671/86‑1 (MQ=255)
                                                                                                              ggttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                       >  2:265671/5‑90 (MQ=255)
                                                                                                              ggttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTTCt      <  2:630037/135‑1 (MQ=255)
                                                                                                              ggttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTTCt      <  2:1033623/135‑1 (MQ=255)
                                                                                                               gttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTTCTg     >  2:250974/4‑139 (MQ=255)
                                                                                                                 tgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              <  1:1545088/79‑1 (MQ=255)
                                                                                                                 tgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              >  2:1545088/2‑80 (MQ=255)
                                                                                                                  gATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTTCTGAtt  >  1:635053/1‑139 (MQ=255)
                                                                                                                      |                                                                                                                                      
GAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTTCTGATT  >  NC_000913/582034‑582286

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: