Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,802,536 C→T 5.9% R12H (CGC→CAC)  thrS ← threonyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,802,5360CT 5.9% 220.2 / 2.8 68R12H (CGC→CAC) thrSthreonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (35/29);  new base T (2/2);  total (37/31)
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.55e-01

TAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCCAACAAATGCTAAATTGGT  >  NC_000913/1802390‑1802680
                                                                                                                                                  |                                                                                                                                                
tAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGt                                                                                                                                                >  4:165245/1‑149 (MQ=255)
    gATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGc                                                                                                                                            >  3:209576/1‑149 (MQ=255)
     aTCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCt                                                                                                                                           >  2:211420/1‑149 (MQ=255)
     aTCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCt                                                                                                                                           <  4:108558/149‑1 (MQ=255)
     aTCGACAGCTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATAGCGTTGGCt                                                                                                                                           >  3:465407/1‑149 (MQ=255)
             ttGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAg                                                                                                                                   <  4:209576/149‑1 (MQ=255)
              tgtgCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAgg                                                                                                                                  <  1:484508/149‑1 (MQ=255)
                  cgtcgtTTTCAAACAGATCGCAAGCATCAAACAGTTCGCCATTAACGCGCCCCGCGATACAGGCTTTCGGCAGACCTGGACCAATGTCCAGGGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAaga                                                                                                                              <  1:286989/149‑1 (MQ=255)
                   gtcgtTTTCAATCAGATCGCAAGCATGAACCCGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATTTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAagag                                                                                                                             >  2:489141/1‑149 (MQ=255)
                    tcgtTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGt                                                                                                                            >  4:7019/1‑149 (MQ=255)
                      gtTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTta                                                                                                                          <  3:427572/149‑1 (MQ=255)
                      gtTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTta                                                                                                                          >  3:57174/1‑149 (MQ=255)
                       ttttCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTtat                                                                                                                         <  2:110214/149‑1 (MQ=255)
                       ttttCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTtat                                                                                                                         >  2:335955/1‑149 (MQ=255)
                         ttCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATaa                                                                                                                       >  3:492461/1‑149 (MQ=255)
                            aaTCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAg                                                                                                                    <  3:629/149‑1 (MQ=255)
                                aGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCAt                                                                                                                <  3:533766/149‑1 (MQ=255)
                                aGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCAt                                                                                                                >  2:159359/1‑149 (MQ=255)
                                 gATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCAtt                                                                                                               <  1:159359/149‑1 (MQ=255)
                                         gCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATcc                                                                                                       >  1:239623/1‑149 (MQ=255)
                                              aaCCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTAtt                                                                                                  >  1:191785/1‑149 (MQ=255)
                                              aaCCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTAtt                                                                                                  <  1:577959/149‑1 (MQ=255)
                                              aaCCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTAtt                                                                                                  >  1:59390/1‑149 (MQ=255)
                                                          aTTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGAc                                                                                      <  4:57174/149‑1 (MQ=255)
                                                             aaCGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCca                                                                                   <  2:211393/149‑1 (MQ=255)
                                                              aCGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAg                                                                                             <  2:55575/138‑1 (MQ=255)
                                                              aCGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAg                                                                                             >  1:55575/1‑138 (MQ=255)
                                                              aCGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCcac                                                                                  <  2:191785/149‑1 (MQ=255)
                                                              aCGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCcac                                                                                  >  2:579468/1‑149 (MQ=255)
                                                                       gCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGAt                                                                         <  1:579468/149‑1 (MQ=255)
                                                                       gCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGAt                                                                         >  3:357334/1‑149 (MQ=255)
                                                                        cGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATc                                                                        <  2:59390/149‑1 (MQ=255)
                                                                           tACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATcaca                                                                     >  2:312462/1‑149 (MQ=255)
                                                                            aCAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACAt                                                                    >  4:404408/1‑149 (MQ=255)
                                                                              aGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATAc                                                                  >  2:381896/1‑149 (MQ=255)
                                                                              aGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATAc                                                                  >  3:114280/1‑149 (MQ=255)
                                                                                  tttCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAg                                                              >  4:3312/1‑149 (MQ=255)
                                                                                        cAGACCTGGACCAACGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAAt                                                        >  3:41600/1‑149 (MQ=255)
                                                                                          gACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAAttt                                                      <  2:239623/149‑1 (MQ=255)
                                                                                          gACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAAttt                                                      <  4:357334/149‑1 (MQ=255)
                                                                                                gACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTttat                                                >  3:69289/1‑149 (MQ=255)
                                                                                                gACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTttat                                                >  2:211874/1‑149 (MQ=255)
                                                                                                gACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAAGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTttat                                                >  4:347492/1‑149 (MQ=255)
                                                                                                  ccAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTttatta                                              <  3:353834/149‑1 (MQ=255)
                                                                                                   cAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTaa                                             >  3:150660/1‑149 (MQ=255)
                                                                                                    aaTGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAAc                                            <  4:438233/149‑1 (MQ=255)
                                                                                                     aTGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACa                                           <  1:312462/149‑1 (MQ=255)
                                                                                                          cAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCa                                      >  2:514692/1‑149 (MQ=255)
                                                                                                           aGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCAt                                     >  3:138476/1‑149 (MQ=255)
                                                                                                            gcgcAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCAtt                                    <  1:514692/149‑1 (MQ=255)
                                                                                                                  cATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGaaa                              <  4:137107/149‑1 (MQ=255)
                                                                                                                   aTCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAAc                             >  3:139318/1‑149 (MQ=255)
                                                                                                                      cATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCAt                          <  1:250288/149‑1 (MQ=255)
                                                                                                                       aTGGGGCTTACAGCGTGATCGTAATGGTGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATc                         >  4:72119/1‑149 (MQ=255)
                                                                                                                       aTGGGGCTTACAGCGTGATCGTAATGGTGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATc                         >  4:253778/1‑149 (MQ=255)
                                                                                                                         ggggCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTa                       <  1:211874/149‑1 (MQ=255)
                                                                                                                           ggCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGc                     >  1:449013/1‑149 (MQ=255)
                                                                                                                            gCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGcc                    <  2:430061/149‑1 (MQ=255)
                                                                                                                            gCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGcc                    >  1:430061/1‑149 (MQ=255)
                                                                                                                                aCAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCcaaca                >  4:248058/1‑149 (MQ=255)
                                                                                                                                aCAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCcaaca                >  3:251607/1‑149 (MQ=255)
                                                                                                                                   gCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCCAACAAAt             <  4:41600/149‑1 (MQ=255)
                                                                                                                                           gTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCCAACAAATGCTAAAtt     <  1:335955/149‑1 (MQ=255)
                                                                                                                                           gTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCCAACAAATGCTAAAtt     <  4:69289/149‑1 (MQ=255)
                                                                                                                                            tAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCCAACAAATGCTAAATTg    >  2:240652/1‑149 (MQ=255)
                                                                                                                                              aTGGTGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCCAACAAATGCTAAATTGGt  <  3:72119/149‑1 (MQ=255)
                                                                                                                                              aTGGTGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCCAACAAATGCTAAATTGGt  <  3:253778/149‑1 (MQ=255)
                                                                                                                                              aTGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCCAACAAATGCTAAATTGGt  >  4:433231/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATTAACGCGCCCTGCGATACAGGCTTTCGCCAGACCTGGACCAATGTCCAGCGCAACATCCATGGGGCTTACAGCGTGATCGTAATGGCGTTGGCTGCCATCAGGAAGAGTTATAACAGGCATTTTATATCCTTATTTGCAGTGGTGACCCACACGAAAGATCACATACAAAGAAAAATTTGTTTATTAACAGTTCATTGCGAAACCATCTAGCCAACAAATGCTAAATTGGT  >  NC_000913/1802390‑1802680

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: