Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,622,965	T→C	C1258C (TGT→TGC)	rhsB →	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,622,965	0	T	C	98.9%	296.7 / ‑4.5	87	C1258C (TGT→TGC)	rhsB	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (45/41);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

CTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGTACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCGCTGGTTTAA  >  NC_000913/3622828‑3623101                                                                                                                                          |                                                                                                                                         cTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCAc                                                                                                                                        >  2:111553/1‑140 (MQ=255) cTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCAc                                                                                                                                        >  2:149025/1‑140 (MQ=255)        ggCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGt                                                                                                                                 >  2:436038/1‑140 (MQ=255)        ggCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGt                                                                                                                                 <  2:2467/140‑1 (MQ=255)            atatatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCttt                                                                                                                             >  1:805147/1‑140 (MQ=255)            atatatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCttt                                                                                                                             >  2:551576/1‑140 (MQ=255)              atatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGc                                                                                                                           >  2:540837/1‑140 (MQ=255)               tatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGTACTAATGGTCTTTGCg                                                                                                                          <  1:192542/140‑1 (MQ=255)               tatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCg                                                                                                                          <  1:463511/140‑1 (MQ=255)                atCTCTCAGGATCCGATTGGACTGAA‑GGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGc                                                                                                                        <  2:86024/140‑1 (MQ=255)                   aCTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGcac                                                                                                                      >  1:644501/1‑140 (MQ=255)                      cAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACt                                                                                                                   >  2:205131/1‑140 (MQ=255)                       aGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTa                                                                                                                  >  1:616883/1‑140 (MQ=255)                       aGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTa                                                                                                                  >  1:641025/1‑140 (MQ=255)                       aGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTa                                                                                                                  >  1:917641/1‑140 (MQ=255)                       aGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTa                                                                                                                  >  1:946632/1‑140 (MQ=255)                       aGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTa                                                                                                                  >  1:377133/1‑140 (MQ=255)                        ggATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTAt                                                                                                                 <  1:205131/140‑1 (MQ=255)                         gATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGc                                                                                                                           <  2:20572/129‑1 (MQ=255)                         gATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGc                                                                                                                           >  1:20572/1‑129 (MQ=255)                         gATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATc                                                                                                                <  1:149025/140‑1 (MQ=255)                          aTCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATcc                                                                                                               <  2:8462/140‑1 (MQ=255)                           tCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCg                                                                                                              <  1:341723/140‑1 (MQ=255)                           tCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCg                                                                                                              <  2:917641/140‑1 (MQ=255)                             cGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCt                                                                                                            <  1:436038/140‑1 (MQ=255)                              gATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCt                                                                                                            >  1:99429/1‑139 (MQ=255)                              gATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCt                                                                                                            <  2:99429/139‑1 (MQ=255)                               aTTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGca                                                                                                                       >  1:603701/1‑127 (MQ=255)                               aTTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGca                                                                                                                       <  2:603701/127‑1 (MQ=255)                                 tGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGc                                                                                                                        <  2:708932/124‑1 (MQ=255)                                 tGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGc                                                                                                                        >  1:708932/1‑124 (MQ=255)                                  ggACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGAt                                                                                                       <  1:979298/140‑1 (MQ=255)                                   gACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATa                                                                                                      >  2:268563/1‑140 (MQ=255)                                              ggATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATaa                                                                                           >  2:385013/1‑140 (MQ=255)                                                 tGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAAtt                                                                                        <  2:805147/140‑1 (MQ=255)                                                 tGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAAtt                                                                                        <  2:644501/140‑1 (MQ=255)                                                 tGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAAtt                                                                                        <  2:378995/140‑1 (MQ=255)                                                   gAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCCCTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTg                                                                                      <  2:430818/140‑1 (MQ=255)                                                    aaCCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGa                                                                                     <  1:820884/140‑1 (MQ=255)                                                      ccTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGAtt                                                                                   <  1:111553/140‑1 (MQ=255)                                                      ccTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGAtt                                                                                   <  2:616883/140‑1 (MQ=255)                                                          tATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCtat                                                                               >  1:689157/1‑140 (MQ=255)                                                             ggATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACCCTATCCGCTATGATAGGTCCAGATAAATTTGATTCTataga                                                                            <  1:551576/140‑1 (MQ=255)                                                                 atCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGc                                                                        >  2:139221/1‑140 (MQ=255)                                                                  tCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCa                                                                       >  1:920767/1‑140 (MQ=255)                                                                  tCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCa                                                                       >  2:712029/1‑140 (MQ=255)                                                                  tCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCa                                                                       <  2:689157/140‑1 (MQ=255)                                                                  tCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCa                                                                       >  1:532761/1‑140 (MQ=255)                                                                   cAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCAc                                                                      >  2:537040/1‑140 (MQ=255)                                                                   cAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCAc                                                                      <  2:159266/140‑1 (MQ=255)                                                                         aaTCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCata                                                                >  2:106224/1‑140 (MQ=255)                                                                         aaTCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCata                                                                >  2:825690/1‑140 (MQ=255)                                                                         aaTCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCata                                                                >  2:244566/1‑140 (MQ=255)                                                                                    aGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTtaa                                                     <  1:825690/140‑1 (MQ=255)                                                                                          cGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATTATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGcc                                                         >  1:330025/1‑130 (MQ=255)                                                                                          cGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATTATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGcc                                                         <  2:330025/130‑1 (MQ=255)                                                                                          cGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTaaataaaa                                               >  1:409498/1‑140 (MQ=255)                                                                                          cGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTaaataaaa                                               >  1:399893/1‑140 (MQ=255)                                                                                            aCCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTaaataaaata                                             <  1:907994/140‑1 (MQ=255)                                                                                              ccccTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAataaaataaa                                           <  2:578396/140‑1 (MQ=255)                                                                                                ccTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCGAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAaataaata                                         >  1:684057/1‑140 (MQ=255)                                                                                                       ttATCTATGTGGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATc                                  >  1:45350/1‑140 (MQ=255)                                                                                                              tgtgGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGc                           >  2:412558/1‑140 (MQ=255)                                                                                                               gtgGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCg                          >  1:555029/1‑140 (MQ=255)                                                                                                               gtgGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCg                          <  2:641025/140‑1 (MQ=255)                                                                                                               gtgGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCg                          <  2:946632/140‑1 (MQ=255)                                                                                                                tgGGAGGATGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGa                         >  2:161777/1‑140 (MQ=255)                                                                                                                  ggaggaTGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCg                          >  1:329078/1‑137 (MQ=255)                                                                                                                  ggaggaTGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCg                          <  2:329078/137‑1 (MQ=255)                                                                                                                   gaggaTGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGa                      >  1:878660/1‑140 (MQ=255)                                                                                                                    aggaTGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCg                          <  2:502206/135‑1 (MQ=255)                                                                                                                    aggaTGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCg                          >  1:502206/1‑135 (MQ=255)                                                                                                                    aggaTGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGAt                     >  1:344358/1‑140 (MQ=255)                                                                                                                        tGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAag                 >  1:875728/1‑140 (MQ=255)                                                                                                                        tGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAag                 >  1:741629/1‑140 (MQ=255)                                                                                                                        tGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAag                 >  2:529851/1‑140 (MQ=255)                                                                                                                        tGCAAAATCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAag                 >  1:452265/1‑140 (MQ=255)                                                                                                                            aaaTCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGtt             <  2:741629/140‑1 (MQ=255)                                                                                                                            aaaTCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGtt             >  2:568727/1‑140 (MQ=255)                                                                                                                             aaTCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTc            <  1:412558/140‑1 (MQ=255)                                                                                                                              aTCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCg           <  2:612770/140‑1 (MQ=255)                                                                                                                              aTCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCg           <  1:529851/140‑1 (MQ=255)                                                                                                                              aTCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCg           <  2:787285/140‑1 (MQ=255)                                                                                                                               tCGGGGGCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCGc          <  1:537040/140‑1 (MQ=255)                                                                                                                                    ggCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCGCTGGtt     <  1:161777/140‑1 (MQ=255)                                                                                                                                     gCATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCGCTGGttt    <  2:878660/140‑1 (MQ=255)                                                                                                                                      cATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCGCTGGTTTa   <  2:875728/140‑1 (MQ=255)                                                                                                                                      cATGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCGCTGGTTTa   <  1:385013/140‑1 (MQ=255)                                                                                                                                       aTGCACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCGCTGGTTTaa  <  2:344358/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         CTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGTACTAATGGTCTTTGCGGCACACTATCCGCTATGATAGGTCCAGATAAATTTGATTCTATAGATAGCACCGCATATGACGCCTTAAATAAAATAAATAGCCAATCTATTTGCGAAGATAAAGAGTTCGCTGGTTTAA  >  NC_000913/3622828‑3623101

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