Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,484,779 A→G 5.8% H391H (CAT→CAC pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,484,7790AG 5.8% 220.0 / 2.4 67H391H (CAT→CACpepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base A (15/48);  new base G (2/2);  total (17/50)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCA  >  NC_000913/4484641‑4484914
                                                                                                                                          |                                                                                                                                       
cATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAatat                                                                                                                                        >  1:4841/1‑140 (MQ=255)
 tgctCGGCATCATTGGACTCCAGTTGGTGCTGTGCCTCGTGCCCCCGCGGTAAGCGCCAGGCGCGGTCACCGGTTTGTTCAGACGCGGCTATCAGTTCATGGGCCAGGGGATTATGGTTCGCCATCAGACCAGTCATctg                                                                                                                                       <  2:1200249/137‑1 (MQ=255)
       gCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATGAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGAcc                                                                                                                                 >  2:1071521/1‑140 (MQ=255)
       gCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGAcc                                                                                                                                 >  2:1071550/1‑140 (MQ=255)
            aTTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAgcg                                                                                                                            <  2:927839/140‑1 (MQ=255)
            aTTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAgcg                                                                                                                            <  2:927838/140‑1 (MQ=255)
            aTTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAgcg                                                                                                                            <  1:134317/140‑1 (MQ=255)
            aTTGGACTCCAGTCGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTACAGACGCGGCAAACAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAgcg                                                                                                                            <  1:134401/140‑1 (MQ=255)
                     cAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGc                                                                                                                   <  1:29916/140‑1 (MQ=255)
                     cAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGc                                                                                                                   <  1:28616/140‑1 (MQ=255)
                     cAGTTGTTCCTGAAACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGATCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGCTGACCCAGCGCGATCACGc                                                                                                                   <  1:29966/140‑1 (MQ=255)
                           ttCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCAc                                                                                                             >  1:749274/1‑140 (MQ=255)
                             ccTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCg                                                                                                           >  2:690101/1‑140 (MQ=255)
                             ccTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCg                                                                                                           >  2:540671/1‑140 (MQ=255)
                                       gTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCACTCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCg                                                                                                 >  2:272557/1‑140 (MQ=255)
                                       gTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCg                                                                                                 >  2:411129/1‑140 (MQ=255)
                                       gTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCg                                                                                                 >  2:272556/1‑140 (MQ=255)
                                         cACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGcc                                                                                               <  1:86152/140‑1 (MQ=255)
                                         cACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGcc                                                                                               <  1:86428/140‑1 (MQ=255)
                                         cACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGcc                                                                                               <  1:377073/140‑1 (MQ=255)
                                         cACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGcc                                                                                               <  2:754852/140‑1 (MQ=255)
                                          aCCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCa                                                                                              <  2:859255/140‑1 (MQ=255)
                                                                      ccGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACg                                                                  <  2:106370/140‑1 (MQ=255)
                                                                      ccGGATTGTTCAGACGCGGCAATCAGTTCAAGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACg                                                                  <  2:106455/140‑1 (MQ=255)
                                                                        ggATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCt                                                                <  1:787276/140‑1 (MQ=255)
                                                                           ttGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCaa                                                             <  2:1186516/140‑1 (MQ=255)
                                                                             gTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACg                                                           <  2:700427/140‑1 (MQ=255)
                                                                                       ggCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAaca                                                 <  2:844635/140‑1 (MQ=255)
                                                                                        gCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCGGTAATGTGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAacac                                                >  1:1093054/1‑140 (MQ=255)
                                                                                        gCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCGGTAATGTGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAacac                                                >  1:1093244/1‑140 (MQ=255)
                                                                                        gCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAacac                                                >  1:22437/1‑140 (MQ=255)
                                                                                        gCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAacac                                                >  1:528076/1‑140 (MQ=255)
                                                                                         cAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGCCCAGTACTATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACg                                               <  1:239321/140‑1 (MQ=255)
                                                                                         cAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACg                                               <  1:341280/140‑1 (MQ=255)
                                                                                         cAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACg                                               <  1:1198920/140‑1 (MQ=255)
                                                                                         cAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACg                                               <  1:1198935/140‑1 (MQ=255)
                                                                                            tCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCg                                            <  1:1034614/140‑1 (MQ=255)
                                                                                            tCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCg                                            <  2:579793/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTACTATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:47764/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:1062911/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  2:324536/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:438568/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:1153807/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:12006/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:175920/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:207030/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:909727/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:486771/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:418466/140‑1 (MQ=255)
                                                                                              aGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGca                                          <  1:438517/140‑1 (MQ=255)
                                                                                                        ccAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAgg                                <  2:791551/140‑1 (MQ=255)
                                                                                                         cAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAggc                               <  2:297742/140‑1 (MQ=255)
                                                                                                             ggATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGcc                           <  2:188522/140‑1 (MQ=255)
                                                                                                                     gTTCGCCATCAGACCGGTAATGTGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGt                   <  2:1093244/140‑1 (MQ=255)
                                                                                                                     gTTCGCCATCAGACCGGTAATGTGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGt                   <  2:1093054/140‑1 (MQ=255)
                                                                                                                     gTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGt                   <  2:964745/140‑1 (MQ=255)
                                                                                                                     gTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGt                   <  1:889649/140‑1 (MQ=255)
                                                                                                                       tCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACATTACCAGGCGGCCTTCAGCGTCg                 <  1:518648/140‑1 (MQ=255)
                                                                                                                       tCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCg                 <  1:508127/140‑1 (MQ=255)
                                                                                                                       tCGCCATCAGACCAGATATATGATGACCCAGCGCGATCACGCAGGCACCGGTCCGCGTCGCCACGACAATCACCGCTTCCGGCTCAAAACGCACAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCg                 <  2:1209315/140‑1 (MQ=255)
                                                                                                                        cGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAgg                                <  2:353243/124‑1 (MQ=255)
                                                                                                                        cGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAgg                                >  1:353243/1‑124 (MQ=255)
                                                                                                                         gCCATCAGACCAGTAATATGTTGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCATCGTAAGTTATCACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGgt               >  2:27820/1‑140 (MQ=255)
                                                                                                                            aTCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGtt            <  2:517930/140‑1 (MQ=255)
                                                                                                                            aTCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGtt            <  2:515725/140‑1 (MQ=255)
                                                                                                                               aGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAg         <  2:1059038/140‑1 (MQ=255)
                                                                                                                                     gTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTc   >  1:486352/1‑140 (MQ=255)
                                                                                                                                      tAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCa  >  2:735667/1‑140 (MQ=255)
                                                                                                                                      tAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGCTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCa  >  2:544396/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                       
CATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCA  >  NC_000913/4484641‑4484914

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