Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	693,258	T→A	8.3%	I65F (ATT→TTT)	USA300HOU_0636 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	693,258	0	T	A	8.3%	178.8 / 7.4	64	I65F (ATT→TTT)	USA300HOU_0636	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base A (2/3);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

GGTCATGTTAAATCACTAATACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATACGACTAGCAACATAGTTTAAAATAGGTTCACTC  >  CP000730/693120‑693397                                                                                                                                           |                                                                                                                                            ggTCATGTTAAATCACTAATACTATATGTTAAACATCAATATCGCCAATATAAATCCAATCTGCCACAGTCACTGTGAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAAtt                                                                                                                                            <  2:78777/140‑1 (MQ=255)      tGTTAAATCACTAATACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATCGACAAAGGTCAGCGTAAATTTTAgg                                                                                                                                       <  2:533642/140‑1 (MQ=255)      tGTTAAATCACTAATACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAgg                                                                                                                                       <  1:96175/140‑1 (MQ=255)        ttAAATCACTCATACTATAGGTTAAACATCAATATCGCCAATATAAATACACCCATCCACCTTCACTTGTAGCATTCTGGCATGTCCTGCTTCACTCAACCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGcc                                                                                                                                     <  2:532159/140‑1 (MQ=255)            aTCACTAATACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGAcc                                                                                                                                 <  2:556595/140‑1 (MQ=255)             tCACTAATACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGTCCTAAAATTTACGCTGACCt                                                                                                                                <  1:440342/140‑1 (MQ=255)             tCACTAATACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTAAAATTTACGCTGACCt                                                                                                                                <  1:592084/140‑1 (MQ=255)             tCACTAATACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTAAAATTTACGCTGACCt                                                                                                                                <  1:592087/140‑1 (MQ=255)             tCACTAATACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTACTACTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATGTTAGGCCGACCt                                                                                                                                >  1:222884/1‑140 (MQ=255)                   aTACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTACAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTAcc                                                                                                                          >  2:141748/1‑140 (MQ=255)                   aTACTATATGTTAAACATCAATATCGCCAATATAAATACAATATTACAACTTCACATGTCGCATAATGCCATGTCCTGCTCCACTCAAGCACATACACAATAGACAAAAACAAGCGTAAAATATAGGCCGACCTTTTAcc                                                                                                                          >  2:62935/1‑140 (MQ=255)                    tACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCg                                                                                                                         <  2:176052/140‑1 (MQ=255)                       tataTGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAt                                                                                                                      <  1:556568/140‑1 (MQ=255)                        tcatGTTAAACATCAATATCGCCAAGATAAATACAATCTTCCAAATTCACTTGTAGCATGCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAtt                                                                                                                     <  2:140578/138‑1 (MQ=255)                        atatGTTAAACATCAATATCGCCAATATGAATACAATCTTCCAAATTCACTTGTAGCATTCTGTAATGTCCTGCTTCACTCAAGCTCATACACAATGAACAAAGGTCAGCGAAAACTTTAGGCCGACCTTTGACAGTagt                                                                                                                     >  2:377560/1‑138 (MQ=255)                        atatGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAAACAAAGGTCAGAGTAAATTTTAGGCCGACCTGTTACCGTAtt                                                                                                                     >  2:377561/1‑140 (MQ=255)                         tatGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATCCACAATAGACAAAGGACAACGTAAACTCTTGTCCTACCATGTACCGTAatt                                                                                                                    >  1:523813/1‑140 (MQ=255)                            gTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCt                                                                                                                 >  1:675042/1‑140 (MQ=255)                            gTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATAAACAATAGAAAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACGGTATTTTCt                                                                                                                 >  1:520573/1‑140 (MQ=255)                                  cATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACg                                                                                                           <  2:125844/140‑1 (MQ=255)                                   ctCAATCTCGCCAATATAAATACCATCTTCCAAATTCACTTGTAGCATTCTGCCAGGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGt                                                                                                          <  2:173723/139‑1 (MQ=255)                                   aTCAATATCGCCAATATAAATACAATCTTCCAAATTCATTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGt                                                                                                          <  1:100806/140‑1 (MQ=255)                                   aTCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGt                                                                                                          >  1:39489/1‑140 (MQ=255)                                    tCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTa                                                                                                         >  2:551065/1‑140 (MQ=255)                                    tCAATATCGACAATATAAATACAATCTTCCAAATTGACTTGTAGCATTCTGCCATGTCCTGCTTCCCTCAAGCTCATACAAAATAGCCAAAGGTCAGCGTAAACTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTa                                                                                                         >  2:232654/1‑140 (MQ=255)                                      aaTATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACGGTATTTTCTTTAACGTAAc                                                                                                       >  2:553533/1‑140 (MQ=255)                                           cGCCAATATAAATACGATCTTCCAATTTCTCTTGTAGCATTCTGCCTTGTCCTGCTTCTCTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCa                                                                                                  <  1:623857/140‑1 (MQ=255)                                                atatAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAGTTTCTTTAACGTAACTATCATTaaa                                                                                             <  1:503782/140‑1 (MQ=255)                                                  ataAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAAtt                                                                                           >  1:587059/1‑140 (MQ=255)                                                   taAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAAttt                                                                                          <  1:485307/140‑1 (MQ=255)                                                      aTACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGa                                                                                       >  1:253472/1‑140 (MQ=255)                                                       tACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACATAACTATCATTAAATTTAGAc                                                                                      >  1:467375/1‑140 (MQ=255)                                                       tACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTACTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTATGCCGACCTTTTACCGTATTTTCTTTTACTTAACTATCATTAATTTTTGAc                                                                                      >  2:64454/1‑140 (MQ=255)                                                             cTTCCAAATTCACTTGTAGGATTCTGCCATGTCCTGCTTCACTCAAGCTCATACGCAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATGTTCTTTAACGTAACTATCATTAAATTTAGACTGTAtt                                                                                >  1:471563/1‑140 (MQ=255)                                                                ccAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCt                                                                             >  2:285600/1‑140 (MQ=255)                                                                         cTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTGATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACTGTATTTTCTTTAACGTAACTATCATTAAATTTACCCTGTATTTCTTTAGCTaaa                                                                    <  1:471359/140‑1 (MQ=255)                                                                           tGTAGCATTCTGCCATGTCCTGCTTCACTCAAGTTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGt                                                                  <  1:639620/140‑1 (MQ=255)                                                                           tGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGt                                                                  <  1:572418/140‑1 (MQ=255)                                                                           tGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGt                                                                  <  1:361629/140‑1 (MQ=255)                                                                           tGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGt                                                                  <  1:659129/140‑1 (MQ=255)                                                                                cATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTa                                                             <  2:181988/140‑1 (MQ=255)                                                                                 aTTCTTCCATGTCCTGCCTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTaa                                                            <  1:203949/140‑1 (MQ=255)                                                                                 aTTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGCTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTaa                                                            <  1:131332/140‑1 (MQ=255)                                                                                     tGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTAATTAAGGaa                                                        >  2:664397/1‑140 (MQ=255)                                                                                      gCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGaaa                                                       <  1:318007/140‑1 (MQ=255)                                                                                       ccATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAAc                                                      <  1:311594/140‑1 (MQ=255)                                                                                         tggTCCTGCTTCACTCAAGCACATACACTATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACtt                                                    <  1:652725/138‑1 (MQ=255)                                                                                                     acccAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAAt                                        <  2:523813/137‑1 (MQ=255)                                                                                                     aCTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAAt                                        >  1:364894/1‑140 (MQ=255)                                                                                                     aCTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAAt                                        <  1:660741/140‑1 (MQ=255)                                                                                                          aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTtata                                   >  2:354353/1‑140 (MQ=255)                                                                                                          aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTtata                                   >  2:561859/1‑140 (MQ=255)                                                                                                           gCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATAc                                  <  1:551065/140‑1 (MQ=255)                                                                                                                 aaaatacggtaaaaGGTCGGCCTAAAATTTACGCTGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATACGACTAg                            >  2:721679/12‑140 (MQ=37)                                                                                                                  cacaAGAGTCAAAGTCCAGCGTTAATTTTAGTCCGACATTTTCCAGTATTTTATTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAATTCATTAATGAAACTTTATCGCCTTAATTTATACGACTAGc                           <  2:219654/140‑1 (MQ=255)                                                                                                                    caATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACATTATCGCCCCAATTTATACGACTAGCaa                         >  1:47034/1‑140 (MQ=255)                                                                                                                       tAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATACGACTAGCAACAt                      <  1:141748/140‑1 (MQ=255)                                                                                                                              aGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTAATTATTAGAACTATATCACCCCAATTTATACGACTAGCAACATAGTTTaa               >  1:377914/1‑140 (MQ=255)                                                                                                                                 tCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCACTTTATACGACTAGCAAAATAGTTTAAAAt            >  2:310860/1‑140 (MQ=255)                                                                                                                                 tCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATACGACTAGCAACATAGTTTAAAAt            >  2:528944/1‑140 (MQ=255)                                                                                                                                  cAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATCGCACTAGCATCATAGTTTAAAATa           >  1:130558/1‑140 (MQ=255)                                                                                                                                     cGTAACTTTTAGGCCGACCTTTTACCGTATTTTCTTTCACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATACGACTAGCAACATAGTTTAAAATAGGt        >  1:177609/1‑140 (MQ=255)                                                                                                                                     cGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATAAAACTTTATCGCCCCAATTTATACGACTAGCAACATAGTTTAAAATAGGt        >  1:205977/1‑140 (MQ=255)                                                                                                                                       tataTTTTAGGCCGACCTTTTACCGTATTTTCTTTAATGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATACGACTAGCAACATAGTTTAAAATAGGTTc      <  2:467030/137‑1 (MQ=255)                                                                                                                                          aTTTTAGGCAGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATACGACTAGCAACATAGTTTAAAATAGGTTCACt   <  1:62935/140‑1 (MQ=255)                                                                                                                                           ttttAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATACGACTAGCAACATAGTTTAAAATAGGTTCACTc  <  1:232654/140‑1 (MQ=255)                                                                                                                                           ttttAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATACGACTAGCAACATAGTTTAAAATAGGTTCACTc  <  2:222884/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                            GGTCATGTTAAATCACTAATACTATATGTTAAACATCAATATCGCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGCCCCAATTTATACGACTAGCAACATAGTTTAAAATAGGTTCACTC  >  CP000730/693120‑693397

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: ‑