Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,258,185 C→T 6.0% G110D (GGC→GAC)  pth ← peptidyl‑tRNA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,258,1850CT 6.0% 220.0 / 2.6 67G110D (GGC→GAC) pthpeptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base T (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TAACAGGCGGTTTGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTGAGATTCATAAATGTAGTCGGGACTAA  >  NC_000913/1258045‑1258330
                                                                                                                                            |                                                                                                                                                 
tAAGAGGCGGTTTGACTAACACTAAAACGAGAACTTTATTTTTATCTCCCGGATGTCCGATTCCGATGCGTATACGGTGAAAGTTAGGGTTATTACCCAATTTACTTATGATGTCTTTCAGTCCATTGTGACCACCATGGCCAcc                                                                                                                                               <  6:282759/145‑1 (MQ=255)
tAACAGGCGGTTTGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGcc                                                                                                                                            >  2:159064/1‑148 (MQ=255)
tAACAGGCGGTTTGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGcc                                                                                                                                            >  4:76423/1‑148 (MQ=255)
 aaCAGGCGGTTTGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGcc                                                                                                                                                  >  1:265368/1‑141 (MQ=255)
 aaCAGGCGGTTTGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGcc                                                                                                                                                  <  2:265368/141‑1 (MQ=255)
 aaCAGGCGGTTTGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGcc                                                                                                                                            >  1:221331/1‑147 (MQ=255)
 aaCAGGCGGTATGCCTAACACAAAACAGACAACTTTATTTTTATCGACCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACAACCAAGGCCACCg                                                                                                                                              >  5:282759/1‑145 (MQ=255)
  aCAGGCGGTTTGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGccc                                                                                                                                           >  8:285781/1‑147 (MQ=255)
  aCAGGCGGTTTGCCTAACACAAAACCGACAACTTTATTTTTATCGACCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGccc                                                                                                                                           >  6:256815/1‑147 (MQ=255)
         gTTTGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTaaa                                                                                                                                   <  8:51500/148‑1 (MQ=255)
          tttGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAAt                                                                                                                                  <  1:159064/148‑1 (MQ=255)
          tttGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAAt                                                                                                                                  >  8:257211/1‑148 (MQ=255)
            tGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAAtt                                                                                                                                 <  3:202250/147‑1 (MQ=255)
                tAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACAGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCg                                                                                                                            <  5:256815/148‑1 (MQ=255)
                 aaCACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGa                                                                                                                           >  4:355925/1‑148 (MQ=255)
                  acacAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCg                                                                                                                            <  3:76423/146‑1 (MQ=255)
                  acacAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCg                                                                                                                            <  7:257211/146‑1 (MQ=255)
                  acacAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCg                                                                                                                            <  7:285781/146‑1 (MQ=255)
                  acacAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGGGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCg                                                                                                                            <  4:182724/146‑1 (MQ=255)
                      aaaaCCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCa                                                                                                                      <  3:355925/148‑1 (MQ=255)
                         agcGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCa                                                                                                                      <  5:334653/143‑1 (MQ=255)
                         aCCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACAATGGCCACCGCCCAATTTAAATTTGGCGACGCCa                                                                                                                      >  6:334653/1‑145 (MQ=255)
                                  tttATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAAt                                                                                                                                  <  5:135642/124‑1 (MQ=255)
                                  tttATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAAt                                                                                                                                  >  6:135642/1‑124 (MQ=255)
                                    tATTTTTATCGCCCGGATGACCGATTCCGATGCGTACACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATACAg                                                                                                        >  5:302550/1‑148 (MQ=255)
                                           aTCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTtcgtc                                                                                                  <  2:221331/147‑1 (MQ=255)
                                            tCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTtcgtcg                                                                                                 <  6:345763/147‑1 (MQ=255)
                                                    aTGCCCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGTCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATcc                                                                                                          <  1:352974/130‑1 (MQ=255)
                                                    aTGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGTCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATcc                                                                                                          <  1:347493/130‑1 (MQ=255)
                                                    aTGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGTCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATcc                                                                                                          >  2:352974/1‑130 (MQ=255)
                                                    aTGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGTCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATcc                                                                                                          >  2:347493/1‑130 (MQ=255)
                                                       aCCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTtcgtc                                                                                                  >  6:289898/1‑135 (MQ=255)
                                                       aCCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTtcgtc                                                                                                  <  5:289898/135‑1 (MQ=255)
                                                          gATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAAtt                                                                                  <  2:12970/148‑1 (MQ=255)
                                                               cGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGGCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGt                                                                              <  2:310919/147‑1 (MQ=255)
                                                                gATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCag                                                                                                                     >  8:2922/1‑107 (MQ=255)
                                                                gATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCag                                                                                                                     <  7:2922/107‑1 (MQ=255)
                                                                     gTAAACGGTGAAAGATAGGCTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGAtt                                                                       <  6:102048/148‑1 (MQ=255)
                                                                                 aGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTTGGCCACCAGAATTTCGTCCTGATTAATGCGGaaaaa                                                           <  6:302550/148‑1 (MQ=255)
                                                                                 aGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGCGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCAACAGAATTTCGTCCGGATTAATGCGGaaaaa                                                           <  5:150515/148‑1 (MQ=255)
                                                                                 aGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGCGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCAACAGAATTTCGTCCGGATTAATGCGGaaaa                                                            >  6:150515/1‑147 (MQ=255)
                                                                                          tattaCCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTa                                                                      >  4:61251/1‑128 (MQ=255)
                                                                                          tattaCCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTa                                                                      <  3:61251/128‑1 (MQ=255)
                                                                                           attaCCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGcc                                                    >  6:77240/1‑145 (MQ=255)
                                                                                           attaCCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGcc                                                    <  5:77240/145‑1 (MQ=255)
                                                                                              aCCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATcgcc                                              >  2:298550/1‑148 (MQ=255)
                                                                                                  aaTTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCaa                                          >  4:106292/1‑148 (MQ=255)
                                                                                                  aaTTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCaa                                          >  6:275398/1‑148 (MQ=255)
                                                                                                    tttACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGaaaaa                                                           >  5:203698/1‑129 (MQ=255)
                                                                                                    tttACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGaaaaa                                                           <  6:203698/129‑1 (MQ=255)
                                                                                                     ttACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGc                                                     >  3:80392/1‑134 (MQ=255)
                                                                                                     ttACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACgg                                       >  3:274291/1‑148 (MQ=255)
                                                                                                     ttACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGGGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGc                                                     <  4:80392/134‑1 (MQ=255)
                                                                                                     ttACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATACACTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACgg                                       >  5:235652/1‑148 (MQ=255)
                                                                                                        cTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCt                                     <  4:274291/147‑1 (MQ=255)
                                                                                                          gatgatGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGCTCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTg                                  >  7:116343/1‑148 (MQ=255)
                                                                                                           agtatgTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCTGAAAAAACTGGCCATCGCCGCAACGGCttt                                   <  4:211552/143‑1 (MQ=255)
                                                                                                                     tCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTgaga                         <  3:467/146‑1 (MQ=255)
                                                                                                                     tCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTgaga                         <  3:106292/146‑1 (MQ=255)
                                                                                                                         tCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTGAGATTCAt                    >  3:89290/1‑147 (MQ=255)
                                                                                                                         tCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTGAGATTCAt                    >  3:180679/1‑147 (MQ=255)
                                                                                                                             ttGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTGAGATTCATAAATg               <  5:275398/148‑1 (MQ=255)
                                                                                                                              tgtgACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGAtt                                                                       <  7:332045/91‑1 (MQ=255)
                                                                                                                              tgtgACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGAtt                                                                       >  8:332045/1‑91 (MQ=255)
                                                                                                                                        aTGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTGAGATTCATAAATGTAGTCGGGAc     >  3:249204/1‑147 (MQ=255)
                                                                                                                                        aTGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTGAGATTCATAAATGTAGTCGGGAc     >  4:46008/1‑147 (MQ=255)
                                                                                                                                         tGGCCACCGCCCAATTTAATTTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGATAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTGAGATTCATAAATGTAGTCGGGACTa   <  6:235652/148‑1 (MQ=255)
                                                                                                                                         tGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCTTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTGAGATTCATAAATGTAGTCGGGACTa   <  8:116343/148‑1 (MQ=255)
                                                                                                                                           gCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTGAGATTCATAAATGTAGTCGGGACTaa  <  1:298550/147‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
TAACAGGCGGTTTGCCTAACACAAAACCGACAACTTTATTTTTATCGCCCGGATGACCGATTCCGATGCGTAAACGGTGAAAGTTAGGGTTATTACCCAATTTACTGATGATGTCTTTCAGTCCATTGTGACCACCATGGCCACCGCCCAATTTAAATTTGGCGACGCCAGGAGGCAGATCCAGTTCGTCGTGGGCCACCAGAATTTCGTCCGGATTAATGCGGAAAAAACTGGCCATCGCCGCAACGGCTTTGCCGCTGAGATTCATAAATGTAGTCGGGACTAA  >  NC_000913/1258045‑1258330

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