Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,431,170	A→C	5.3%	R82R (CGT→CGG)	tsaC ←	threonylcarbamoyl‑AMP synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,431,170	0	A	C	5.3%	228.6 / 2.6	71	R82R (CGT→CGG)	tsaC	threonylcarbamoyl‑AMP synthase
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base C (2/2);  total (34/37)
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

CATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGATCGCACCCA  >  NC_000913/3431022‑3431317                                                                                                                                                     |                                                                                                                                                    cATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCAc                                                                                                                                                    <  2:1060986/150‑1 (MQ=255)    aaGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCt                                                                                                                                                 >  2:469494/1‑150 (MQ=255)    aaGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCt                                                                                                                                                 >  1:513955/1‑150 (MQ=255)     aGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTg                                                                                                                                                <  1:49803/150‑1 (MQ=255)      gCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGTGAAACCATTTTTTCCCGCTgg                                                                                                                                               <  1:463636/150‑2 (MQ=255)      gCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGc                                                                                                                                               >  2:481787/1‑150 (MQ=255)          ggCACAAAGCAACCACCAACGGGTGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACgt                                                                                                                                           >  1:918851/1‑150 (MQ=255)           gCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACgtc                                                                                                                                          >  2:217908/1‑150 (MQ=255)            cacaAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACgtca                                                                                                                                         >  1:61075/1‑150 (MQ=255)               aaaGCAACCACCAACGGGTGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACgtcagtc                                                                                                                                      <  2:918851/150‑1 (MQ=255)                       caccaACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGgt                                                                                                                              >  1:571102/1‑150 (MQ=255)                         ccaACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGgtgt                                                                                                                            <  2:1086315/150‑1 (MQ=255)                          caACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGtc                                                                                                                           <  2:61075/150‑1 (MQ=255)                             cGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGtcatc                                                                                                                        <  1:988665/150‑1 (MQ=255)                              ggATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGtcatca                                                                                                                       <  2:987748/150‑1 (MQ=255)                               gATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCaa                                                                                                                      >  1:236687/1‑150 (MQ=255)                                 tGGTCTGTTACTCGTACATCAAGCGAATCAATGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTTACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAata                                                                                                                    <  2:82811/150‑1 (MQ=255)                                   gTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACACGGTGTCATCAAtata                                                                                                                  >  2:615391/1‑150 (MQ=255)                                    tCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAg                                                                                                                 >  2:1274092/1‑150 (MQ=255)                                      ggTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAggg                                                                                                               <  1:871615/150‑1 (MQ=255)                                      ggTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAggg                                                                                                               >  2:1074837/1‑150 (MQ=255)                                       gTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGt                                                                                                              >  1:1378799/1‑150 (MQ=255)                                        tGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCTAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGtt                                                                                                             <  2:54302/150‑1 (MQ=255)                                          aCTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTa                                                                                                           <  1:1074617/150‑1 (MQ=255)                                                 cAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctgct                                                                                                    >  1:1229484/1‑150 (MQ=255)                                                 cAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctgct                                                                                                    <  1:217908/150‑1 (MQ=255)                                                     aaGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTa                                                                                                <  1:616258/150‑1 (MQ=255)                                                         gAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTGGTAAttt                                                                                            >  1:1310024/1‑150 (MQ=255)                                                                aGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCCCGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCACGTAATTTGCTGCGa                                                                                     >  2:1349084/1‑150 (MQ=255)                                                                   ggCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTa                                                                                  <  2:513955/150‑1 (MQ=255)                                                                   ggCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTa                                                                                  >  2:1075688/1‑150 (MQ=255)                                                                     cccGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTaaa                                                                                <  2:537009/150‑1 (MQ=255)                                                                       cGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAAt                                                                              <  1:565143/150‑1 (MQ=255)                                                                            aCCAGCGCGGTGTTGTCGCAGGTGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGcc                                                                         >  1:1051143/1‑150 (MQ=255)                                                                             ccAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAggg                                                                                                               >  1:1160399/1‑111 (MQ=255)                                                                             ccAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAggg                                                                                                               <  2:1160399/111‑1 (MQ=255)                                                                               aGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCCGCCCCt                                                                      >  2:1205448/1‑150 (MQ=255)                                                                                gcgcGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCtt                                                                     <  1:1075688/150‑1 (MQ=255)                                                                                 cgcgGTGTTGTCGCAGGAGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATTGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTa                                                                    <  2:1383252/150‑1 (MQ=255)                                                                                   cgGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATc                                                                  >  2:471619/1‑150 (MQ=255)                                                                                   cgGTGTTCTCGCAGGCAAGCGAAAAACAAAGGTGACAGGAACTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATc                                                                  >  2:1416611/1‑150 (MQ=255)                                                                                     gtgtTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCaa                                                                >  1:1250232/1‑150 (MQ=255)                                                                                         tgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTa                                                                                                           >  1:53326/1‑103 (MQ=255)                                                                                         tgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTa                                                                                                           <  2:53326/103‑1 (MQ=255)                                                                                         tgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCgg                                                            <  1:469494/150‑1 (MQ=255)                                                                                         tgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCgg                                                            >  1:863822/1‑150 (MQ=255)                                                                                            cGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTa                                                                                                           <  1:698748/100‑1 (MQ=255)                                                                                            cGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTa                                                                                                           >  2:698748/1‑100 (MQ=255)                                                                                              cAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCt                                                       <  1:1074837/150‑1 (MQ=255)                                                                                                ggCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTGGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGt                                                     <  2:1310024/150‑1 (MQ=255)                                                                                                      ggAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                               <  2:887673/150‑1 (MQ=255)                                                                                                       gAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACt                                              <  2:86172/150‑1 (MQ=255)                                                                                                           gACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCaa                                          >  1:1196701/1‑150 (MQ=255)                                                                                                             cAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACa                                        >  1:754198/1‑150 (MQ=255)                                                                                                                    gACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCAt                                 <  1:481787/150‑1 (MQ=255)                                                                                                                      cAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCa                               >  2:263546/1‑150 (MQ=255)                                                                                                                      cAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCa                               >  2:136981/1‑150 (MQ=255)                                                                                                                      cAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCa                               >  1:726324/1‑150 (MQ=255)                                                                                                                          aCCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCActgc                           <  2:1378799/150‑1 (MQ=255)                                                                                                                            acgtgcAGCGTGAAACCATTTTTTCCCGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGc                                                                          <  2:389329/96‑1 (MQ=37)                                                                                                                            acgtgcAGCGTGAAACCATTTTTTCCCGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGc                                                                          >  1:389329/6‑101 (MQ=37)                                                                                                                                 cAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCg                    >  1:776075/1‑150 (MQ=255)                                                                                                                                 cAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCg                    <  1:263546/150‑1 (MQ=255)                                                                                                                                   gCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCt                  <  2:726324/150‑1 (MQ=255)                                                                                                                                      ggAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATc               <  2:776075/150‑1 (MQ=255)                                                                                                                                       gAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCa              <  2:571102/150‑1 (MQ=255)                                                                                                                                       gAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCa              >  2:719258/1‑150 (MQ=255)                                                                                                                                            aaTGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGATCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGATc         <  2:754198/150‑1 (MQ=255)                                                                                                                                             aTGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGGTTAAGCTGCTCGTAATTTACTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGATCg        <  1:1274092/150‑1 (MQ=255)                                                                                                                                                   tCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGATCGCACCCa  >  1:197939/1‑150 (MQ=255)                                                                                                                                                   tCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGATCGCACCCa  <  1:136981/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGATCGCACCCA  >  NC_000913/3431022‑3431317

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