Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,620,134	C→T	8.1%	A163A (GCG→GCA)	rnz ←	Ribonuclease Z

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,620,134	0	C	T	8.1%	202.2 / 4.7	69	A163A (GCG→GCA)	rnz	Ribonuclease Z
Reads supporting (aligned to +/- strand):  ref base C (22/41);  new base T (3/3);  total (25/44)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AAATTGCTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGATCATCGATTTCAATAT  >  USA300TCH1516_ALE/1619998‑1620270                                                                                                                                         |                                                                                                                                         aaaTTGCTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTATTGATATTTAGGACCTGGTTCATGACCAATGTTTTTTAACGCtt                                                                                                                                      >  2:627340/1‑141 (MQ=255)       cTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAATGCTTCAACat                                                                                                                                >  2:725935/1‑141 (MQ=255)       cTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAATGCTTCAACat                                                                                                                                >  2:725781/1‑141 (MQ=255)       cTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACat                                                                                                                                <  1:301623/141‑1 (MQ=255)          cAACTGGCCCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACattga                                                                                                                             <  1:756707/141‑1 (MQ=255)                 aCCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTAcc                                                                                                                      <  2:110038/141‑1 (MQ=255)                   ccTGTTTGGACTCATCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGGATTGTAACGCTTCAACATTGATTGTACCTg                                                                                                                    <  2:602886/141‑1 (MQ=255)                        ttGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATAAATCCAGAGAGTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGtt                                                                                                               <  1:210629/141‑1 (MQ=255)                          ggACTGAACTCTAAATTCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGACATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTCCCTGTAGTTTc                                                                                                             <  1:797500/141‑1 (MQ=255)                            aCTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAg                                                                                                           >  2:633564/1‑141 (MQ=255)                            aCTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAg                                                                                                           >  2:574908/1‑141 (MQ=255)                              tCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGt                                                                                                         <  1:764700/141‑1 (MQ=255)                               cACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTg                                                                                                        <  2:334080/141‑1 (MQ=255)                               cACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTg                                                                                                        <  1:730154/141‑1 (MQ=255)                               cACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACACTGATTCTACCTGTAGTTTCAGGTg                                                                                                        <  2:670487/141‑1 (MQ=255)                               cACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATCAAAATTAAGTTCTTGATATTTAGGACCTCGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTg                                                                                                        <  1:898050/141‑1 (MQ=255)                                aCCTCTAAAATCTTTTGATTGATAAACTTGGCCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGc                                                                                                       <  2:131472/141‑1 (MQ=255)                                aCCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGc                                                                                                       <  1:453154/141‑1 (MQ=255)                                     tAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTGGTTTCAGGTGCCTTTa                                                                                                  >  2:103043/1‑141 (MQ=255)                                     tAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTa                                                                                                  >  2:487692/1‑141 (MQ=255)                                     tAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTa                                                                                                  >  1:719336/1‑141 (MQ=255)                                            ttttGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAATGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCtata                                                                                           <  1:725781/141‑1 (MQ=255)                                            ttttGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAATGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCtata                                                                                           <  1:725935/141‑1 (MQ=255)                                             tttGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTatat                                                                                          <  1:82851/141‑1 (MQ=255)                                             tttGATTGATAAACTTGACCAATATGTTCAAATGCATCACGAGATGTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTatat                                                                                          <  2:399806/141‑1 (MQ=255)                                               tgattgatAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATcc                                                                                        >  2:800362/1‑141 (MQ=255)                                                gattgatAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCa                                                                                       >  1:297462/1‑141 (MQ=255)                                                gattgatAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCGTGAAATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCa                                                                                       >  1:604803/1‑141 (MQ=255)                                                 attgatAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCat                                                                                      <  1:879139/141‑1 (MQ=255)                                                     atAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATatga                                                                                  >  2:275131/1‑141 (MQ=255)                                                     aaaaaCTTGACCATTATGTTGTGATGTAGCATGAGATTTTACTTCTTAATAGTAAGGGTCTGGTTCCAGACCAATGTTTATTAACGCTTCAACTTTGATTGGACCTGTAGTTTCAGGCGCCATTACTCTATATCCATatga                                                                                  <  2:383932/139‑1 (MQ=255)                                                      tAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATatgat                                                                                 >  2:143567/1‑141 (MQ=255)                                                              aCCCTTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATAcc                                                                         <  1:512612/141‑1 (MQ=255)                                                               ccATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCa                                                                        <  1:4032/141‑1 (MQ=255)                                                                  ttATGTTCAAATGTATCATGAGATTTTACTTCTTGATATATAGGACATGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGa                                                                     >  2:373541/1‑141 (MQ=255)                                                                          agaTGTATTATGAGATTTTAATTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTaa                                                             <  1:397729/139‑1 (MQ=255)                                                                             tGTTTCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATg                                                          <  2:578098/141‑1 (MQ=255)                                                                             tGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATTTGATGGTATACCATGATTAAGTATATg                                                          <  2:260730/141‑1 (MQ=255)                                                                               tATCATGAGATTTTACTTCTTGCTATTTAGGACCTGGTTCAAGACCAATGTTTTTGAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATgcg                                                        >  2:664641/1‑141 (MQ=255)                                                                               tATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATgcg                                                        >  2:703315/1‑141 (MQ=255)                                                                                aTCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATgcgc                                                       >  1:778913/1‑141 (MQ=255)                                                                                   aTGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCctc                                                    >  1:159375/1‑141 (MQ=255)                                                                                        aTTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAg                                               <  1:330940/141‑1 (MQ=255)                                                                                        aTTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAg                                               <  1:231696/141‑1 (MQ=255)                                                                                        aTTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAg                                               <  1:137733/141‑1 (MQ=255)                                                                                         ttttACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAATGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGt                                              <  1:120932/141‑1 (MQ=255)                                                                                            tACATCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTATATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTaaacaaaa                                           >  1:654629/1‑134 (MQ=255)                                                                                             aCTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCCGCACCTTGGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTaaaa                                          <  1:121177/141‑1 (MQ=255)                                                                                                tcttGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGAACGTGCCAGTGGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACca                                       <  1:662726/141‑1 (MQ=255)                                                                                                 tttGATATTTAGGACCTTGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTTTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACcat                                      <  2:146389/140‑1 (MQ=255)                                                                                                 tttGATATGTAGGAGCTTGTTCAACGCCATTGTCTGTTAACGCTTCAACTTTTATTCTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGACTAAGTAAATGCGCCTCTACAGTAAAACcat                                      <  2:829280/140‑1 (MQ=255)                                                                                                 cttGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACcat                                      <  2:332383/141‑1 (MQ=255)                                                                                                      tatTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCatga                                 <  1:69171/141‑1 (MQ=255)                                                                                                        tttAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGCTGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCatgatg                               >  1:895073/1‑141 (MQ=255)                                                                                                           aGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGata                            >  1:748135/1‑141 (MQ=255)                                                                                                           aGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGata                            >  2:708661/1‑141 (MQ=255)                                                                                                           aGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACAGTGAGTATGGCTGGAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGata                            >  1:436911/1‑141 (MQ=255)                                                                                                                 tGGTTCAAGACCAATGTTTTTTTTTTTTTTTACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTTTATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATTATGATATTTCaa                      >  2:471876/1‑140 (MQ=255)                                                                                                                   gTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGCTTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAgat                    <  1:861326/141‑1 (MQ=255)                                                                                                                    ttCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGTGCCTCTACAGTAAAACCATCATGATGATATGTCAgatg                   <  1:824932/141‑1 (MQ=255)                                                                                                                       aaGCCCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGatc                <  2:12157/141‑1 (MQ=255)                                                                                                                         gACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGatcat              <  1:275131/141‑1 (MQ=255)                                                                                                                          gccAATGTTTTTTAACGCTTCAACATGGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCAGCATGATGATATGTCAGATGatcatc             <  2:736275/140‑1 (MQ=255)                                                                                                                          aCCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGatcatc             <  1:143567/141‑1 (MQ=255)                                                                                                                          aCCAATGCTTTTTAACGCTTAAACATTGATTGTACCTGTAGTTTCAGGGTCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGatcctc             <  2:298960/141‑1 (MQ=255)                                                                                                                               tGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGATCATCGAttt        <  1:543406/141‑1 (MQ=255)                                                                                                                                  tttttAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGATCATCGATTTCaa     <  1:184096/141‑1 (MQ=255)                                                                                                                                   ttttAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAATAAAACCATCATGAAGATAAGTCAGATGATCATCAATTTCaaa    >  2:375885/1‑140 (MQ=255)                                                                                                                                     ttAACGCTTCGACATTGACTGTACCTGTAGTTTCACGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGATCATCGATTTCAatat  <  2:106228/141‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         AAATTGCTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCAAGACCAATGTTTTTTAACGCTTCAACATTGATTGTACCTGTAGTTTCAGGTGCCATTACTCTATATCCATATGATGGAATACCATGATTAAGTAAATGCGCCTCTACAGTAAAACCATCATGATGATATGTCAGATGATCATCGATTTCAATAT  >  USA300TCH1516_ALE/1619998‑1620270

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