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_000913	582,147	0	G	.	100.0%	293.8 / NA	72	intergenic (+50/+5)	tfaD/ybcY	pseudogene, 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 (0/0);  new base . (33/39);  total (33/39)

CTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTT  >  NC_000913/582010‑582280                                                                                                                                           |                                                                                                                                      cTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTga                                                                                                                                      >  2:1709785/1‑138 (MQ=255)    gAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTgatgt                                                                                                                                   >  2:1066641/1‑135 (MQ=255)              cGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTc                                                                                                                        >  1:573252/1‑139 (MQ=255)                 gtgtTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTg                                                                                                                     <  2:573252/139‑1 (MQ=255)                   gtTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGaa                                                                                                                   <  2:486316/139‑1 (MQ=255)                   gtTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGaa                                                                                                                   >  1:593417/1‑139 (MQ=255)                    ttGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAg                                                                                                                  >  2:831882/1‑139 (MQ=255)                         gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             <  2:1484474/139‑1 (MQ=255)                         gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             <  1:831882/139‑1 (MQ=255)                         gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             <  2:1006108/139‑1 (MQ=255)                              gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        <  1:453008/139‑1 (MQ=255)                              gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        <  2:1007125/139‑1 (MQ=255)                               tgtTGATACATCAACTGCACCTGATATTGCGTGGCCTACGAACCCTGTCAGGGAGTAATCAGT‑AGGAGTATGCCAGCGCACGTCTTAAGCAAGAACATGCTGCGGGT‑GATGTTATTTTTTTCCTGCACCGGAAAACAtt                                                                                                       >  2:444387/1‑139 (MQ=255)                               tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAtt                                                                                                       >  2:1100143/1‑139 (MQ=255)                               tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGCAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAtt                                                                                                       <  1:1145256/139‑1 (MQ=255)                                gtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                      <  2:1627600/139‑1 (MQ=255)                                         tCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                 >  2:876144/1‑135 (MQ=255)                                         tCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                 <  1:876144/135‑1 (MQ=255)                                                   cTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCtt                                                                                   <  1:1709785/139‑1 (MQ=255)                                                                gCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTAAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACa                                                                         <  2:67563/136‑1 (MQ=255)                                                                gCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTAAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACa                                                                         >  1:67563/1‑136 (MQ=255)                                                                 ccTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             >  2:1197664/1‑99 (MQ=255)                                                                 ccTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             <  1:1197664/99‑1 (MQ=255)                                                                 ccTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtct                                                                     >  1:323190/1‑139 (MQ=255)                                                                      gAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                                <  1:444387/139‑1 (MQ=255)                                                                           cTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                           >  2:1486282/1‑123 (MQ=255)                                                                           cTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                           <  1:1486282/123‑1 (MQ=255)                                                                           cTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctc                                                           >  2:1505105/1‑139 (MQ=255)                                                                           cTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                            >  2:1582529/1‑138 (MQ=255)                                                                           cTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                            <  1:1582529/138‑1 (MQ=255)                                                                            tGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          <  2:1048535/139‑1 (MQ=255)                                                                            tGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          >  1:895058/1‑139 (MQ=255)                                                                            tGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          <  1:1066641/139‑1 (MQ=255)                                                                                aGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTgaga                                                      <  2:323190/139‑1 (MQ=255)                                                                                 gggAGTAATCATTGGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTgaga                                                      <  1:1630801/139‑1 (MQ=255)                                                                                 gggAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGagag                                                     <  2:895058/139‑1 (MQ=255)                                                                                  ggAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAGGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGt                                                    >  2:141497/1‑139 (MQ=255)                                                                                  ggAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTAttt                                                                                                                             >  1:1229219/1‑66 (MQ=255)                                                                                  ggAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTAttt                                                                                                                             <  2:1229219/66‑1 (MQ=255)                                                                                    aGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtat                                                  <  1:1505105/139‑1 (MQ=255)                                                                                      tAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                 <  2:197750/90‑1 (MQ=255)                                                                                      tAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                 >  1:197750/1‑90 (MQ=255)                                                                                      tAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                 >  2:736518/1‑138 (MQ=255)                                                                                      tAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                 <  1:736518/138‑1 (MQ=255)                                                                                      tAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTatat                                                >  1:104750/1‑139 (MQ=255)                                                                                        aTCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                                >  1:1704290/1‑121 (MQ=255)                                                                                        aTCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                                <  2:1704290/121‑1 (MQ=255)                                                                                            tt‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTgatgt                                                                                                                                   >  1:1048948/1‑47 (MQ=255)                                                                                            tt‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTgatgt                                                                                                                                   <  2:1048948/51‑5 (MQ=255)                                                                                             t‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTaa                                         >  1:805127/1‑139 (MQ=255)                                                                                               gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAg                                        >  2:1166496/1‑139 (MQ=255)                                                                                               gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAg                                        >  2:345164/1‑139 (MQ=255)                                                                                                ggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGc                                       <  2:805127/139‑1 (MQ=255)                                                                                                 gcttATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGcc                                      >  2:831435/3‑139 (MQ=255)                                                                                                   ttATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCctt                                    <  1:141497/139‑1 (MQ=255)                                                                                                      tGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          <  2:1209156/114‑1 (MQ=255)                                                                                                      tGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          >  1:1209156/1‑114 (MQ=255)                                                                                                          gcagcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGa                             <  1:831435/139‑1 (MQ=255)                                                                                                                cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             >  1:1574772/1‑53 (MQ=38)                                                                                                                cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             <  2:1574772/53‑1 (MQ=38)                                                                                                                cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        >  2:1712286/1‑58 (MQ=255)                                                                                                                cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        <  1:1712286/58‑1 (MQ=255)                                                                                                                 gTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCtttt                      >  2:54968/1‑139 (MQ=255)                                                                                                                   cTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                           <  1:1355014/52‑1 (MQ=38)                                                                                                                   cTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                           >  2:1355014/1‑52 (MQ=38)                                                                                                                            gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGaaaa           >  2:1400375/1‑139 (MQ=255)                                                                                                                              aCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGa         <  1:1100143/139‑1 (MQ=255)                                                                                                                                  gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCt                                                                                           <  2:622367/53‑1 (MQ=38)                                                                                                                                  gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCt                                                                                           >  1:622367/1‑53 (MQ=38)                                                                                                                                  gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGcc     >  1:210146/1‑139 (MQ=255)                                                                                                                                    tgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                           <  1:7536/67‑1 (MQ=255)                                                                                                                                    tgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                           >  2:7536/1‑67 (MQ=255)                                                                                                                                     gcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCttt  <  1:1400375/139‑1 (MQ=255)                                                                                                                                     gcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCttt  <  1:640921/139‑1 (MQ=255)                                                                                                                                      cGGTT‑GATGTTATTTTTTTCCTGAAGCGG‑AAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                 >  1:106624/1‑90 (MQ=255)                                                                                                                                      cGGTT‑GATGTTATTTTTTTCCTGAAGCGG‑AAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                 <  2:106624/90‑1 (MQ=255)                                                                                                                                           |                                                                                                                                      CTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTT  >  NC_000913/582010‑582280

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: ‑