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
                                                                                                                      |                                                                                                                                      
gAACCGTGTTGATACATCAACTGCACGTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   <  2:1055694/139‑1 (MQ=255)
gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   >  1:492668/1‑139 (MQ=255)
gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   <  1:357870/139‑1 (MQ=255)
     gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              <  1:1370425/139‑1 (MQ=255)
     gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              <  2:492668/139‑1 (MQ=255)
        ttGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATAATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATtac                                                                                                           >  2:210193/1‑139 (MQ=255)
            tACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                       <  2:1526278/139‑1 (MQ=255)
            tACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                       <  1:19864/139‑1 (MQ=255)
             aCATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAg                                                                                                      <  2:1444612/139‑1 (MQ=255)
                tCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTAc                                                                                                   >  1:1190481/1‑139 (MQ=255)
                        aCCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                           <  1:210003/139‑1 (MQ=255)
                        aCCTGAGATTGAGTGGCCTACGAACCCTGACAGGGAGTAATAATTGGGACTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑AAGGTTATTTTTTTCCTGTAGTGGAAAAAAATCCTATAGTACTTTGaata                                                                                           >  2:1561179/1‑137 (MQ=255)
                         ccTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                          <  2:525685/139‑1 (MQ=255)
                                  gAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                                 <  1:145293/139‑1 (MQ=255)
                                   aGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATtact                                                                                                          >  1:1495790/1‑113 (MQ=255)
                                   aGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATtact                                                                                                          <  2:1495790/113‑1 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTg                                                                                                                           >  2:1448327/1‑92 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTg                                                                                                                           <  1:1448327/92‑1 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCgg                                                                                                                     >  2:775761/1‑98 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCgg                                                                                                                     <  1:775761/98‑1 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc                                                                            >  2:675685/1‑139 (MQ=255)
                                       gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc                                                                            >  1:278168/1‑139 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   <  1:401001/99‑1 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   >  2:401001/1‑99 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              >  1:751621/1‑104 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              <  2:751621/104‑1 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              <  1:1542832/104‑1 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtct                                                                           >  2:1420897/1‑139 (MQ=255)
                                        ccTACGAACCCTGTCAGGGAGTAATAATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATACTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              >  2:1542832/1‑104 (MQ=255)
                                         cTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtctc                                                                          <  1:210193/139‑1 (MQ=255)
                                          tACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCg                                                                         >  1:497506/1‑139 (MQ=255)
                                                 ccTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                                  <  1:1533607/139‑1 (MQ=255)
                                                   tGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                          <  1:1436785/113‑1 (MQ=255)
                                                   tGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                          >  2:1436785/1‑113 (MQ=255)
                                                               aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGa                                                                        >  2:385476/1‑119 (MQ=255)
                                                               aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGa                                                                        <  1:385476/119‑1 (MQ=255)
                                                               aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATAtgt                                                    >  2:1139343/1‑139 (MQ=255)
                                                                     gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAg                                              >  1:1216352/1‑139 (MQ=255)
                                                                       gATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGagag                                                           <  1:1266943/124‑1 (MQ=255)
                                                                       gATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGagag                                                           >  2:1266943/1‑124 (MQ=255)
                                                                        aTTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                                 >  2:797734/1‑69 (MQ=255)
                                                                        aTTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                                 <  1:797734/69‑1 (MQ=255)
                                                                          tATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCcttc                                         >  1:1604654/1‑139 (MQ=255)
                                                                            tGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                       <  1:721346/75‑1 (MQ=255)
                                                                            tGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                       >  2:721346/1‑75 (MQ=255)
                                                                                 cagcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   >  2:167037/1‑58 (MQ=255)
                                                                                 cagcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                                   <  1:167037/58‑1 (MQ=255)
                                                                                 cagcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              <  1:924125/111‑1 (MQ=255)
                                                                                 cagcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              >  2:924125/1‑111 (MQ=255)
                                                                                  agcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                                 <  2:966522/91‑1 (MQ=255)
                                                                                  agcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                                 >  1:966522/1‑91 (MQ=255)
                                                                                     aCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCtt                              <  1:776677/139‑1 (MQ=255)
                                                                                      cGTCTTTAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTAGAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCctt                                          <  2:1083477/126‑1 (MQ=255)
                                                                                      cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                            <  2:573514/60‑1 (MQ=255)
                                                                                      cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                            >  1:573514/1‑60 (MQ=255)
                                                                                      cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCctt                                          >  1:1083477/1‑126 (MQ=255)
                                                                                      cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCttt                             <  2:497506/139‑1 (MQ=255)
                                                                                        tCTTAAGCAAGAACATGCTCCGGTT‑GATTCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTAGTTTTAACATTCTCGAAATGCTCTGACAGTATATGTGTTAAGCCTTCTTCGGACTCTTTtg                           <  2:1069368/139‑1 (MQ=255)
                                                                                         cTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTtgt                          >  2:1253091/1‑139 (MQ=255)
                                                                                              gCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTg                     >  1:407/1‑139 (MQ=255)
                                                                                                 aGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGaaa                  <  1:1253091/139‑1 (MQ=255)
                                                                                                  gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTg                                                               >  1:648379/1‑93 (MQ=255)
                                                                                                  gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTg                                                               <  2:648379/93‑1 (MQ=255)
                                                                                                  gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGaaaa                 >  1:1263745/1‑139 (MQ=255)
                                                                                                   aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                                      >  1:1239086/1‑85 (MQ=255)
                                                                                                   aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                                      <  2:1239086/85‑1 (MQ=255)
                                                                                                   aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              <  1:496262/93‑1 (MQ=255)
                                                                                                   aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              <  2:715670/93‑1 (MQ=255)
                                                                                                   aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              >  1:715670/1‑93 (MQ=255)
                                                                                                   aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                              >  2:496262/1‑93 (MQ=255)
                                                                                                   aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAg                >  1:630037/1‑139 (MQ=255)
                                                                                                   aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAg                >  2:487229/1‑139 (MQ=255)
                                                                                                     cATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGAt              <  2:407/139‑1 (MQ=255)
                                                                                                         ctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCctt                                          >  1:1435148/1‑107 (MQ=255)
                                                                                                         ctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCctt                                          <  2:1435148/107‑1 (MQ=255)
                                                                                                         ctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCt          >  2:1108698/1‑139 (MQ=255)
                                                                                                         ctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCt          >  1:153936/1‑139 (MQ=255)
                                                                                                          tgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              <  1:152699/38‑1 (MQ=25)
                                                                                                          tgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                              >  2:152699/1‑38 (MQ=25)
                                                                                                            cGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATtac                                                                                                           >  1:234760/1‑39 (MQ=25)
                                                                                                            cGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATtac                                                                                                           <  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: ‑