Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,402	+A	coding (596/1128 nt)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,402	1	.	A	97.9%	336.6 / ‑0.4	97	coding (596/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (53/42);  total (54/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTA  >  NC_000913/3379307‑3379494                                                                                                 |                                                                                             cATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGttataa                                                                                           >  2:149554/1‑96 (MQ=255)  aTCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGttataat                                                                                          >  2:788316/1‑95 (MQ=255)  aTCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGttataat                                                                                          >  1:1800077/1‑95 (MQ=255)     gTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGc                                                                                       >  1:2019016/1‑100 (MQ=255)      tGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCa                                                                                      <  1:2178407/100‑1 (MQ=255)       gAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCAt                                                                                     <  1:247662/100‑1 (MQ=255)        aaGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATc                                                                                    <  2:355141/100‑1 (MQ=255)          gTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCTCAATGCTGTTATAATGGCATCGa                                                                                  >  2:796631/1‑100 (MQ=255)          gTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGa                                                                                  >  2:952475/1‑100 (MQ=255)          gTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGa                                                                                  <  1:1182465/100‑1 (MQ=255)            ggCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGAtt                                                                                >  1:1207833/1‑100 (MQ=255)              cGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgc                                                                              >  2:2427259/1‑100 (MQ=255)              cGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgc                                                                              <  1:1591992/100‑1 (MQ=255)               gAAAGCCACAGATGCGCCTGAGTGAGTGTACGCCGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgca                                                                             >  1:2268160/1‑100 (MQ=255)               gAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgca                                                                             <  2:497728/100‑1 (MQ=255)               gAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgca                                                                             >  1:26151/1‑100 (MQ=255)               gAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgca                                                                             <  1:1102078/100‑1 (MQ=255)                aaaGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgcag                                                                            <  2:340497/100‑1 (MQ=255)                aaaGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgcag                                                                            <  1:2167895/100‑1 (MQ=255)                   gCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgcaggca                                                                         >  1:2044027/1‑100 (MQ=255)                   gCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgcaggca                                                                         <  1:1040427/100‑1 (MQ=255)                     cacaGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGa                                                                       >  1:1541397/1‑100 (MQ=255)                       caGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGaaa                                                                     >  1:562236/1‑100 (MQ=255)                        aGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGaaaa                                                                    <  2:806263/100‑1 (MQ=255)                         gATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAAc                                                                   >  1:879472/1‑100 (MQ=255)                         gATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAAc                                                                   >  1:2258517/1‑100 (MQ=255)                            gcgcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAAcgcg                                                                >  1:661252/1‑100 (MQ=255)                            gcgcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAAcgcg                                                                >  1:1019811/1‑100 (MQ=255)                             ngcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACgcgc                                                               >  2:343958/2‑100 (MQ=255)                              gcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCa                                                              >  1:1400129/1‑100 (MQ=255)                              gcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCa                                                              <  1:952483/100‑1 (MQ=255)                               ccTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCAc                                                             >  2:785303/1‑100 (MQ=255)                               ccTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCAc                                                             >  1:163215/1‑100 (MQ=255)                                cTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACg                                                            <  1:802465/100‑1 (MQ=255)                                 tgagtgagTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGt                                                           >  2:140019/1‑100 (MQ=255)                                    gtgagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTGATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTg                                                         >  2:1163399/1‑99 (MQ=255)                                    gtgagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGc                                                        <  2:1694320/100‑1 (MQ=255)                                    gtgagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCa                                                       >  2:2621008/1‑100 (MQ=255)                                     tgagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCa                                                       <  1:510662/100‑1 (MQ=255)                                      gagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAg                                                      <  1:68675/100‑1 (MQ=255)                                       agtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAgg                                                     >  1:747595/1‑100 (MQ=255)                                       agtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAgg                                                     <  2:152348/100‑1 (MQ=255)                                        gtgtACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAggn                                                    <  1:343962/100‑2 (MQ=255)                                        gtgtACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGc                                                    >  1:365775/1‑100 (MQ=255)                                        gtgtACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGc                                                    <  2:2547656/100‑1 (MQ=255)                                          gtACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCa                                                  <  1:464221/100‑1 (MQ=255)                                           tACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCAt                                                 >  1:1575280/1‑100 (MQ=255)                                           tACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCAt                                                 <  1:2585122/100‑1 (MQ=255)                                           tACGCAGACGATAATCAACACCAACGTCCACGTTCATTACCTCACAATGCTGTTATAATGGCATCGAGTGCAGGCAGAACACGCGCACGTTGCAAGCCa                                                  >  2:1367826/1‑99 (MQ=255)                                             cGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCAttt                                               >  1:39770/1‑100 (MQ=255)                                             cGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCAttt                                               >  1:2680778/1‑100 (MQ=255)                                                   cGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATaa                                         <  2:58402/100‑1 (MQ=255)                                                   cGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATaa                                         >  1:1047523/1‑100 (MQ=255)                                                     aTAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAg                                       <  1:1361336/100‑1 (MQ=255)                                                      tAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCGATTTCGATAAAGt                                      <  2:1663941/100‑1 (MQ=255)                                                      tAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGt                                      <  2:251948/100‑1 (MQ=255)                                                      tAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGt                                      <  2:2018994/100‑1 (MQ=255)                                                       aaTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGtt                                     <  1:1312583/100‑1 (MQ=255)                                                       aaTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGtt                                     >  1:270985/1‑100 (MQ=255)                                                       aaTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGtt                                     >  1:1551486/1‑100 (MQ=255)                                                        aTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTc                                    >  1:1936670/1‑100 (MQ=255)                                                        aTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTc                                    >  2:2485748/1‑100 (MQ=255)                                                        aTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTc                                    >  2:73562/1‑100 (MQ=255)                                                         tCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCa                                   >  1:842834/1‑100 (MQ=255)                                                          cAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCAt                                  >  2:2585405/1‑100 (MQ=255)                                                           aaCACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATc                                 <  1:2090417/100‑1 (MQ=255)                                                           aaCACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATc                                 >  1:2185680/1‑100 (MQ=255)                                                           aaCACCAGCGTCAACGTTCATTACATCACACTGCTGTTATAATGGCCTCGCTCCCAGGCAGCAAACGCGCCTGCTGCCGGCCAATTCGATACAGCTCATc                                 >  2:1838390/1‑100 (MQ=255)                                                            acacCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATcc                                <  1:350432/100‑1 (MQ=255)                                                             cacCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCcg                               >  2:2093297/1‑100 (MQ=255)                                                               ccAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCcggc                             >  2:1015652/1‑100 (MQ=255)                                                               ccAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCcggc                             >  1:1272518/1‑100 (MQ=255)                                                               caaGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCcggc                             >  1:1809212/3‑100 (MQ=255)                                                                    gTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGa                          >  1:548137/1‑98 (MQ=255)                                                                     tCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAata                       <  2:1311206/100‑1 (MQ=255)                                                                       cACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATAtt                     >  1:159584/1‑100 (MQ=255)                                                                        aCGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATAttt                    <  2:1805600/100‑1 (MQ=255)                                                                         cGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTg                   <  1:472850/100‑1 (MQ=255)                                                                          gTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGa                  >  1:33519/1‑100 (MQ=255)                                                                           ttCATTACATCACAATTCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGAc                 <  2:1575266/100‑1 (MQ=255)                                                                           ttCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGAc                 >  1:2287103/1‑100 (MQ=255)                                                                           ttCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGAc                 <  2:2044005/100‑1 (MQ=255)                                                                             cATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGt               <  1:302433/100‑1 (MQ=255)                                                                              aTTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTc              >  1:2347078/1‑100 (MQ=255)                                                                              aTTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTc              >  2:526613/1‑100 (MQ=255)                                                                              aTTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCg             <  1:2621036/100‑1 (MQ=255)                                                                               tcaCATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCAATAAAGTTCATCCGGCGGAATATTTGACGTCg             >  1:1825784/3‑100 (MQ=255)                                                                                tACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGc            <  1:2239747/100‑1 (MQ=255)                                                                                 aCATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCt           >  1:1079686/1‑100 (MQ=255)                                                                                  cATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTa          >  1:2028421/1‑100 (MQ=255)                                                                                    tCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTAcc        >  2:2660727/1‑100 (MQ=255)                                                                                    tCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTAcc        >  2:1034858/1‑100 (MQ=255)                                                                                        aaTGGTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAggg    <  2:1541383/100‑1 (MQ=255)                                                                                        aaTGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAggg    >  2:1965667/1‑100 (MQ=255)                                                                                        aaTGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAggg    <  1:2050242/100‑1 (MQ=255)                                                                                         aTGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGATGTCGCTACCAGGGt   <  2:1215359/100‑1 (MQ=255)                                                                                         aTGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGt   <  2:432591/100‑1 (MQ=255)                                                                                         aTGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGt   <  2:1331336/100‑1 (MQ=255)                                                                                          tGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTa  <  2:2185658/100‑1 (MQ=255)                                                                                          tGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTa  >  1:2488937/1‑100 (MQ=255)                                                                                                 |                                                                                             CATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTA  >  NC_000913/3379307‑3379494

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