Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,728,444	C→A	100%	T33N (ACC→AAC)	rnt →	RNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,728,444	0	C	A	98.9%	295.1 / ‑5.9	89	T33N (ACC→AAC)	rnt	RNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (60/27);  total (61/28)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCA  >  NC_000913/1728348‑1728538                                                                                                  |                                                                                               tGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACACA‑AAAAACGa                                                                                             >  1:1446040/1‑100 (MQ=255) tGTCCGACAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGa                                                                                             >  2:1257102/1‑100 (MQ=255)  gTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGAt                                                                                            <  1:769591/100‑1 (MQ=255)  gTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGAt                                                                                            >  2:2192152/1‑100 (MQ=255)   tCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATg                                                                                           <  1:2568006/100‑1 (MQ=255)    ccGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATgc                                                                                          >  2:2148075/1‑100 (MQ=255)      gATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATgcgc                                                                                        >  2:1361982/1‑100 (MQ=255)      gATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATgcgc                                                                                        >  2:239683/1‑100 (MQ=255)          aCGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTATGAAAGATGTTGAAACAACACGATTTAACGCC‑AAAACCGATGCgctgct                                                                                    >  1:248121/1‑100 (MQ=255)           cGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGTTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCtt                                                                                   <  1:2758117/100‑1 (MQ=255)           cGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAACACGATGCGCTGCtt                                                                                   >  2:2789957/1‑100 (MQ=255)           cGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCtt                                                                                   >  1:1779862/1‑100 (MQ=255)            gCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGCTGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCtt                                                                                   >  2:2093316/1‑99 (MQ=255)            gCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTg                                                                                  >  1:870068/1‑100 (MQ=255)            gCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTg                                                                                  <  1:297231/100‑1 (MQ=255)               cAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTgaga                                                                               <  1:1624642/100‑1 (MQ=255)                 aCTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAAACGGATAAAACGCC‑AAAAACGATGCGCTGCTTGAGAtt                                                                             >  1:1580586/1‑100 (MQ=255)                      ccGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTgccgc                                                                        >  1:532630/1‑100 (MQ=255)                      ccGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTgccgc                                                                        >  1:2568893/1‑100 (MQ=255)                      ccGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTgccgc                                                                        <  1:282328/100‑1 (MQ=255)                         gTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCAt                                                                     >  1:1731504/1‑100 (MQ=255)                         gTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCAt                                                                     <  1:1391526/100‑1 (MQ=255)                          tCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATc                                                                    <  1:2947215/100‑1 (MQ=255)                           cTGTGCGACCGTTTTCGTGGTTTTTATCCTTTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCa                                                                   <  1:236258/100‑1 (MQ=255)                           cTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCa                                                                   >  2:447421/1‑100 (MQ=255)                            tgtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCAc                                                                  >  2:336897/1‑100 (MQ=255)                            tgtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCAc                                                                  >  1:1020801/1‑100 (MQ=255)                               gCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCt                                                               >  2:1609920/1‑100 (MQ=255)                               gCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCt                                                               <  2:1800543/100‑1 (MQ=255)                                   ccGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGaaa                                                           >  2:577675/1‑100 (MQ=255)                                   ccGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGaaa                                                           >  1:3030658/1‑100 (MQ=255)                                       tttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATgg                                                       <  2:2422390/100‑1 (MQ=255)                                        ttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGa                                                      >  1:2217563/1‑100 (MQ=255)                                        ttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGATGCTTGAGATTGCCGCCATCACCCTGAAAATGGa                                                      >  1:245170/1‑100 (MQ=255)                                         tCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGAt                                                     <  1:1257114/100‑1 (MQ=255)                                         tCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGAt                                                     >  2:1914378/1‑100 (MQ=255)                                          cGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATg                                                    <  1:457472/100‑1 (MQ=255)                                           gTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATg                                                    <  2:243691/100‑1 (MQ=255)                                             ggTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTACTTAAGATTGCCGCCATCACCCTGAAAATGGATGAAc                                                 >  1:809/1‑100 (MQ=255)                                              gTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACa                                                <  2:2404810/100‑1 (MQ=255)                                                ttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAg                                              <  2:1066420/100‑1 (MQ=255)                                                      ccTGTTGTGCTCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTACCGCCATCACCCTGAAAATGGATGAACAAggctggc                                        >  2:1552426/1‑100 (MQ=255)                                                      ccTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggc                                        <  2:1673579/100‑1 (MQ=255)                                                      ccTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCAGCCATCACCCTGAAAATGGATGAACAAggctggc                                        >  2:1171762/1‑100 (MQ=255)                                                       ntgttgtGATCGATGTTGAAGCAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggct                                       >  2:1848699/2‑100 (MQ=255)                                                       cTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggct                                       <  1:451696/100‑1 (MQ=255)                                                       cTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGAGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggct                                       >  2:1803996/1‑100 (MQ=255)                                                        tgttgtGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGgctggctg                                      >  1:544780/1‑100 (MQ=255)                                                        tgttgtGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGgctggctg                                      <  2:487358/100‑1 (MQ=255)                                                         gttgtGATCGATGTTGAAACAGCTGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCGATCACCCTGAAAATGGATGAACAAGgctggctg                                      <  2:2039931/99‑1 (MQ=255)                                                          ttgtGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGgctggctggt                                    >  2:927951/1‑98 (MQ=255)                                                           tgtgATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATg                                   >  2:2391456/1‑100 (MQ=255)                                                           tgtgATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATg                                   >  2:580678/1‑100 (MQ=255)                                                           tgtgATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATg                                   >  1:2996992/1‑100 (MQ=255)                                                            gtgATCGATGTTGAATCAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGc                                  <  1:2995072/100‑1 (MQ=255)                                                            gtgATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAATGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGc                                  >  1:2360913/1‑100 (MQ=255)                                                             tgATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGcc                                 >  1:2184580/1‑100 (MQ=255)                                                              gATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGACTGGCTGATGCCg                                >  2:2613903/1‑100 (MQ=255)                                                                 cGATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGa                              >  1:234357/1‑99 (MQ=255)                                                                  gATGTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGaca                            <  2:79660/100‑1 (MQ=255)                                                                    tGTTGAAACAGCCGGATTTAACGCC‑AAAAAAGATGCGCTGCTTGAGAATGACGCAAACACCCTGAAAACGGATGAAAAAGGCTGGCAGATGCACGaca                            >  2:882657/1‑98 (MQ=255)                                                                     gTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcacca                         <  1:447424/100‑1 (MQ=255)                                                                     gTTGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcacca                         >  1:1625741/1‑100 (MQ=255)                                                                      ttGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcaccac                        >  1:1296948/1‑100 (MQ=255)                                                                       tGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACcaca                       >  1:2703373/1‑100 (MQ=255)                                                                       tGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACcaca                       >  1:2703813/1‑100 (MQ=255)                                                                       tGAAACAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACcaca                       >  1:1090819/1‑100 (MQ=255)                                                                         aaaCAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacatt                     >  2:2639001/1‑100 (MQ=255)                                                                          aaCAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacatta                    >  2:64408/1‑100 (MQ=255)                                                                           aCAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattac                   <  1:1654832/100‑1 (MQ=255)                                                                            cAGCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattaca                  >  1:931406/1‑100 (MQ=255)                                                                              gCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                >  2:115747/1‑100 (MQ=255)                                                                              gCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                >  2:112400/1‑100 (MQ=255)                                                                              gCCGGATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                >  1:2963405/1‑100 (MQ=255)                                                                                 ggATTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTcc             >  2:2216009/1‑100 (MQ=255)                                                                                   aTTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGTTGGCTGATGCCGGACACCACATTACATTTCCan           <  2:221260/100‑2 (MQ=255)                                                                                   aTTTAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCAc           >  2:1335602/1‑100 (MQ=255)                                                                                    tttAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACg          >  1:1676277/1‑100 (MQ=255)                                                                                      tAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc        >  1:1946628/1‑100 (MQ=255)                                                                                      tAACGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCAGCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACAACACATTACATTTCCACGTc        >  2:1863462/1‑100 (MQ=255)                                                                                       aaCGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCg       >  1:1148045/1‑100 (MQ=255)                                                                                       aaCGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCg       <  2:1537355/100‑1 (MQ=255)                                                                                         cGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGaa     >  2:1552608/1‑100 (MQ=255)                                                                                         cGCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCAGGCTGATACCGGACACCACATTACATTTCCACGTCGaa     >  2:595789/1‑100 (MQ=255)                                                                                          gCC‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACATCGAAc    >  2:146224/1‑100 (MQ=255)                                                                                            c‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGGCGGACACCACATTACATTTCCACGTCGAACCa  <  1:115748/100‑1 (MQ=255)                                                                                            c‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCa  <  1:112401/100‑1 (MQ=255)                                                                                            c‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCa  >  2:2037138/1‑100 (MQ=255)                                                                                            c‑AAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCa  >  1:2345345/1‑100 (MQ=255)                                                                                                  |                                                                                               TGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCC‑AAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCA  >  NC_000913/1728348‑1728538

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