Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,738,554 G→A 6.2% P237P (CCG→CCA)  purR → transcriptional repressor, hypoxanthine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,738,5540GA 6.2% 202.1 / 2.5 64P237P (CCG→CCA) purRtranscriptional repressor, hypoxanthine‑binding
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base A (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

TGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGC  >  NC_000913/1738406‑1738702
                                                                                                                                                    |                                                                                                                                                    
tGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCg                                                                                                                                                      >  1:19978/1‑149 (MQ=255)
 ggAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCg                                                                                                                                                      >  8:51084/1‑148 (MQ=255)
    aCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATc                                                                                                                                                  <  4:370872/149‑1 (MQ=255)
      gTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGCTGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGc                                                                                                                                                >  3:335063/1‑149 (MQ=255)
          cacCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTa                                                                                                                                            <  6:332814/149‑1 (MQ=255)
                     gccgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttc                                                                                                                                 >  5:321070/1‑149 (MQ=255)
                     gccgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttc                                                                                                                                 >  4:95081/1‑149 (MQ=255)
                     gccgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttc                                                                                                                                 >  2:333271/1‑149 (MQ=255)
                       cgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGACTCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtg                                                                                                                               >  2:251800/1‑149 (MQ=255)
                             gCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggtg                                                                                                                          <  5:79105/148‑1 (MQ=255)
                             gCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGtggtgg                                                                                                                         >  5:93273/1‑149 (MQ=255)
                              cgggTTTTATGTAGGCGTTGGAAGAAGAGATGATCATGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTTTACCTGAATGCGGTTATCGCGCCATGGAGCAATTCCTGTCGCAGCCGCTTGGCCCTACTGCAGTCTTCTGTGGTGGc                                                                                                                        <  5:192587/147‑1 (MQ=255)
                              ccGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGtggtgg                                                                                                                         <  6:321070/148‑1 (MQ=255)
                                 gtgttATGAAGGCGATGGAAGAAGCGATGAGCACGGTGCCGGAAAGCTGGAGTGTGCAGGGGGACTTTGAACCTGAATCCGGTTCTCGCGCCATGCAGCAAATCCTGGCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGat                                                                                                                     <  6:371603/146‑1 (MQ=255)
                                       tGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATg                                                                                                               <  1:141767/149‑1 (MQ=255)
                                         aaGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATgg                                                                                                              <  8:227869/148‑1 (MQ=255)
                                             cGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTGATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACAGCCGTCTTCTGTGGTGGCGATATCATGGCAATg                                                                                                         <  2:127634/149‑1 (MQ=255)
                                              gATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtgg                                                                                                                            >  4:267502/1‑129 (MQ=255)
                                              gATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtgg                                                                                                                            <  3:267502/129‑1 (MQ=255)
                                                   aagaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCa                                                                                                   >  5:67057/1‑149 (MQ=255)
                                                   aagaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCa                                                                                                   >  3:119047/1‑149 (MQ=255)
                                                       agCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGcc                                                                                                                                       >  8:103543/1‑109 (MQ=255)
                                                       agCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGcc                                                                                                                                       <  7:103543/109‑1 (MQ=255)
                                                         cGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTtgt                                                                                             >  1:357943/1‑149 (MQ=255)
                                                           atgatCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGtgc                                                                                           >  1:49291/1‑149 (MQ=255)
                                                            tgatCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGtgct                                                                                          >  7:98971/1‑149 (MQ=255)
                                                             gatCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGtgctg                                                                                         >  8:293127/1‑149 (MQ=255)
                                                              atCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGtgctgc                                                                                        <  7:51084/149‑1 (MQ=255)
                                                              atCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGtgctg                                                                                         <  3:319027/148‑1 (MQ=255)
                                                                       ccGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAAt                                                                               <  7:53420/149‑1 (MQ=255)
                                                                        cGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATg                                                                              <  6:93273/149‑1 (MQ=255)
                                                                        cGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATg                                                                              >  5:40879/1‑149 (MQ=255)
                                                                             aGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGAGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCt                                                                         <  2:161994/149‑1 (MQ=255)
                                                                               cTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGtgctgc                                                                                        <  7:126856/132‑1 (MQ=255)
                                                                               cTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGtgctgc                                                                                        >  8:126856/1‑132 (MQ=255)
                                                                               cTGGATTGTGCAGGGGGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTgc                                                                       <  1:348911/149‑1 (MQ=255)
                                                                                       tgCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCg                                                               >  7:105909/1‑149 (MQ=255)
                                                                                                  tttGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTTGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTc                                                    >  3:211450/1‑149 (MQ=255)
                                                                                                  tttGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTc                                                    >  2:18928/1‑149 (MQ=255)
                                                                                                    tGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCg                                                   <  6:67057/148‑1 (MQ=255)
                                                                                                      aaCCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCt                                                 <  1:215965/148‑1 (MQ=255)
                                                                                                         cTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGa                                                                                                                      >  6:166659/1‑76 (MQ=255)
                                                                                                         cTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGa                                                                                                                      <  5:166659/76‑1 (MQ=255)
                                                                                                         cTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGTCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATc                                             <  2:142883/149‑1 (MQ=255)
                                                                                                              tCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTa                                        <  7:293127/149‑1 (MQ=255)
                                                                                                              tCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTa                                        <  4:335063/149‑1 (MQ=255)
                                                                                                               ccGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTAt                                       >  5:202301/1‑149 (MQ=255)
                                                                                                                     aTCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATaa                                  <  2:357943/148‑1 (MQ=255)
                                                                                                                     aTCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAAc                                 >  5:120574/1‑149 (MQ=255)
                                                                                                                      tCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACg                                >  5:88562/1‑149 (MQ=255)
                                                                                                                       cgcgCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGt                               >  7:323344/1‑149 (MQ=255)
                                                                                                                        gcgcCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTg                              <  5:306107/149‑1 (MQ=255)
                                                                                                                        gcgcCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTg                              >  3:153336/1‑149 (MQ=255)
                                                                                                                          gcCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTTGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCCGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTgcg                            <  6:310997/149‑1 (MQ=255)
                                                                                                                              tGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAc                        <  2:19978/149‑1 (MQ=255)
                                                                                                                               gcagcaAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGACTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAcg                       >  4:144369/1‑149 (MQ=255)
                                                                                                                                cagcaAATCCTGTCGCAGCCACATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAcgc                      >  6:329320/1‑149 (MQ=255)
                                                                                                                                cagcaAATCCTGTCGCAGCCACATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAcgc                      >  6:329317/1‑149 (MQ=255)
                                                                                                                                 agcaAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAcgcg                     >  3:334501/1‑149 (MQ=255)
                                                                                                                                    aaaTCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTa                                        >  6:13529/1‑127 (MQ=255)
                                                                                                                                    aaaTCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTa                                        <  5:13529/127‑1 (MQ=255)
                                                                                                                                           gTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGGTATGATAACGTGCGCATCGCGCGCTATTTTa           >  1:321454/1‑149 (MQ=255)
                                                                                                                                            tCGCAGCCACATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTa           <  5:329320/148‑1 (MQ=255)
                                                                                                                                            tCGCAGCCACATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTa           <  5:329317/148‑1 (MQ=255)
                                                                                                                                                aGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCg      <  3:144369/149‑1 (MQ=255)
                                                                                                                                                   cGCATCGCCCTACTGCCCTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGgcg   >  7:337478/1‑149 (MQ=255)
                                                                                                                                                    gCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCGGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGgcgc  >  5:159884/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGC  >  NC_000913/1738406‑1738702

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