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	96.9%	320.3 / 0.3	98	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/2);  new base A (55/39);  total (57/41)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATT  >  NC_000913/1728348‑1728540                                                                                                 |                                                                                                 tGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGa                                                                                               >  2:1453182/1‑100 (MQ=255)     cGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATgcg                                                                                           >  1:1928242/1‑100 (MQ=255)     cGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATgcg                                                                                           >  2:1881156/1‑100 (MQ=255)       aTAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCgct                                                                                         >  1:88690/1‑100 (MQ=255)            gCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACCGCCGGATTTACCGCCAAAAACGATGCGCTGCTTg                                                                                    >  2:1226127/1‑100 (MQ=255)            gCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTg                                                                                    <  1:978625/100‑1 (MQ=255)            gCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTg                                                                                    >  2:2018918/1‑100 (MQ=255)            gCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTg                                                                                    <  1:1575007/100‑1 (MQ=255)             cTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTga                                                                                   <  2:1763717/100‑1 (MQ=255)                    tACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgcc                                                                            <  1:1982058/100‑1 (MQ=255)                     aCCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccg                                                                           <  1:1358910/100‑1 (MQ=255)                     aCCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccg                                                                           <  1:2397522/100‑1 (MQ=255)                      ccGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccgc                                                                          >  1:9493/1‑100 (MQ=255)                        ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                        <  1:1183219/100‑1 (MQ=255)                        ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                        <  1:588864/100‑1 (MQ=255)                        ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATATAAAGAAAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                        <  2:232128/100‑1 (MQ=255)                            tgtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAc                                                                    <  2:1952013/100‑1 (MQ=255)                            tgtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAc                                                                    >  2:87007/1‑100 (MQ=255)                             gtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAcc                                                                   >  1:507103/1‑100 (MQ=255)                             gtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAcc                                                                   >  1:1900127/1‑100 (MQ=255)                              tgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAccc                                                                  <  1:1119283/100‑1 (MQ=255)                              tgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAccc                                                                  >  1:202048/1‑100 (MQ=255)                               gCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCt                                                                 <  1:1626444/100‑1 (MQ=255)                                gatctgTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTg                                                                <  2:232602/95‑1 (MQ=255)                                cGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTg                                                                >  1:1416159/1‑100 (MQ=255)                                cGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTg                                                                >  1:489768/1‑100 (MQ=255)                                cGACCGATTTCGTGGGTTTTATCCTGTTGTGATCGCTGTTGCATCAGCCGGTTTTTACGCCAAAACCGATGCGCTGATTGAGATTGCCGCCATCACCCt                                                                 >  2:2386820/1‑99 (MQ=255)                                 gACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGa                                                               >  1:805494/1‑100 (MQ=255)                                  aCCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGaa                                                              <  1:1453193/100‑1 (MQ=255)                                     gTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGagatc                                                           >  1:232604/1‑96 (MQ=255)                                      ttttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATg                                                          <  2:255950/100‑1 (MQ=255)                                      ttttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATg                                                          <  1:2415095/100‑1 (MQ=255)                                       tttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATgg                                                         >  1:299419/1‑100 (MQ=255)                                       nttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATgg                                                         >  2:1420476/2‑100 (MQ=255)                                        ttCGTGGTTTTTATCTTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGa                                                        <  1:713210/100‑1 (MQ=255)                                        ttCGTCGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGa                                                        <  2:55461/100‑1 (MQ=255)                                         tCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGAt                                                       <  2:1570003/100‑1 (MQ=255)                                          cGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATg                                                      <  1:2018932/100‑1 (MQ=255)                                          cGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATg                                                      <  1:1259070/100‑1 (MQ=255)                                               tttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACaa                                                 >  2:1814124/1‑100 (MQ=255)                                                ttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGCCATTGCCGCCATCACCCTGAAAATGGATGAACAAg                                                >  1:650340/1‑100 (MQ=255)                                                ttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAg                                                >  1:977533/1‑100 (MQ=255)                                                ttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAg                                                >  2:625067/1‑100 (MQ=255)                                                  ttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggc                                              <  1:2338362/100‑1 (MQ=255)                                                   tATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggct                                             <  1:391700/100‑1 (MQ=255)                                                    aTCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctg                                            <  1:863448/100‑1 (MQ=255)                                                      ccTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggc                                          >  1:174064/1‑100 (MQ=255)                                                       cTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggct                                         >  2:2110939/1‑100 (MQ=255)                                                         gttgtGATCGATGTTGAAACCGCCGGCTTTAACGCCAAAAACGATGGGCTGCTTGAGATTGCCGGCATCACCCTGAAAACGGATGAACAGgactgact                                         >  2:2202694/1‑98 (MQ=255)                                                          ttgtGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGAt                                      >  1:154449/1‑100 (MQ=255)                                                            gtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGc                                    >  2:1056899/1‑100 (MQ=255)                                                            gtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGACGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGc                                    >  2:1059364/1‑100 (MQ=255)                                                              gATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCg                                  >  2:592059/1‑100 (MQ=255)                                                              gATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAAGGGATGAACAAGGCTGGCTGATGCCg                                  >  2:2404565/1‑100 (MQ=255)                                                               attgATGTTGAAACAGCCGGATGTAACGCCAAAAACGATGCGGTGCTTGAGATTGCGGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                 <  2:1803410/97‑1 (MQ=255)                                                               aTCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                 <  2:1620606/100‑1 (MQ=255)                                                               aTCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                 <  2:2103631/100‑1 (MQ=255)                                                                tCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGa                                >  1:1027415/1‑100 (MQ=255)                                                                  gATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGaca                              <  1:2276505/100‑1 (MQ=255)                                                                  gATGTTGAAACAGCCGGATTCAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGaca                              <  1:389583/100‑1 (MQ=255)                                                                    tGTTGAAACAGCCGGATTTAGCGCCAAAAACGATCCGCTGCTTGAGATTGCGGCCATCGCGCGGAAGATGGATGAACAAGGCTGGCTGATGCCgg                                 <  2:1162423/95‑1 (MQ=255)                                                                      ttGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTCATGCCGGAcaccac                          >  2:2172996/1‑100 (MQ=255)                                                                       tGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACcaca                         >  2:727435/1‑100 (MQ=255)                                                                        gAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacat                        >  1:1049416/1‑100 (MQ=255)                                                                         aaaCAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacatt                       >  2:230390/1‑100 (MQ=255)                                                                           aCAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattac                     <  2:272550/100‑1 (MQ=255)                                                                           aCAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattac                     <  2:2499704/100‑1 (MQ=255)                                                                              gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                  >  1:1168854/1‑100 (MQ=255)                                                                              gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                  >  1:1257376/1‑100 (MQ=255)                                                                              gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                  >  2:1816865/1‑100 (MQ=255)                                                                              gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                  <  2:2540382/100‑1 (MQ=255)                                                                              gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                  >  2:2473751/1‑100 (MQ=255)                                                                              gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattac                     <  2:832353/97‑1 (MQ=255)                                                                              gCCGGATCTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                  >  1:1343521/1‑100 (MQ=255)                                                                               ccGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                  >  2:182028/1‑99 (MQ=255)                                                                               ccGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACAttt                 >  2:336438/1‑100 (MQ=255)                                                                               ccGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACAttt                 >  1:2541282/1‑100 (MQ=255)                                                                                cGGATTTAACGCCAAAAACGGTGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTc                <  2:1748456/100‑1 (MQ=255)                                                                                 ggATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTcc               <  1:1885262/100‑1 (MQ=255)                                                                                 ggATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTcc               <  1:947637/100‑1 (MQ=255)                                                                                  gATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCa              <  2:2467160/100‑1 (MQ=255)                                                                                   aTTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGTCTGGCTGATGCCGGACACCACATTACATTTCCAc             >  2:1351907/1‑100 (MQ=255)                                                                                    tttAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACg            >  1:1711761/1‑100 (MQ=255)                                                                                    tttAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACg            >  1:2259476/1‑100 (MQ=255)                                                                                    tttAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACg            <  1:1088360/100‑1 (MQ=255)                                                                                    tttAACGCCAAAAACGATGCGCTGCTTGAGATTGCCCCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCa              <  2:2560944/98‑1 (MQ=255)                                                                                     ttAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGt           >  1:139305/1‑100 (MQ=255)                                                                                      tAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc          >  2:230193/1‑100 (MQ=255)                                                                                      tAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc          >  2:1092899/1‑100 (MQ=255)                                                                                      tAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc          <  1:143692/100‑1 (MQ=255)                                                                                      tAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc          >  1:2302148/1‑100 (MQ=255)                                                                                        aCGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGa        >  1:1873862/1‑100 (MQ=255)                                                                                          gCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAc      >  2:2205616/1‑100 (MQ=255)                                                                                           ccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAcc     >  1:721069/1‑100 (MQ=255)                                                                                           ccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAcc     >  1:2148918/1‑100 (MQ=255)                                                                                           ccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAcc     >  1:2112607/1‑100 (MQ=255)                                                                                           ccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAcc     >  1:2112435/1‑100 (MQ=255)                                                                                            cAAAATCGATGCGATGCTTGAGATTGCCGCCATCACACTGAAAATGGATGCACAAGTCTGGCTAATACCGAACACTACATTACATATCAACCTCGAAc      >  2:1495397/1‑98 (MQ=255)                                                                                              aaaaCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCAtt  <  2:507097/100‑1 (MQ=255)                                                                                              aaaaCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCAtt  <  2:9493/100‑1 (MQ=255)                                                                                                 |                                                                                                 TGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATT  >  NC_000913/1728348‑1728540

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