Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,605,286	(A)6→5	coding (264/291 nt)	b1518 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,605,281	0	A	.	98.7%	345.8 / ‑1.6	80	coding (259/291 nt)	b1518	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (40/38);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACGGGAGTAA  >  NC_000913/1605189‑1605365                                                                                             |                                                                                     ataAAGATGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTa                                                                                      >  2:890618/1‑93 (MQ=255) ataAAGATGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTa                                                                                      >  2:539127/1‑93 (MQ=255) ataAAGATGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTa                                                                                      >  2:3046698/1‑93 (MQ=255) ataAAGATGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTa                                                                                      <  2:1816357/93‑1 (MQ=255)  taAAGATGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaa                                                                                     >  2:1462860/1‑93 (MQ=255)   aaaGATGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaa                                                                                    <  1:3506740/93‑1 (MQ=255)    aaGATGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaaa                                                                                   >  2:4629295/1‑93 (MQ=255)    aaGATGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaaa                                                                                   >  2:3324431/1‑93 (MQ=255)     aGATGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaaaa                                                                                  <  2:1018521/93‑1 (MQ=255)      gATGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaaaac                                                                                 >  1:4522985/1‑92 (MQ=255)       aTGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaaaacg                                                                                >  2:3039820/1‑91 (MQ=255)       aTGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaaaacg                                                                                >  1:2905036/1‑91 (MQ=255)        tGAAGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaaaacgt                                                                               <  1:2919166/93‑4 (MQ=255)          aaGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCt                                                                            <  2:1354732/93‑1 (MQ=255)          aaGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCt                                                                            >  2:2995820/1‑93 (MQ=255)          aaGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCt                                                                            <  2:2670512/93‑1 (MQ=255)          aaGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCt                                                                            >  1:2173/1‑93 (MQ=255)           aGACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTg                                                                           >  2:2588111/1‑93 (MQ=255)            gACGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGt                                                                          >  1:990734/1‑93 (MQ=255)              cGCAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCG‑AGAAAACGTCTGTTc                                                                        >  1:2677421/1‑93 (MQ=255)                cAGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCaa                                                                      <  2:3445225/93‑1 (MQ=255)                 aGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAAt                                                                     <  2:4023916/93‑1 (MQ=255)                 aGTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAAt                                                                     <  2:1892868/93‑1 (MQ=255)                  gTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATg                                                                    >  2:2956349/1‑93 (MQ=255)                  gTGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATg                                                                    >  1:3160439/1‑93 (MQ=255)                   tGGCGCTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATgg                                                                   <  1:4530207/93‑1 (MQ=255)                     gcgcTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGtt                                                                 <  2:181399/93‑1 (MQ=255)                     gcgcTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGtt                                                                 <  2:1548706/93‑1 (MQ=255)                       gctCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTg                                                               <  1:4520000/91‑1 (MQ=255)                       gctCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTg                                                               >  1:1832118/3‑93 (MQ=255)                       gctCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTg                                                               >  2:1472691/3‑93 (MQ=255)                        ctCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGa                                                              >  2:1816339/2‑93 (MQ=255)                        ctCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAAAGGTTTGa                                                              <  2:3164190/92‑1 (MQ=255)                         tCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGAt                                                             <  2:865006/93‑1 (MQ=255)                         tCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGAt                                                             <  1:2367570/93‑1 (MQ=255)                           cATAAACCCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGc                                                           <  1:2987351/93‑1 (MQ=255)                              aaaaCCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGt                                                        <  2:2253344/93‑1 (MQ=255)                                  ccACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAgg                                                    >  2:2123616/1‑93 (MQ=255)                                   cACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGc                                                   >  1:2342685/1‑93 (MQ=255)                                   cACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGc                                                   <  1:198098/93‑1 (MQ=255)                                    aCGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCg                                                  <  1:581409/93‑1 (MQ=255)                                     cGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGa                                                 <  2:2935492/93‑1 (MQ=255)                                     cGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGa                                                 >  1:3041636/1‑93 (MQ=255)                                      gCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGaa                                                <  1:2492168/93‑1 (MQ=255)                                      gCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGaa                                                <  1:3556966/93‑1 (MQ=255)                                      gCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGaa                                                <  2:3300943/93‑1 (MQ=255)                                        cccACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAAtt                                              <  1:289907/93‑1 (MQ=255)                                        cccACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAAtt                                              <  1:4459739/93‑1 (MQ=255)                                         ccACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAAttt                                             <  2:1505181/93‑1 (MQ=255)                                           acaaCAAATCCTGTGTCGCGAATCTGGAATCTTTAATGGCCGGGCCGC‑TTAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTAt                                           <  2:4076/90‑1 (MQ=255)                                           aCTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTAt                                           >  1:1329276/1‑93 (MQ=255)                                            cTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATc                                          >  1:606250/1‑93 (MQ=255)                                              aCAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCaa                                        <  1:3537117/93‑1 (MQ=255)                                               cAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAAt                                       >  2:1421248/1‑93 (MQ=255)                                               cAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAAt                                       <  1:2697719/93‑1 (MQ=255)                                               cAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAAt                                       >  1:1143374/1‑93 (MQ=255)                                                aaaaCCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAAtt                                      <  2:3585299/93‑1 (MQ=255)                                                aaaaCCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAAtt                                      >  2:4200642/1‑93 (MQ=255)                                                 aaaCCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAAttt                                     >  1:1258631/1‑93 (MQ=255)                                                  aaCCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAAtttt                                    >  1:2795380/1‑93 (MQ=255)                                                      tgtgTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGCATTTATCAATTTTATCt                                >  1:925156/1‑93 (MQ=255)                                                      tgtgTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCt                                <  2:2886260/93‑1 (MQ=255)                                                       gtgtCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTa                               >  2:2469661/1‑93 (MQ=255)                                                        tgtCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTAc                              <  2:2803724/93‑1 (MQ=255)                                                         gtCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACa                             >  2:1629620/1‑93 (MQ=255)                                                          tCGCGAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACaa                            <  1:177914/93‑1 (MQ=255)                                                              gAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTgg                        <  1:2533114/93‑1 (MQ=255)                                                              gAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTgg                        >  2:4537349/1‑93 (MQ=255)                                                              gAAACTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAAGTgg                        >  1:4678999/1‑93 (MQ=255)                                                               aaaCTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTggg                       <  2:1488693/93‑1 (MQ=255)                                                               aaaCTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGCGAGGCGAATTTATCAATTTTATCTACAATTggg                       >  1:1067139/1‑93 (MQ=255)                                                                aaCTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTgggg                      >  1:1167193/1‑93 (MQ=255)                                                                 aCTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGt                     >  2:4224391/1‑93 (MQ=255)                                                                  cTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTa                    <  2:2628856/93‑1 (MQ=255)                                                                  cTGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTa                    >  2:3346989/1‑93 (MQ=255)                                                                   tGGAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTaa                   >  1:3678794/1‑93 (MQ=255)                                                                    ggAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAAc                  >  1:2449369/1‑93 (MQ=255)                                                                    ggAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAAc                  <  1:2716453/93‑1 (MQ=255)                                                                    ggAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAAc                  >  1:2183638/1‑93 (MQ=255)                                                                    ggAATCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAAc                  >  1:1552789/1‑93 (MQ=255)                                                                        tCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACgcgc              >  1:3524988/1‑93 (MQ=255)                                                                        tCTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACgcgc              >  1:1815106/1‑93 (MQ=255)                                                                         cTTTAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCt             >  2:1759002/1‑93 (MQ=255)                                                                          tttAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTg            >  2:2185958/1‑93 (MQ=255)                                                                           ttAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGa           >  1:402301/1‑93 (MQ=255)                                                                           ttAATGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGa           <  2:1529934/93‑1 (MQ=255)                                                                             aaTGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACg         >  1:2365625/1‑93 (MQ=255)                                                                             aaTGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACg         <  2:2858846/93‑1 (MQ=255)                                                                              aTGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACgg        <  2:758873/93‑1 (MQ=255)                                                                               tGACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACggg       >  2:724866/1‑93 (MQ=255)                                                                                gACCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACGGGa      <  1:1506487/93‑1 (MQ=255)                                                                                 aCCGGGCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACGGGAg     <  2:4181625/93‑1 (MQ=255)                                                                                    gggCCGCGT‑AAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACGGGAGTaa  >  1:177456/1‑93 (MQ=255)                                                                                             |                                                                                     ATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACGGGAGTAA  >  NC_000913/1605189‑1605365

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑