Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 582,148 Δ1 bp intergenic (+51/+4) tfaD → / ← ybcY pseudogene, DLP12 prophage, tail fiber assembly protein family,Phage or Prophage Related/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,1470G.97.5% 308.7 / 1.1 79intergenic (+50/+5)tfaD/ybcYpseudogene, DLP12 prophage, tail fiber assembly protein family,Phage or Prophage Related/pseudogene, DLP12 prophage, methyltransferase homology,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (38/39);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

CTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTT  >  NC_000913/582010‑582279
                                                                                                                                         |                                                                                                                                    
cTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTga                                                                                                                                     >  2:638148/1‑138 (MQ=255)
 tGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTgat                                                                                                                                    <  1:1660278/139‑3 (MQ=255)
  ggAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTgatg                                                                                                                                   >  1:660862/1‑136 (MQ=255)
   gAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTgatgt                                                                                                                                  <  1:1758866/139‑5 (MQ=255)
   gAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTgatgt                                                                                                                                  >  1:640314/1‑135 (MQ=255)
                   ttGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAg                                                                                                                 >  2:297365/1‑139 (MQ=255)
                      cTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCGGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCgg                                                                                                              >  1:702875/1‑139 (MQ=255)
                        gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                            <  2:640314/139‑1 (MQ=255)
                        gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                            <  1:287905/139‑1 (MQ=255)
                              tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAtt                                                                                                      >  2:1866214/1‑139 (MQ=255)
                              tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAtt                                                                                                      >  2:1360663/1‑139 (MQ=255)
                                  gATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATtacta                                                                                                  <  1:297365/139‑1 (MQ=255)
                                    tACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                <  1:1805554/139‑1 (MQ=255)
                                        tCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAATCGGAAAACATTACTACAGTAc                                                                                            >  2:1106835/1‑139 (MQ=255)
                                                aCCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                    <  1:1360663/139‑1 (MQ=255)
                                                aCCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                    <  1:1877493/139‑1 (MQ=255)
                                                 ccTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAc                                                                                     >  1:142705/1‑137 (MQ=255)
                                                 ccTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAc                                                                                     <  2:142705/137‑1 (MQ=255)
                                                            gTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                   >  1:1873051/1‑128 (MQ=255)
                                                            gTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                   <  2:1873051/128‑1 (MQ=255)
                                                              ggCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGG‑TGTTATTTTTTTCCTGAAGCGGa                                                                                                             <  1:16642/100‑1 (MQ=255)
                                                              ggCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGG‑TGTTATTTTTTTCCTGAAGCGGa                                                                                                             >  2:16642/1‑100 (MQ=255)
                                                              ggCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACAtt                                                                      >  1:1013718/1‑139 (MQ=255)
                                                               gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                            <  1:206519/100‑1 (MQ=255)
                                                               gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                            >  2:206519/1‑100 (MQ=255)
                                                               gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc                                                                     <  2:1013718/139‑1 (MQ=255)
                                                               gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAACACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc                                                                     >  2:387597/1‑139 (MQ=255)
                                                                       aCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATg                                                             >  1:354681/1‑139 (MQ=255)
                                                                         ccTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                           <  2:1622320/139‑1 (MQ=255)
                                                                           tGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                         <  1:387597/139‑1 (MQ=255)
                                                                                  gAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                            >  2:612424/1‑81 (MQ=255)
                                                                                  gAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                            <  1:612424/81‑1 (MQ=255)
                                                                                    gTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                <  1:1866214/139‑1 (MQ=255)
                                                                                      aaTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                            >  2:685740/1‑77 (MQ=255)
                                                                                      aaTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                            <  1:685740/77‑1 (MQ=255)
                                                                                       aTCATTGGGATTATGCTGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCtt                                                                                  >  2:778406/1‑102 (MQ=255)
                                                                                       aTCATTGGGATTATGCTGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCtt                                                                                  <  1:778406/102‑1 (MQ=255)
                                                                                             gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                          >  2:286862/1‑72 (MQ=255)
                                                                                             gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                          <  1:286862/72‑1 (MQ=255)
                                                                                             gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                               >  2:1461790/1‑115 (MQ=255)
                                                                                             gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                               <  1:1461790/115‑1 (MQ=255)
                                                                                             gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAg                                       >  2:914197/1‑139 (MQ=255)
                                                                                               gcttATGCCTCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCCGAGAGTATATGTGTTAAGcc                                     <  1:1106835/137‑1 (MQ=255)
                                                                                               gATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                       <  1:1398424/73‑1 (MQ=255)
                                                                                               gATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                       >  2:1398424/1‑73 (MQ=255)
                                                                                                aTTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCct                                    >  2:141346/1‑139 (MQ=255)
                                                                                                  tATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCcttc                                  >  1:1837747/1‑139 (MQ=255)
                                                                                                  tATGCCGCAGCACCTCTTAAGCAAGCACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGCACCTTGGTTTTAACATTCTCGACATGCTCTGAGAGTATATGTGTTAAGCctta                                  >  1:292819/1‑138 (MQ=255)
                                                                                                   aTGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTGTTTTTTTCCTGAAGCGGaaaa                                                                                                          >  1:1327720/1‑66 (MQ=255)
                                                                                                   aTGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                          <  2:1327720/66‑1 (MQ=255)
                                                                                                   aTGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                         >  1:1115152/1‑115 (MQ=255)
                                                                                                   aTGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                         <  2:1115152/115‑1 (MQ=255)
                                                                                                    tGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCcttctt                                >  2:1419971/1‑139 (MQ=255)
                                                                                                       cGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCgg                             <  2:354681/139‑1 (MQ=255)
                                                                                                         cagcaCGTCTTAAGCAAGAACATGCTGCGGTTgatg                                                                                                                                   <  2:1707203/36‑4 (MQ=255)
                                                                                                         cagcaCGTCTTAAGCAAGAACATGCTGCGGTTgatg                                                                                                                                   >  1:1707203/1‑33 (MQ=255)
                                                                                                         cagcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAg                                                                                               >  1:1484044/1‑71 (MQ=255)
                                                                                                         cagcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAg                                                                                               <  2:1484044/71‑1 (MQ=255)
                                                                                                           gcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                          >  2:111359/1‑90 (MQ=255)
                                                                                                           gcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                          <  1:111359/90‑1 (MQ=255)
                                                                                                              cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCttt                      <  1:1419971/139‑1 (MQ=255)
                                                                                                                         aGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                       <  1:711678/95‑1 (MQ=255)
                                                                                                                         aGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                       >  2:711678/1‑95 (MQ=255)
                                                                                                                          gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATtacta                                                                                                  >  2:1671134/1‑51 (MQ=38)
                                                                                                                          gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATtacta                                                                                                  <  1:1671134/51‑1 (MQ=38)
                                                                                                                          gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                         <  2:288177/92‑1 (MQ=255)
                                                                                                                          gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                         >  1:288177/1‑92 (MQ=255)
                                                                                                                           aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                >  2:979700/1‑52 (MQ=38)
                                                                                                                           aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                <  1:979700/52‑1 (MQ=38)
                                                                                                                           aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                          >  2:1706990/1‑74 (MQ=255)
                                                                                                                           aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                          <  1:1706990/74‑1 (MQ=255)
                                                                                                                           aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                       >  1:645318/1‑93 (MQ=255)
                                                                                                                           aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                       <  2:645318/93‑1 (MQ=255)
                                                                                                                           aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGaaatg         >  1:1103982/1‑137 (MQ=255)
                                                                                                                           aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAg         <  2:1837747/139‑1 (MQ=255)
                                                                                                                             cATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGAt       <  1:911604/139‑1 (MQ=255)
                                                                                                                                gctgcGGTT‑GTTGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                    <  2:1444174/59‑1 (MQ=255)
                                                                                                                                gctgcGGTT‑GTTGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                    >  1:1444174/1‑59 (MQ=255)
                                                                                                                                gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAg         >  2:250689/1‑134 (MQ=255)
                                                                                                                                gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAg         <  1:250689/134‑1 (MQ=255)
                                                                                                                                 ctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATgtgt                                           <  1:1635617/99‑1 (MQ=255)
                                                                                                                                 ctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATgtgt                                           >  2:1635617/1‑99 (MQ=255)
                                                                                                                                 ctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCt   <  2:702875/139‑1 (MQ=255)
                                                                                                                                  tgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCtt  >  1:145695/1‑139 (MQ=255)
                                                                                                                                    cGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                       >  1:404725/1‑84 (MQ=255)
                                                                                                                                    cGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                       <  2:404725/84‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
CTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTT  >  NC_000913/582010‑582279

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