Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,743	A→G	D189G (GAC→GGC)	creD →	putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,743	0	A	G	98.5%	235.0 / ‑4.4	67	D189G (GAC→GGC)	creD	putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (40/26);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

TGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAT  >  NC_000913/4638608‑4638879                                                                                                                                        |                                                                                                                                         tgATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTgg                                                                                                                                      >  2:1273123/1‑140 (MQ=255)   aTTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCg                                                                                                                                    <  1:662607/140‑1 (MQ=255)   aTTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCg                                                                                                                                    <  2:119615/140‑1 (MQ=255)             gggATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAAc                                                                                                                          >  1:190397/1‑140 (MQ=255)              ggATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAAcc                                                                                                                         <  1:240277/140‑1 (MQ=255)              ggATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAAcc                                                                                                                         <  1:200575/140‑1 (MQ=255)                   cgcgTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAg                                                                                                                    >  2:935896/1‑140 (MQ=255)                             tggtgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATTCCTTTGCCTGTAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATgg                                                                                                          <  2:1462494/140‑1 (MQ=255)                               gtgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGcc                                                                                                        <  1:1273123/140‑1 (MQ=255)                               gtgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGcc                                                                                                        <  2:981852/140‑1 (MQ=255)                                  tggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTg                                                                                                     >  1:1215284/1‑140 (MQ=255)                                    gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGaa                                                                                                   >  1:868756/1‑140 (MQ=255)                                    gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGaa                                                                                                   >  2:1539594/1‑140 (MQ=255)                                    gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGaa                                                                                                   >  1:61252/1‑140 (MQ=255)                                    gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGaa                                                                                                   >  1:599671/1‑140 (MQ=255)                                    gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGaa                                                                                                   >  2:853690/1‑140 (MQ=255)                                    gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGaa                                                                                                   >  2:132902/1‑140 (MQ=255)                                      gAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAAtt                                                                                                 <  1:1539594/140‑1 (MQ=255)                                        aaGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTa                                                                                               <  2:61252/140‑1 (MQ=255)                                         aGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAAt                                                                                                  <  1:1394368/136‑1 (MQ=255)                                         aGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAAt                                                                                                  >  2:1394368/1‑136 (MQ=255)                                          gcgcCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAg                                                                                             <  1:73215/140‑1 (MQ=255)                                                 aaGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTAcc                                                                                      <  1:375944/140‑1 (MQ=255)                                                           aacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCttt                                                                            >  2:1320714/1‑140 (MQ=255)                                                           aacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCttt                                                                            <  1:132902/140‑1 (MQ=255)                                                           aacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCttt                                                                            >  2:1422934/1‑140 (MQ=255)                                                           aacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCttt                                                                            <  2:190397/140‑1 (MQ=255)                                                            acggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTc                                                                           >  2:362082/1‑140 (MQ=255)                                                                cgcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATACCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCtgtg                                                                       >  1:1140572/1‑140 (MQ=255)                                                                 gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTgtgg                                                                      <  1:853690/140‑1 (MQ=255)                                                                     aCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCATAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGcc                                                                  >  1:356107/1‑140 (MQ=255)                                                                       cATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTg                                                                <  1:959913/140‑1 (MQ=255)                                                                         ttGAACCCGGCGCCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCAGCGATCTTTCTGTGGTGCCTggc                                                              <  2:883535/140‑1 (MQ=255)                                                                         ttGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTggc                                                              <  1:1128892/140‑1 (MQ=255)                                                                          tGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTggcg                                                             >  2:438570/1‑140 (MQ=255)                                                                              cccGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTggcgggcg                                                         <  2:868756/140‑1 (MQ=255)                                                                                 ggCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTaa                                                      >  1:184486/1‑140 (MQ=255)                                                                                 ggCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTaa                                                      >  2:418493/1‑140 (MQ=255)                                                                                   cgcaggggTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTAACTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATa                                                    >  2:250921/5‑140 (MQ=255)                                                                                   cACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATa                                                    <  2:1291868/140‑1 (MQ=255)                                                                                            agagCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGa                                           >  1:671509/1‑140 (MQ=255)                                                                                              agCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGAcc                                         <  1:438570/140‑1 (MQ=255)                                                                                               gCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCt                                        <  1:418493/140‑1 (MQ=255)                                                                                                cAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCtt                                       >  2:1251795/1‑140 (MQ=255)                                                                                                   ggCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAAc                                    >  1:63812/1‑140 (MQ=255)                                                                                                   ggCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCAGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAAc                                    >  2:972248/1‑140 (MQ=255)                                                                                                       ggCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGGGTAATAGCGAAATGACCTTAAACAGc                                >  2:1685754/1‑140 (MQ=255)                                                                                                        gCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCa                               >  2:194331/1‑140 (MQ=255)                                                                                                        gCAGGGCGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCa                               >  1:1512552/1‑140 (MQ=255)                                                                                                              cGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGc                         >  1:550954/1‑140 (MQ=255)                                                                                                              cGTGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGc                         >  1:1649453/1‑140 (MQ=255)                                                                                                                tGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCg                       <  1:1290504/140‑1 (MQ=255)                                                                                                                  cATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCa                     <  1:1422934/140‑1 (MQ=255)                                                                                                                      tCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATccc                 >  2:144512/1‑140 (MQ=255)                                                                                                                      tCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATccc                 >  2:1089009/1‑140 (MQ=255)                                                                                                                      tCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATccc                 >  2:106015/1‑140 (MQ=255)                                                                                                                       cccTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCCTcccg                >  1:65042/1‑139 (MQ=255)                                                                                                                       cccTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCa                >  1:742456/1‑140 (MQ=255)                                                                                                                           ttACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGttt            >  1:263903/1‑140 (MQ=255)                                                                                                                           ttACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGttt            >  2:343587/1‑140 (MQ=255)                                                                                                                            tACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGtttt           >  2:1647032/1‑140 (MQ=255)                                                                                                                             aCCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGttttt          <  2:1215284/140‑1 (MQ=255)                                                                                                                              ccTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGtttttt         >  1:1300212/1‑140 (MQ=255)                                                                                                                               cTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTa        >  1:1500551/1‑140 (MQ=255)                                                                                                                                 gAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAgg      <  1:1089009/140‑1 (MQ=255)                                                                                                                                  aaGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGt     <  1:194331/140‑1 (MQ=255)                                                                                                                                  aaGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGt     >  2:1282953/1‑140 (MQ=255)                                                                                                                                  aaGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGt     >  1:675826/1‑140 (MQ=255)                                                                                                                                     ggggCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAt  <  2:1500551/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         TGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAT  >  NC_000913/4638608‑4638879

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