Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,161 G→A 6.4% noncoding (76/171 nt)
Q6* (CAA→TAA) 
rttR ←
tpr ←
rtT sRNA, processed from tyrT transcript
protamine‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,1610GA 6.4% 191.9 / 2.7 73noncoding (76/171 nt)
Q6* (CAA→TAA) 
rttR
tpr
rtT sRNA, processed from tyrT transcript
protamine‑like protein
Reads supporting (aligned to +/- strand):  ref base G (36/27);  new base A (2/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGG  >  NC_000913/1287012‑1287310
                                                                                                                                                     |                                                                                                                                                     
ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                                                                                                                                       <  2:1788919/150‑1 (MQ=255)
 gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                                                                      <  1:1604634/150‑1 (MQ=255)
       tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCACGCACGCGGGCGCTCGCGTAGAACCGTGGTCGAAg                                                                                                                                                >  1:1188411/1‑150 (MQ=14)
       tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACTCTCGCGGGCGCGCGAGTCGAACCTTGGTCGaac                                                                                                                                                >  1:1188577/1‑149 (MQ=14)
       tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACTCACGCTGGCGCTCGAGTCGAACCATGGTCGaaa                                                                                                                                                >  1:1188574/1‑149 (MQ=14)
         tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAAAGCTCACACGCTTTCGCTCGATTAGAACCATGGTCGAAGCt                                                                                                                                              >  1:1033545/1‑150 (MQ=14)
         tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGCGCTGCGCGCGTCGCCCGACGCGTCGATGCCACCGGCCACGCTCACTCGCTCGCGCTCGAGGCGAACCGTGGACGAAGca                                                                                                                                              >  2:1882822/1‑149 (MQ=14)
          ttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGGCCACGCTCTATAGCGGGCGCTCGACTCGAACCATTCTCGaatcat                                                                                                                                             >  1:2547256/1‑146 (MQ=14)
          ttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCAATCGCTGGCGCGCGAGTCGAACCTTGGTCGAAGcat                                                                                                                                             >  1:859783/1‑148 (MQ=255)
           tGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                                                            <  1:1745524/150‑1 (MQ=255)
           tGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCGCTCTCGCTGGCGCGCGAGTCGAAACATGGTCGAAGCGtc                                                                                                                                            >  2:448224/1‑150 (MQ=14)
           tGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGACCTTTCCTGCGGCAACGCTCTCTAGCTGGCGCGCGAGTCGAACCATGTTCGAAGCGtc                                                                                                                                            >  2:1072070/1‑150 (MQ=14)
                 cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                                                      <  1:2258431/150‑1 (MQ=255)
                          ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCGCGCATCGCCCTTCGGGTCGTTGCCTGAGGCAACGCTCTCTCGCTGTCGCTCGAATCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCa                                                                                                                             >  2:2453245/1‑150 (MQ=14)
                            gagaGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCACGCTGGCGCGCGAGTCGAACCATGGTAGAAGCATCTCAACCATACCCCCaag                                                                                                                           >  1:2387331/1‑148 (MQ=14)
                                 ttCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGGCGCTCGCCTCGAAACATGGTAGAAGCTTAACAAACTACCCCCCAAGGGCAg                                                                                                                      >  1:265017/1‑150 (MQ=14)
                                        gcTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                                             <  1:167587/135‑1 (MQ=37)
                                        gcTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                                             >  2:167587/2‑136 (MQ=37)
                                          ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAAAGCACTCAACCTGGCGCACGAGTCGAACCGTGGTCCAAACATCAAATCCATCCCCCCATTGGCAGAAAATTTaa                                                                                                             >  1:2358284/1‑150 (MQ=11)
                                          gtttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                                                             <  1:1147929/149‑1 (MQ=38)
                                           tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCACACGCTGGCGCTCGAGTCGAACCTTGGTCGAAACATATCAAACATACCCCCATGGGCAGAATATTTaaa                                                                                                            >  1:1692112/1‑149 (MQ=11)
                                              gacaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                                                         <  1:816168/147‑1 (MQ=255)
                                                           gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAATTAATTGCGGATTCGTTGGGa                                                                                            >  2:2538457/1‑150 (MQ=255)
                                                             gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAg                                                                                          <  1:1311259/150‑1 (MQ=255)
                                                                 gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCATCGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                                                                      <  1:1397502/150‑1 (MQ=255)
                                                                 gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                                                         >  2:326332/1‑99 (MQ=33)
                                                                 gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                                                         >  2:326345/1‑99 (MQ=33)
                                                                 gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                                                         <  1:326332/99‑1 (MQ=33)
                                                                 gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                                                         <  1:326345/99‑1 (MQ=33)
                                                                   aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGTAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                                                                    <  1:359641/150‑1 (MQ=255)
                                                                   aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                                                                    <  1:2459511/150‑1 (MQ=255)
                                                                     gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                                                      <  1:2487599/98‑1 (MQ=18)
                                                                      aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGAc                                                                                 >  2:1572221/1‑150 (MQ=255)
                                                                        tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtt                                                                               <  1:2448723/150‑1 (MQ=255)
                                                                        tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtt                                                                               <  2:880745/150‑1 (MQ=255)
                                                                           tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCGCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAGTCGTTGGGGAGTTCAGGGACTTTTg                                                                            >  2:121210/1‑150 (MQ=11)
                                                                                  tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                     <  1:584514/150‑1 (MQ=255)
                                                                                    cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCAGGGGCAGAAAATTGAAAATCGGAAACGTTGTGAAGTTAAGGGACGTGTGAGAGTGAtg                                                                   >  2:1077411/1‑150 (MQ=14)
                                                                                           cggccgtCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                                            <  1:1847229/144‑1 (MQ=255)
                                                                                           cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                                            <  1:1848165/150‑1 (MQ=255)
                                                                                             cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                                          >  2:870827/1‑150 (MQ=255)
                                                                                               cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACATTTGAAAGTGATGGTGGTGGGGGa                                                        >  2:1034266/1‑150 (MQ=255)
                                                                                                         ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa                                              >  2:1520989/1‑150 (MQ=255)
                                                                                                         ttGCCCGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCCTCTCTTCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa                                              <  2:689805/150‑1 (MQ=11)
                                                                                                         gtGCCGGCGGCAACTCTCTCTCGCTGCCGCTCGAGCCGAACCCTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa                                              <  2:734995/149‑1 (MQ=14)
                                                                                                            ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCt                                           >  2:1415448/1‑150 (MQ=255)
                                                                                                            ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGGGGTGGTGGTGGGGGAAGGATTCGAACCt                                           >  2:478323/1‑150 (MQ=255)
                                                                                                            ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATAGCGGATGCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCt                                           >  2:1415425/1‑150 (MQ=255)
                                                                                                             cTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtttt                                                                             <  1:1175578/115‑1 (MQ=255)
                                                                                                             cTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtttt                                                                             >  2:1175578/1‑115 (MQ=255)
                                                                                                               gCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCg                                        <  2:1188574/150‑1 (MQ=255)
                                                                                                               gCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCg                                        <  2:1188577/150‑1 (MQ=255)
                                                                                                               gCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCg                                        <  2:1188411/150‑1 (MQ=255)
                                                                                                                 ggCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGaa                                      <  2:1078269/150‑1 (MQ=255)
                                                                                                                    aaCGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTc                                   <  1:2543569/150‑1 (MQ=255)
                                                                                                                        ctctctCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTATCAACCTTCCCCCCAAGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATTGTGGTGGGGGAAGGATTCGAAACATCGAAGTAGaag                               >  1:2302925/1‑148 (MQ=17)
                                                                                                                         tctctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGa                              >  2:484455/1‑150 (MQ=255)
                                                                                                                         tctctcGCTGCCGCTCGAGACGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGa                              <  1:852734/150‑1 (MQ=255)
                                                                                                                           tctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                                                             <  2:1583391/69‑1 (MQ=18)
                                                                                                                             tcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCATCGAAGTCGATTACGGc                          >  1:1453212/1‑150 (MQ=255)
                                                                                                                              cGCTGCCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCa                         <  2:35628/150‑1 (MQ=255)
                                                                                                                              cGCCGGAGCTCGAGTCGATCCTTGGTCGAATGTCCTCAGCCTTCGCAGCACGGGTATAATATTTTATTGCGGATTCGTTTGGTATTTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCa                         <  2:35652/150‑1 (MQ=17)
                                                                                                                                     gcTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAc                  <  1:715544/150‑1 (MQ=255)
                                                                                                                                       tCGAGTCGAACCTTGGTCGAAGCTTCTCATACTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGGTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAg                >  1:139079/1‑150 (MQ=255)
                                                                                                                                         gAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTc              <  1:1164289/150‑1 (MQ=255)
                                                                                                                                            tCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGc           >  2:1713600/1‑150 (MQ=255)
                                                                                                                                              gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACttt                                                                              >  1:960404/1‑81 (MQ=255)
                                                                                                                                              gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACttt                                                                              <  2:960404/81‑1 (MQ=255)
                                                                                                                                               aaCCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                       >  2:2104057/1‑39 (MQ=40)
                                                                                                                                               aaCCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                       <  1:2104057/39‑1 (MQ=40)
                                                                                                                                               aaCCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                       <  1:1376375/39‑1 (MQ=40)
                                                                                                                                               aaCCTTGGTCGAAGCTTCTCATCCATACCCGCATGGGCAGAAAATTTAATTTCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATTGTGGTGGGGGAAGGATTCGAACCATCGAATTCGATGACGGCAGATTTAAAGTATTCTcc        >  1:1842318/1‑150 (MQ=18)
                                                                                                                                               aaCCTTGGTCGAAGCATCTCATACTTCCCCGCATGGGCAGAATAATTAATTGAGGATTCGTTGGGAAGTTCAGGGACTTTTTAAAGTGATTGTGGTGGGGGAAGGATTAGAACCATAGAAGTCGATGACGGCAGATTTACAGTCTGCAcc        >  1:1008425/1‑150 (MQ=21)
                                                                                                                                                   ttGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCttt    >  1:1885484/1‑150 (MQ=255)
                                                                                                                                                   ttGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGGTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCtgt    >  1:1885209/1‑148 (MQ=255)
                                                                                                                                                     ggTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTgg  <  2:662839/150‑1 (MQ=255)
                                                                                                                                                     ggTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTgg  <  1:1355743/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGG  >  NC_000913/1287012‑1287310

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