Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,132	(C)8→9	pseudogene (3856/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,124	1	.	C	96.1%	278.8 / 3.4	77	pseudogene (2/1461 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (36/38);  total (37/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 = 9.96e-01

AGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTC  >  NC_000913/1593974‑1594268                                                                                                                                                        |                                                                                                                                                    aGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTg                                                                                                                                                      <  2:757886/151‑1 (MQ=255)    gTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cc‑c                                                                                                                                                >  1:379593/1‑151 (MQ=255)      cgcgTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cc‑ccc                                                                                                                                              <  2:173143/151‑1 (MQ=255)        cgTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cc‑ccccc                                                                                                                                            <  4:530236/151‑1 (MQ=255)        cgTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cc‑ccccc                                                                                                                                            <  4:790661/151‑1 (MQ=255)            gAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cc‑ccccc‑cca                                                                                                                                        <  4:580692/150‑3 (MQ=255)              aCGTCATCAGGCCATTGTTGCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cc‑ccccc‑ccaggc                                                                                                                                     <  2:1613987/151‑6 (MQ=255)                   atcaGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAAc                                                                                                                                 >  4:848260/1‑151 (MQ=255)                       ggCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTc                                                                                                                             >  4:52005/1‑151 (MQ=255)                       ggCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATGCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGtt                                                                                                                              <  2:936215/150‑1 (MQ=255)                           aTTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCatat                                                                                                                         <  4:863740/151‑1 (MQ=255)                           aTTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCatat                                                                                                                         <  4:886260/151‑1 (MQ=255)                              gTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCa                                                                                                                       <  2:1496111/150‑1 (MQ=255)                              gTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCAc                                                                                                                      >  4:430148/1‑151 (MQ=255)                               tttCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCAcc                                                                                                                     >  1:32260/1‑151 (MQ=255)                                ttCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCAcc                                                                                                                     >  4:581631/1‑150 (MQ=255)                                   cGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCa                                                                                                                 <  4:1553580/151‑1 (MQ=255)                                     gTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAgt                                                                                                               <  4:431093/151‑1 (MQ=255)                                       ataGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAgtg                                                                                                              <  1:331451/150‑1 (MQ=255)                                         aGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAgtggtg                                                                                                           <  4:1100083/151‑1 (MQ=255)                                         aGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAgtggt                                                                                                            <  3:344590/150‑1 (MQ=255)                                          gTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAgtggtg                                                                                                           <  2:404616/150‑1 (MQ=255)                                              ccGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACAtt                                                                                                      >  4:1212009/1‑151 (MQ=255)                                               cGTTAACGGTCAGTCCATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACAttt                                                                                                     <  2:628393/151‑1 (MQ=255)                                                 ttAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                                                                                   >  3:42575/1‑151 (MQ=255)                                                  tAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTccc                                                                                                  >  4:918759/1‑151 (MQ=255)                                                   aaCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTccc                                                                                                  <  2:492740/150‑1 (MQ=255)                                                   aaCGGTCAGTACATTCCCCGGGTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCg                                                                                                 >  1:64156/1‑151 (MQ=255)                                                    aCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGa                                                                                                <  1:318125/151‑1 (MQ=255)                                                          aGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcac                                                                                          >  1:448934/1‑151 (MQ=255)                                                           gTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcac                                                                                          >  2:1168292/1‑150 (MQ=255)                                                             aCATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGACCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACAcc                                                                                       <  2:1199487/151‑1 (MQ=255)                                                             aCATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                                                        >  1:32928/1‑151 (MQ=255)                                                             aCATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                                                        <  4:1245676/150‑1 (MQ=255)                                                              cATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                                                      <  3:348871/151‑1 (MQ=255)                                                                 tCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCt                                                                                    <  3:158385/150‑1 (MQ=255)                                                                  ccccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCtt                                                                                   <  2:237974/150‑2 (MQ=255)                                                                  ccccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                                                                 <  1:111143/151‑1 (MQ=255)                                                                   cccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                                                                 <  3:330450/151‑1 (MQ=255)                                                                       gTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                                                                              <  3:302761/150‑1 (MQ=255)                                                                        tttCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATgg                                                                             >  2:125692/1‑150 (MQ=255)                                                                          tCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGat                                                                           <  2:1341100/150‑1 (MQ=255)                                                                           ccccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGatat                                                                         <  4:579248/151‑1 (MQ=255)                                                                              cccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGatat                                                                         >  4:35207/1‑148 (MQ=255)                                                                              cccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTaa                                                                      <  2:675103/151‑1 (MQ=255)                                                                                    gCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGtt                                                                <  2:379132/151‑1 (MQ=255)                                                                                      tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGtt                                                                <  1:348804/149‑1 (MQ=255)                                                                                      tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg                                                               >  4:326716/1‑150 (MQ=255)                                                                                      tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg                                                               >  4:165195/1‑150 (MQ=255)                                                                                       aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg                                                               <  2:625659/149‑1 (MQ=255)                                                                                       aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg                                                               <  2:1677396/149‑1 (MQ=255)                                                                                       aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGa                                                              >  1:135379/1‑150 (MQ=255)                                                                                          cATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACtt                                                           >  1:103060/1‑150 (MQ=255)                                                                                          cATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACtt                                                           <  4:3842/150‑1 (MQ=255)                                                                                          cATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTg                                                          >  2:836654/1‑151 (MQ=255)                                                                                           aTTTGAACCATTCCGCCATCATCCAGGTTGTTCCCGATAGTGGTTTTAGCGAAACGTAAAG‑‑CCAACCCC‑CAGGTATTACGTTTATATCACCTGCAGTGGTGACATTTCCCGACATTACACCGCTGTTATGGATATTAATGTGTTGatt                                                            <  2:430136/148‑3 (MQ=255)                                                                                              tGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa                                                       <  2:561856/150‑1 (MQ=255)                                                                                              tGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa                                                       <  2:541946/150‑1 (MQ=255)                                                                                              tGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa                                                       <  4:1546101/150‑1 (MQ=255)                                                                                                 aCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATAtt                                                   >  4:395257/1‑151 (MQ=255)                                                                                                  ccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATAtt                                                   >  2:566996/1‑150 (MQ=255)                                                                                                   cGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTaa                                                 <  2:1141346/151‑1 (MQ=255)                                                                                                    gTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGCCATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTaa                                                 >  1:124382/1‑150 (MQ=255)                                                                                                     tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAAt                                                <  1:328623/150‑1 (MQ=255)                                                                                                      gccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATgg                                              >  1:90666/1‑151 (MQ=255)                                                                                                       ccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc                                             >  1:359102/1‑151 (MQ=255)                                                                                                       ccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑CC‑CCCCCTCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc                                             >  2:1576600/1‑151 (MQ=255)                                                                                                        cgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc                                             >  1:308264/1‑150 (MQ=255)                                                                                                          ccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc                                             >  2:482589/1‑148 (MQ=255)                                                                                                          ccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc                                             >  4:381734/1‑148 (MQ=255)                                                                                                               ctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCgg                                     >  4:942810/1‑151 (MQ=255)                                                                                                                     gTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCg                                >  4:1016654/1‑150 (MQ=255)                                                                                                                       tGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGaa                              >  2:486977/1‑150 (MQ=255)                                                                                                                        gccgccGATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAg                            <  4:812401/151‑1 (MQ=255)                                                                                                                          cgccgATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCa                          <  4:1092709/151‑1 (MQ=255)                                                                                                                           gccgATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCaa                         <  4:976614/151‑1 (MQ=255)                                                                                                                            ccgATAGTGGTTTTAGCGACACGCAGTGCCCC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCaa                         <  1:63028/151‑1 (MQ=255)                                                                                                                             cgATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAAtt                       >  4:820954/1‑151 (MQ=255)                                                                                                                             cgATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAAtt                       >  4:857501/1‑151 (MQ=255)                                                                                                                              gATAGTGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAAttt                      >  4:515451/1‑151 (MQ=255)                                                                                                                                   tGGTTTTAGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccc                 >  4:498997/1‑151 (MQ=255)                                                                                                                                          aGCGACACGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACg           >  2:184709/1‑150 (MQ=255)                                                                                                                                                acGCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCcg    >  2:1560872/1‑151 (MQ=255)                                                                                                                                                  gCAGTGC‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCcgtc  <  1:277234/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                    AGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑CC‑CCCCC‑CAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTC  >  NC_000913/1593974‑1594268

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