Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,926,894 T→G 8.1% D649A (GAC→GCC)  ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,926,8940TG 8.1% 289.8 / 17.5 86D649A (GAC→GCC) ptrBprotease II
Reads supporting (aligned to +/- strand):  ref base T (41/38);  new base G (4/3);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

GTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTT  >  NC_000913/1926748‑1927042
                                                                                                                                                  |                                                                                                                                                    
gTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTcc                                                                                                                                                    <  2:1073713/149‑1 (MQ=255)
         tGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCAt                                                                                                                                                  >  1:690560/1‑142 (MQ=255)
         tGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCAt                                                                                                                                                  <  2:690560/142‑1 (MQ=255)
         tGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                                                                           >  4:698236/1‑149 (MQ=255)
                   tCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATa                                                                                                                                 >  2:1142403/1‑149 (MQ=255)
                        aaaaTACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCa                                                                                                                                                   >  2:680803/1‑126 (MQ=255)
                        aaaaTACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCa                                                                                                                                                   <  1:680803/126‑1 (MQ=255)
                               ttAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                                                                           >  2:582064/1‑127 (MQ=255)
                               ttAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                                                                           <  1:582064/127‑1 (MQ=255)
                               ttAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAt                                                                                                                     <  3:698236/149‑1 (MQ=255)
                                    ccGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGtt                                                                                                                <  1:1142403/149‑1 (MQ=255)
                                          gCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                                                                                          >  4:1022279/1‑149 (MQ=255)
                                                cGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                                                                                                          <  1:523126/127‑1 (MQ=255)
                                                cGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTCCGTACGATTGAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                                                                                                          >  2:523126/1‑127 (MQ=255)
                                                    tAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAtt                                                                                                <  3:1022279/149‑1 (MQ=255)
                                                    tAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATCTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCAt                                                                                                                                                  >  2:782105/1‑99 (MQ=255)
                                                    tAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATCTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCAt                                                                                                                                                  <  1:782105/99‑1 (MQ=255)
                                                         ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg                                                                                           >  3:497611/1‑149 (MQ=255)
                                                            ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                                                                           <  1:210008/98‑1 (MQ=255)
                                                            ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                                                                           >  2:210008/1‑98 (MQ=255)
                                                                     gcgcGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGcc                                                                               <  4:497611/149‑1 (MQ=255)
                                                                         gACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGtt                                                                           >  1:30459/1‑149 (MQ=255)
                                                                            cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                                                                        >  3:1145199/1‑149 (MQ=255)
                                                                                aaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGc                                                                                            >  1:1069200/1‑125 (MQ=255)
                                                                                aaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGc                                                                                            <  2:1069200/125‑1 (MQ=255)
                                                                                 aaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGcc                                                                               <  2:595230/137‑1 (MQ=255)
                                                                                 aaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGcc                                                                               >  1:595230/1‑137 (MQ=255)
                                                                                 aaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                                                                        <  3:1330462/144‑1 (MQ=255)
                                                                                 aaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                                                                        >  4:1330462/1‑144 (MQ=255)
                                                                                   gCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg                                                                                           <  3:153084/123‑1 (MQ=255)
                                                                                   gCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg                                                                                           >  4:153084/1‑123 (MQ=255)
                                                                                    cATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATc                                                                                                                    <  4:497758/97‑1 (MQ=255)
                                                                                    cATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATc                                                                                                                    >  3:497758/1‑97 (MQ=255)
                                                                                    cATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc                                                                >  2:1465600/1‑149 (MQ=255)
                                                                                     atatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCa                                                               >  4:1357141/1‑149 (MQ=255)
                                                                                     atatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCa                                                               <  2:30459/149‑1 (MQ=255)
                                                                                        ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                                                                              <  1:207446/83‑1 (MQ=255)
                                                                                        ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                                                                              >  2:207446/1‑83 (MQ=255)
                                                                                        ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCt                                                            >  2:1254568/1‑149 (MQ=255)
                                                                                          ccATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTga                                                          >  1:515493/1‑149 (MQ=255)
                                                                                            aTCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGcc                                                                               <  1:1127129/126‑1 (MQ=255)
                                                                                            aTCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGcc                                                                               >  2:1127129/1‑126 (MQ=255)
                                                                                               gCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATc                                                     <  3:1357141/149‑1 (MQ=255)
                                                                                                 tACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGt                                                   <  2:515493/149‑1 (MQ=255)
                                                                                                 tACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGt                                                   <  1:1465600/149‑1 (MQ=255)
                                                                                                    gCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCa                                                >  1:1399175/1‑149 (MQ=255)
                                                                                                      cTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTa                                                                                              <  2:205645/101‑1 (MQ=255)
                                                                                                      cTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTa                                                                                              >  1:205645/1‑101 (MQ=255)
                                                                                                       ttCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAAc                                             >  2:954428/1‑149 (MQ=255)
                                                                                                           tACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                                                                                                          <  1:1287444/68‑1 (MQ=255)
                                                                                                           tACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                                                                                                          >  2:1287444/1‑68 (MQ=255)
                                                                                                            aCGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTc                                        <  4:1145199/149‑1 (MQ=255)
                                                                                                               aTTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACa                                                                                                                                         >  3:578699/1‑49 (MQ=255)
                                                                                                               aTTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACa                                                                                                                                         <  4:578699/49‑1 (MQ=255)
                                                                                                                 ttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAc                                   <  2:1399175/149‑1 (MQ=255)
                                                                                                                         aCCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAc                                   <  2:672025/141‑1 (MQ=255)
                                                                                                                         aCCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAc                                   >  1:672025/1‑141 (MQ=255)
                                                                                                                          ccAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTga                                                          <  3:805551/117‑1 (MQ=255)
                                                                                                                          ccAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTga                                                          >  4:805551/1‑117 (MQ=255)
                                                                                                                             tctgtACCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCgg                       >  1:963461/5‑149 (MQ=255)
                                                                                                                              aTTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGa                      <  1:954428/149‑1 (MQ=255)
                                                                                                                                    cgacATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGcc                >  1:1226621/4‑149 (MQ=255)
                                                                                                                                     gacATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTga                                                          >  4:1035007/3‑106 (MQ=255)
                                                                                                                                     gacATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTga                                                          <  3:1035007/104‑1 (MQ=255)
                                                                                                                                     gacATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTg                                                  >  3:1001573/3‑114 (MQ=255)
                                                                                                                                     gacATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTg                                                  <  4:1001573/112‑1 (MQ=255)
                                                                                                                                      ccATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCaaa                                              <  2:1154736/117‑1 (MQ=255)
                                                                                                                                      ccATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCaaa                                              >  1:1154736/1‑117 (MQ=255)
                                                                                                                                      ccATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCa                                 >  3:385673/1‑130 (MQ=255)
                                                                                                                                      ccATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCa                                 <  4:385673/130‑1 (MQ=255)
                                                                                                                                       cATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCtgt             <  2:963461/149‑1 (MQ=255)
                                                                                                                                          ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTTCACCTGAGAATCGt                                                   <  4:945404/108‑1 (MQ=255)
                                                                                                                                          ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTTCACCTGAGAATCGt                                                   >  3:945404/1‑108 (MQ=255)
                                                                                                                                          ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCg          <  1:1254568/149‑1 (MQ=255)
                                                                                                                                           gCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTc                                        >  4:1156922/1‑118 (MQ=255)
                                                                                                                                           gCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCgg         >  3:629923/1‑149 (MQ=255)
                                                                                                                                           gCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGCGCAAACCGGTc                                        <  3:1156922/118‑1 (MQ=255)
                                                                                                                                             cTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTg       >  1:1125417/1‑149 (MQ=255)
                                                                                                                                             cTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTg       >  1:882244/1‑149 (MQ=255)
                                                                                                                                             cTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTg       >  1:865645/1‑149 (MQ=255)
                                                                                                                                              tcaGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCg                                                                              <  4:956530/75‑1 (MQ=255)
                                                                                                                                              tGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCg                                                                              >  3:956530/1‑77 (MQ=255)
                                                                                                                                                 gTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCGTTTTGCCGGTTCCCAata                                                                    >  3:720883/1‑84 (MQ=255)
                                                                                                                                                 gTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCGTTTTGCCGGTTCCCAata                                                                    <  4:720883/84‑1 (MQ=255)
                                                                                                                                                 gTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCa                                                               <  4:1288739/89‑1 (MQ=255)
                                                                                                                                                 gTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCa                                                               >  3:1288739/1‑89 (MQ=255)
                                                                                                                                                  tCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGtt  <  2:1125417/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTT  >  NC_000913/1926748‑1927042

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