Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,894,996	C→T	T115T (ACC→ACT)	yieF →	chromate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,894,996	0	C	T	97.6%	275.5 / ‑2.1	82	T115T (ACC→ACT)	yieF	chromate reductase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (40/40);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

TCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCG  >  NC_000913/3894860‑3895132                                                                                                                                         |                                                                                                                                         tCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGc                                                                                                                                       >  1:453077/1‑140 (MQ=255)  acGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCt                                                                                                                                      <  2:507647/139‑1 (MQ=255)    gTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCa                                                                                                                                    <  1:1100607/140‑1 (MQ=255)             aCGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTg                                                                                                                           >  1:1000101/1‑140 (MQ=255)              cGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGa                                                                                                                          >  2:277355/1‑140 (MQ=255)                   tggtgATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTggc                                                                                                                     >  1:481292/1‑140 (MQ=255)                       gATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGgcggcg                                                                                                                 >  2:629035/1‑140 (MQ=255)                       gATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGgcggcg                                                                                                                 <  2:254864/140‑1 (MQ=255)                        aTCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgc                                                                                                                >  1:66821/1‑140 (MQ=255)                         tCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcg                                                                                                               <  1:103067/140‑1 (MQ=255)                         tCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCAGCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcg                                                                                                               <  2:258879/140‑1 (MQ=255)                          cGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgc                                                                                                              <  2:453077/140‑1 (MQ=255)                          cGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgc                                                                                                              >  1:85392/1‑140 (MQ=255)                          cGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgc                                                                                                              <  2:85392/140‑1 (MQ=255)                           gTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgcg                                                                                                             >  2:231493/1‑140 (MQ=255)                           gTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgcg                                                                                                             >  2:164565/1‑140 (MQ=255)                           gTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgcg                                                                                                             >  2:885811/1‑140 (MQ=255)                           gTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgcg                                                                                                             >  1:944388/1‑140 (MQ=255)                           gTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGA‑CAGCTCAATGGGCGTGATTGGCGgcgcgcac                                                                                                            >  1:851105/1‑138 (MQ=255)                             cACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCt                                                                                                           >  1:818905/1‑140 (MQ=255)                             cACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGA‑CAGCTCAATGGGCGTGATTGGCGGCGCGCACTg                                                                                                          <  2:851105/140‑1 (MQ=255)                              aCGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTg                                                                                                          <  2:948749/140‑1 (MQ=255)                              aCGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTg                                                                                                          <  1:629035/140‑1 (MQ=255)                               cGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgcgc                                                                                                            <  2:598711/137‑1 (MQ=255)                               cGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgcgc                                                                                                            >  1:598711/1‑137 (MQ=255)                               cGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGt                                                                                                         >  1:614432/1‑140 (MQ=255)                                gCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgcgc                                                                                                            >  1:738446/1‑136 (MQ=255)                                gCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTc                                                                                                        >  1:1136817/1‑140 (MQ=255)                                gCCGGAATATAACTACTCGGTACCGGGGGGGCTGAAAAATGCCATCGCCTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGgcgcgcgc                                                                                                            <  2:738446/136‑1 (MQ=255)                                     aaTATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTAt                                                                                                   >  2:1047292/1‑140 (MQ=255)                                            ctactCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTgc                                                                                            >  1:924036/1‑140 (MQ=255)                                             tactCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTgcg                                                                                           <  1:231493/140‑1 (MQ=255)                                                 cGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAg                                                                                       >  2:894459/1‑140 (MQ=255)                                                   gTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGAt                                                                                     >  1:488587/1‑140 (MQ=255)                                                    tACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGAtt                                                                                    >  1:649137/1‑140 (MQ=255)                                                       cGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGAtt                                                                                    <  2:60484/137‑1 (MQ=255)                                                       cGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGAtt                                                                                    >  1:60484/1‑137 (MQ=255)                                                       cGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTg                                                                                 >  2:1178291/1‑140 (MQ=255)                                                       cGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTg                                                                                 >  2:584222/1‑140 (MQ=255)                                                             ggCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTc                                                                           <  1:584222/140‑1 (MQ=255)                                                              gCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTcc                                                                          >  2:16472/1‑140 (MQ=255)                                                              gCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTcc                                                                          >  1:871433/1‑140 (MQ=255)                                                              gCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTcc                                                                          >  1:1077486/1‑140 (MQ=255)                                                              gCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTcc                                                                          >  2:553206/1‑140 (MQ=255)                                                               cTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCt                                                                         <  1:1084705/140‑1 (MQ=255)                                                               cTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCt                                                                         <  2:614432/140‑1 (MQ=255)                                                               cTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCt                                                                         <  1:734600/140‑1 (MQ=255)                                                               cTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCt                                                                         <  2:66821/140‑1 (MQ=255)                                                                  aaaaaTGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGa                                                                      <  2:481292/140‑1 (MQ=255)                                                                  aaaaaTGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGa                                                                      <  1:777446/140‑1 (MQ=255)                                                                  aaaaaTGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGa                                                                      >  1:759149/1‑140 (MQ=255)                                                                   aaaaTGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGAt                                                                     >  2:1080745/1‑140 (MQ=255)                                                                      aTGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCa                                                                  <  2:1077486/140‑1 (MQ=255)                                                                      aTGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCa                                                                  <  1:1047292/140‑1 (MQ=255)                                                                       tGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCaa                                                                 <  1:16472/140‑1 (MQ=255)                                                                       tGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCaa                                                                 <  2:649137/140‑1 (MQ=255)                                                                       tGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCaa                                                                 <  2:871433/140‑1 (MQ=255)                                                                       tGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCaa                                                                 <  1:1080745/140‑1 (MQ=255)                                                                       tGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCaa                                                                 <  2:488587/140‑1 (MQ=255)                                                                        gCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCt                                                                         <  2:291990/131‑1 (MQ=255)                                                                        gCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCt                                                                         >  1:291990/1‑131 (MQ=255)                                                                         ccATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATg                                                               <  1:885811/140‑1 (MQ=255)                                                                         ccATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATg                                                               <  2:944388/140‑1 (MQ=255)                                                                         ccATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATg                                                               <  2:527761/140‑1 (MQ=255)                                                                                              tGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAAttt                                          >  2:112377/1‑140 (MQ=255)                                                                                              tGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGCCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAAttt                                          >  1:821191/1‑140 (MQ=255)                                                                                                ccGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTAt                                        >  2:1185987/1‑140 (MQ=255)                                                                                                  ggATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATgg                                      <  2:693320/140‑1 (MQ=255)                                                                                                   gATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATggg                                     <  1:1178291/140‑1 (MQ=255)                                                                                                   gATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATggg                                     <  2:1000101/140‑1 (MQ=255)                                                                                                        aCCGCTGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGgcggcg                                <  1:874269/140‑1 (MQ=255)                                                                                                             tGGCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGAtt                           <  1:553206/140‑1 (MQ=255)                                                                                                               gCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGgcggcg                                <  2:1062100/133‑1 (MQ=255)                                                                                                               gCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGgcggcg                                >  1:1062100/1‑133 (MQ=255)                                                                                                               gCCGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCa                         >  2:233005/1‑140 (MQ=255)                                                                                                                ccGGTAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAg                        <  2:689215/140‑1 (MQ=255)                                                                                                                    tAAACCGGTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACa                    >  1:1070024/1‑140 (MQ=255)                                                                                                                           gTATTGATTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGa             >  1:438771/1‑140 (MQ=255)                                                                                                                                 aTTCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGCGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCa       >  1:458982/1‑140 (MQ=255)                                                                                                                                   tCAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCaaa     >  1:241964/1‑140 (MQ=255)                                                                                                                                    cAGACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAAc    <  1:1185987/140‑1 (MQ=255)                                                                                                                                      gACTAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCg  <  2:1070024/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         TCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCG  >  NC_000913/3894860‑3895132

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