Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	401,548	(G)5→4	coding (163/261 nt)	iraP →	anti‑RssB factor, RpoS stabilzer during Pi starvation, anti‑adapter protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	401,544	0	G	.	100.0%	253.0 / NA	62	coding (159/261 nt)	iraP	anti‑RssB factor, RpoS stabilzer during Pi starvation, anti‑adapter protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (31/31);  total (31/31)

TATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCT  >  NC_000913/401411‑401668                                                                                                                                      |                                                                                                                             tATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         <  2:862112/139‑3 (MQ=255) tATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         <  1:562107/139‑3 (MQ=255) tATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         <  1:1612486/139‑3 (MQ=255) tATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         <  1:813469/139‑3 (MQ=255) tATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         <  2:302512/139‑3 (MQ=255)  aTTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcgc                                                                                                                        >  2:217994/1‑136 (MQ=255)       aGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTGTCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTAc                                                                                                                  <  1:1154528/139‑1 (MQ=255)          ttGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAg                                                                                                               >  1:1056286/1‑139 (MQ=255)            gCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTAc                                                                                                                  >  1:1195106/1‑134 (MQ=255)            gCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTAc                                                                                                                  <  2:1195106/134‑1 (MQ=255)             cccAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTa                                                                                                            >  1:310800/1‑139 (MQ=255)             cccAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTa                                                                                                            >  2:1004114/1‑139 (MQ=255)              ccAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTaa                                                                                                           <  2:1429077/139‑1 (MQ=255)                     gaagaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGa                                                                                                    <  1:642733/139‑1 (MQ=255)                                  aaGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCt                                                                                       <  1:624757/139‑1 (MQ=255)                                                   gTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgct                                                                      >  2:1109212/1‑139 (MQ=255)                                                    tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCt                                                                                       <  2:767894/121‑1 (MQ=255)                                                    tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCt                                                                                       >  1:767894/1‑121 (MQ=255)                                                    tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctg                                                                     >  2:372200/1‑139 (MQ=255)                                                            ttGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGa                                                             <  1:1058427/139‑1 (MQ=255)                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGAcgacga                                                                               <  2:258236/120‑1 (MQ=255)                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGAcgacga                                                                               >  1:258236/1‑120 (MQ=255)                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                            >  2:571939/1‑139 (MQ=255)                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                            >  2:307968/1‑139 (MQ=255)                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                            >  2:324309/1‑139 (MQ=255)                                                              ggAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                           <  1:1037136/139‑1 (MQ=255)                                                                   ttATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTa                                                      >  2:922946/1‑139 (MQ=255)                                                                   ttATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTa                                                      <  1:736629/139‑1 (MQ=255)                                                                    tATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaa                                                     >  2:343915/1‑139 (MQ=255)                                                                     aTCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                    <  2:1374873/139‑1 (MQ=255)                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         >  1:340454/1‑67 (MQ=255)                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         <  2:340454/69‑3 (MQ=255)                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaaga                                                  <  1:307968/139‑1 (MQ=255)                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaaga                                                  <  1:1290662/139‑1 (MQ=255)                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaaga                                                  <  1:1109212/139‑1 (MQ=255)                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaaga                                                  >  1:109643/1‑139 (MQ=255)                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaaga                                                  <  1:343915/139‑1 (MQ=255)                                                                         tCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                >  1:142833/1‑139 (MQ=255)                                                                             tGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTAt                                            <  2:1056286/139‑1 (MQ=255)                                                                               cAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         >  2:1462081/1‑59 (MQ=255)                                                                               cAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         <  1:1462081/61‑3 (MQ=255)                                                                                 aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                        >  2:163304/1‑139 (MQ=255)                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACg                                                        <  2:1501680/120‑1 (MQ=255)                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACg                                                        >  1:1501680/1‑120 (MQ=255)                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCa                                     >  2:375414/1‑139 (MQ=255)                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCa                                     >  1:354877/1‑139 (MQ=255)                                                                                       cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATc                                   >  2:967966/1‑138 (MQ=255)                                                                                       cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACATAAAGAAGTTATTGAAGCacc                                   <  1:967966/138‑3 (MQ=255)                                                                                           atatGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGt                              <  1:375414/139‑1 (MQ=255)                                                                                              tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         <  2:1029217/46‑3 (MQ=255)                                                                                              tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                         >  1:1029217/1‑44 (MQ=255)                                                                                              tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAg                           >  1:1519534/1‑139 (MQ=255)                                                                                              tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAg                           >  1:1305982/1‑139 (MQ=255)                                                                                                gcgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTa                         <  1:1215578/139‑1 (MQ=255)                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCa                                                                                            >  1:1266478/1‑71 (MQ=255)                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCa                                                                                            <  2:1266478/71‑1 (MQ=255)                                                                                                   cAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGAt                                                                                                   <  2:1177971/62‑1 (MQ=255)                                                                                                   cAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGAt                                                                                                   >  1:1177971/1‑62 (MQ=255)                                                                                                      aaTGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTga                                                                                     >  1:454657/1‑73 (MQ=255)                                                                                                      aaTGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTga                                                                                     <  2:454657/73‑1 (MQ=255)                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCCGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAAt               <  2:310800/139‑1 (MQ=255)                                                                                                                    gTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacga                                                                                  >  1:776096/1‑62 (MQ=255)                                                                                                                    gTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacga                                                                                  <  2:776096/62‑1 (MQ=255)                                                                                                                     ttgattgaTCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg    <  1:1407258/139‑1 (MQ=255)                                                                                                                     ttgattgaTCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg    <  1:217994/139‑1 (MQ=255)                                                                                                                       gattgatCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGAAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCt  >  1:1543043/1‑139 (MQ=255)                                                                                                                            atCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                <  1:201012/88‑1 (MQ=255)                                                                                                                            atCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                >  2:201012/1‑88 (MQ=255)                                                                                                                            atCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg    <  2:322714/132‑1 (MQ=255)                                                                                                                            atCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg    >  1:322714/1‑132 (MQ=255)                                                                                                                                 gTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaag                                                   >  2:1513486/1‑80 (MQ=255)                                                                                                                                 gTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaag                                                   <  1:1513486/80‑1 (MQ=255)                                                                                                                                 gTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCg                               <  1:164486/100‑1 (MQ=255)                                                                                                                                 gTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCg                               >  2:164486/1‑100 (MQ=255)                                                                                                                                      |                                                                                                                             TATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCT  >  NC_000913/401411‑401668

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