Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,734	Δ1 bp	noncoding (409/2906 nt)	rrlC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,941,734	0	A	.	98.8%	142.0 / ‑3.1	86	noncoding (409/2906 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (46/39);  minor base G (0/1);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAG  >  NC_000913/3941644‑3941784                                                                                           |                                                   gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                  <  1:596130/93‑1 (MQ=18) gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                  >  1:5610215/1‑93 (MQ=25) gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                  >  1:1853332/1‑93 (MQ=18) gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                  <  1:1016485/93‑1 (MQ=18) gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                  >  1:3026416/1‑93 (MQ=18)  aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                 >  1:410227/1‑92 (MQ=21)  aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                 >  1:5614513/1‑92 (MQ=21)  aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                 >  2:2407835/1‑92 (MQ=21)   tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                                               >  2:5567089/1‑93 (MQ=16)    aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAcc                                              >  2:755550/1‑93 (MQ=12)    aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAcc                                              >  2:14588/1‑93 (MQ=9)     cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCa                                             <  2:1898548/93‑1 (MQ=9)      aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAt                                            >  1:5982236/1‑93 (MQ=16)      aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAt                                            <  2:1126955/93‑1 (MQ=9)      aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAt                                            >  2:45667/1‑93 (MQ=9)      aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAt                                            <  2:1787691/93‑1 (MQ=9)       ggggGAAAGCCCCCGACACAAAAATGGTCATTTTGTGAGCTCGATGAGTAGGGCGGGACACGTGGGATTCTGTCTGCATATGGGGGGGAcatc                                            >  2:1592314/1‑90 (MQ=2)       gggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtc                                           >  2:221987/1‑93 (MQ=16)       gggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtc                                           >  1:85404/1‑93 (MQ=16)       gggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtc                                           <  2:724900/93‑1 (MQ=9)       gggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtc                                           >  2:1743091/1‑93 (MQ=16)       gggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtc                                           >  1:4454761/1‑93 (MQ=9)        ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcc                                          <  2:1868404/93‑1 (MQ=9)         gTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcct                                         >  2:315490/1‑93 (MQ=9)          tGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctc                                        >  1:5173509/1‑93 (MQ=12)           gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctcc                                       <  1:2096950/93‑1 (MQ=9)           gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctcc                                       >  2:5471873/1‑93 (MQ=9)           gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctcc                                       >  2:6267074/1‑93 (MQ=16)           gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctcc                                       <  1:3111791/93‑1 (MQ=9)            aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCa                                      >  1:2940798/1‑93 (MQ=16)             cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                     <  1:5775636/93‑1 (MQ=9)              aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                    >  1:461408/1‑93 (MQ=9)              aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                    >  1:3731753/1‑93 (MQ=16)              aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                    <  2:5400673/93‑1 (MQ=9)              aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                    >  2:2298504/1‑93 (MQ=16)              aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                    <  1:2890374/93‑1 (MQ=9)              aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                    <  1:2890373/93‑1 (MQ=9)               gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                   >  2:1944289/1‑93 (MQ=9)               gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                   >  2:5628423/1‑93 (MQ=16)               gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                   >  2:2413111/1‑93 (MQ=16)               gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                   <  2:5698004/93‑1 (MQ=9)                ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGc                                  <  2:4312851/93‑1 (MQ=9)                 cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt                                 >  1:916875/1‑93 (MQ=9)                 cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt                                 >  2:954142/1‑93 (MQ=16)                 cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt                                 <  1:4528074/93‑1 (MQ=9)                  ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                >  2:5810043/1‑93 (MQ=12)                  ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                <  1:381055/93‑1 (MQ=9)                  ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                >  1:2794414/1‑93 (MQ=16)                   cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaa                               <  1:4285811/93‑1 (MQ=9)                   cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaa                               >  1:4285199/1‑93 (MQ=16)                   cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaa                               >  2:5428635/1‑93 (MQ=16)                    gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                              <  2:1543994/93‑1 (MQ=9)                    gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                              <  2:375580/93‑1 (MQ=9)                    gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                              >  1:1626718/1‑93 (MQ=9)                     tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                             <  1:3404247/93‑1 (MQ=9)                     tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                             >  1:4620191/1‑93 (MQ=12)                     tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                             <  1:5957467/93‑1 (MQ=9)                      acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                            <  2:3723540/93‑1 (MQ=9)                      acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                            <  1:2353402/93‑1 (MQ=9)                      acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                            >  1:1706133/1‑93 (MQ=16)                      acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGGATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                            >  1:4954932/1‑93 (MQ=9)                      acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGCCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAtc                            <  2:1179537/93‑2 (MQ=255)                       cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATAc                           >  1:6291200/1‑93 (MQ=16)                       cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATAc                           >  1:768874/1‑93 (MQ=12)                        acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACt                          <  1:5809219/93‑1 (MQ=9)                        acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGACTGGCGATGGG‑GGGACCATCCTACAAGGCTAAATACt                          >  1:557172/1‑93 (MQ=11)                         caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTc                         >  2:3352749/1‑93 (MQ=16)                         caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTc                         <  1:5347879/93‑1 (MQ=9)                         caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTc                         <  2:3915086/93‑1 (MQ=9)                          aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                        <  2:5358985/93‑1 (MQ=9)                           aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                       <  2:651682/93‑1 (MQ=9)                           aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                       <  2:5786396/93‑1 (MQ=9)                           aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                       >  1:2510881/1‑93 (MQ=9)                           aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                       <  1:2739264/93‑1 (MQ=9)                            aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctg                      <  2:76802/93‑1 (MQ=9)                            aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctg                      >  1:1699175/1‑93 (MQ=12)                             aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctga                     <  1:5237373/93‑1 (MQ=9)                             aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctga                     >  2:669008/1‑93 (MQ=9)                             aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctga                     >  2:1249515/1‑93 (MQ=16)                              aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgac                    >  2:4590979/1‑93 (MQ=16)                              aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgac                    <  2:1566712/93‑1 (MQ=9)                               tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                   >  2:2276523/1‑93 (MQ=16)                               tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                   <  1:37577/93‑1 (MQ=9)                               tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                   >  2:3624838/1‑93 (MQ=16)                               tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                   <  2:3717237/93‑1 (MQ=9)                                gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactg                  <  1:4723507/93‑1 (MQ=9)                                 cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactgc                 >  1:2040853/1‑92 (MQ=12)                                 cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactga                 >  1:518202/1‑93 (MQ=12)                                 cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactga                 <  2:632383/93‑1 (MQ=9)                                 cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGAATCCTGTCTGAATATGGG‑GGGACCATCGTCCAAGGCTAAATACTCctcactga                 >  2:5787065/1‑93 (MQ=11)                                  acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATAGTCCtgactgac                <  2:3906640/93‑1 (MQ=255)                                  acaTATTGTGAGCTCGATGAGTACGGCGGGACACGTGGGATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTTCtgactgac                >  1:3187509/1‑93 (MQ=11)                                   caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc               >  1:4863460/1‑93 (MQ=0)                                   caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc               >  2:1797931/1‑93 (MQ=0)                                   caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc               >  1:4736235/1‑93 (MQ=0)                                   caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc               <  1:3152013/93‑1 (MQ=0)                                                 gATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAGCCAg  <  2:5607765/93‑1 (MQ=16)                                                                                           |                                                   GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAG  >  NC_000913/3941644‑3941784

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