Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,194,516	A→G	K551K (AAA→AAG)	yqiK →	PHB family membrane protein, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,194,516	0	A	G	100.0%	365.7 / NA	98	K551K (AAA→AAG)	yqiK	PHB family membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (54/44);  total (54/44)

CTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGT  >  NC_000913/3194421‑3194612                                                                                                |                                                                                                 cTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCag                                                                                              <  2:2001354/100‑1 (MQ=255)   gATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCGTTAACCTCAACAACTCCCGTCGAAGAAAAGGCagag                                                                                            <  2:1569078/100‑1 (MQ=255)   gATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCagag                                                                                            <  2:418696/100‑1 (MQ=255)   gATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCagag                                                                                            >  2:1714047/1‑100 (MQ=255)    aTTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCagag                                                                                            >  1:1748634/1‑99 (MQ=255)    aTTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGt                                                                                           >  2:426137/1‑100 (MQ=255)    aTTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGt                                                                                           >  2:1216958/1‑100 (MQ=255)    aTTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGt                                                                                           >  2:1820249/1‑100 (MQ=255)     ttGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTa                                                                                          >  2:1233640/1‑100 (MQ=255)        aCTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAc                                                                                       <  2:49498/100‑1 (MQ=255)           ccTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                    >  2:758784/1‑100 (MQ=255)           ccTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                    >  1:1568849/1‑100 (MQ=255)            cTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCctc                                                                                   >  2:246691/1‑100 (MQ=255)            cTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTAGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                    <  2:1231586/99‑1 (MQ=255)             ttGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCctct                                                                                  <  2:1634118/100‑1 (MQ=255)               gCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTACCGTCGAAGAAAAGGCAGAGTAAACCCTCTTc                                                                                >  2:1744573/1‑100 (MQ=255)                cTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTcc                                                                               >  2:877087/1‑100 (MQ=255)                cTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGACAGAGTAAACCCTCTTcc                                                                               <  2:104188/100‑1 (MQ=255)                  cAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                             <  1:1021495/100‑1 (MQ=255)                  cAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                             >  2:1798675/1‑100 (MQ=255)                   aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                            <  2:846740/100‑1 (MQ=255)                   aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                            <  2:1309366/100‑1 (MQ=255)                    aTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGa                                                                           >  1:160990/1‑100 (MQ=255)                    aTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGa                                                                           >  1:2050377/1‑100 (MQ=255)                     tGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaa                                                                          >  2:1312902/1‑100 (MQ=255)                       aaaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaata                                                                        >  2:129733/1‑100 (MQ=255)                        aaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaataa                                                                       >  1:856962/1‑100 (MQ=255)                        aaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaataa                                                                       >  1:715925/1‑100 (MQ=255)                         aTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaataat                                                                      <  1:170142/100‑1 (MQ=255)                          ttGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATc                                                                     >  1:2162331/1‑100 (MQ=255)                          ttGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATc                                                                     <  1:240705/100‑1 (MQ=255)                           tGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATcc                                                                    <  1:271987/100‑1 (MQ=255)                           tGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATcc                                                                    <  1:699609/100‑1 (MQ=255)                            ggCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATccc                                                                   >  2:1213734/1‑100 (MQ=255)                             gCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCt                                                                  >  2:2035249/1‑100 (MQ=255)                              cGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTg                                                                 <  2:1469511/100‑1 (MQ=255)                                   cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCa                                                            >  1:224792/1‑100 (MQ=255)                                     ggcggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaa                                                          >  2:42436/1‑100 (MQ=255)                                     ggcggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaa                                                          <  1:1177158/100‑1 (MQ=255)                                       cggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaaaa                                                        >  2:1088272/1‑100 (MQ=255)                                        ggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAg                                                       >  1:1947578/1‑100 (MQ=255)                                          cTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGGAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGcc                                                     >  1:1048249/1‑100 (MQ=255)                                          cTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGcc                                                     >  1:1702965/1‑100 (MQ=255)                                           tCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCg                                                    <  2:1908036/100‑1 (MQ=255)                                            cACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCgg                                                   <  2:904842/100‑1 (MQ=255)                                            cACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCgg                                                   >  1:2049490/1‑100 (MQ=255)                                            cACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCg                                                    >  1:1441078/1‑99 (MQ=255)                                             aCTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCggg                                                  <  1:1308740/100‑1 (MQ=255)                                                ggcggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAt                                               >  1:1171498/1‑100 (MQ=255)                                                  cggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAttt                                             >  2:267738/1‑100 (MQ=255)                                                    gcATTGACTTCACCCTTAACCTCAACAGCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAttttt                                           <  1:577277/100‑1 (MQ=255)                                                    gcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAttttt                                           >  2:1600542/1‑100 (MQ=255)                                                     cATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                          <  1:775383/100‑1 (MQ=255)                                                     cATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                          <  1:1721870/100‑1 (MQ=255)                                                     cATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                          <  1:1404717/100‑1 (MQ=255)                                                      aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                         <  1:802782/100‑1 (MQ=255)                                                        tGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTtata                                       >  1:999685/1‑100 (MQ=255)                                                         gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTatat                                      >  2:2202089/1‑100 (MQ=255)                                                         gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTatat                                      >  2:1575972/1‑100 (MQ=255)                                                         gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTatat                                      >  1:1541650/1‑100 (MQ=255)                                                          aCTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATc                                     <  2:872943/100‑1 (MQ=255)                                                           cTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCt                                    >  1:660579/1‑100 (MQ=255)                                                            ttCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTg                                   <  2:2162343/100‑1 (MQ=255)                                                            ttCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTg                                   <  1:1820239/100‑1 (MQ=255)                                                            ttCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTg                                   >  1:801699/1‑100 (MQ=255)                                                              cACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGTCCAAAAAGCCGGGGATTTTTTATATCTGCg                                 <  1:2146174/100‑1 (MQ=255)                                                              cACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCg                                 <  1:136866/100‑1 (MQ=255)                                                                cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGtt                               >  1:2124566/1‑100 (MQ=255)                                                                cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGtt                               <  2:455087/100‑1 (MQ=255)                                                                  cTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTcc                             >  1:1520991/1‑100 (MQ=255)                                                                  cTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTcc                             <  2:1008903/100‑1 (MQ=255)                                                                  cTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTcc                             >  2:1345438/1‑100 (MQ=255)                                                                    tAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGc                           >  2:806001/1‑100 (MQ=255)                                                                     aaCCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCt                          <  2:337662/100‑1 (MQ=255)                                                                      aCCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTa                         >  2:253898/1‑100 (MQ=255)                                                                       ccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCGCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaa                        >  1:480723/1‑100 (MQ=255)                                                                        cTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaa                       >  1:1758436/1‑100 (MQ=255)                                                                         tCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCTCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaa                      <  1:1074638/100‑1 (MQ=255)                                                                         tCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaa                      <  2:2071157/100‑1 (MQ=255)                                                                         tCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaa                      >  2:981318/1‑100 (MQ=255)                                                                             ccacTCCCGATGAAGAAAAGGCAGGGTAAACCCTCTTCCAGCAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAg                    <  2:1384575/96‑1 (MQ=21)                                                                              aacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGc                 <  2:1208295/100‑1 (MQ=255)                                                                                 tCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCaaa              <  1:899528/100‑1 (MQ=255)                                                                                  cccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt             >  1:1445322/1‑100 (MQ=255)                                                                                  cccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt             >  1:1468368/1‑100 (MQ=255)                                                                                  cccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt             >  1:678006/1‑100 (MQ=255)                                                                                   ccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg            <  2:391715/100‑1 (MQ=255)                                                                                   ccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg            >  2:169273/1‑100 (MQ=255)                                                                                   ccGTCGAAGAAAAGGCAGAGAAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg            >  1:1359972/1‑100 (MQ=255)                                                                                     gTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGAt          <  2:241940/100‑1 (MQ=255)                                                                                      tCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATc         <  1:1088271/100‑1 (MQ=255)                                                                                       cGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCa        <  1:1831970/100‑1 (MQ=255)                                                                                       cGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCa        <  2:2049502/100‑1 (MQ=255)                                                                                        gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg       >  2:1259346/1‑100 (MQ=255)                                                                                        gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg       >  1:1148275/1‑100 (MQ=255)                                                                                         aagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAgg      <  2:660580/100‑1 (MQ=255)                                                                                         aagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAgg      <  1:877086/100‑1 (MQ=255)                                                                                          agaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgt     >  2:902050/1‑100 (MQ=255)                                                                                             aaagGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAACGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGcttg   <  1:1823560/99‑1 (MQ=255)                                                                                             aaagGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgt  >  2:1735840/1‑100 (MQ=255)                                                                                                |                                                                                                 CTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGT  >  NC_000913/3194421‑3194612

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