Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 741,486 G→A 5.5% A357A (GCC→GCT dtpD ← dipeptide and tripeptide permease D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913741,4860GA 5.5% 226.0 / 2.6 72A357A (GCC→GCTdtpDdipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base A (2/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

CGTAATTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATCGCGGT  >  NC_000913/741348‑741629
                                                                                                                                          |                                                                                                                                               
cgTAATTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCtcatc                                                                                                                                        <  7:392565/148‑1 (MQ=255)
cgTAATTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCtcat                                                                                                                                         >  7:332348/1‑147 (MQ=255)
cgTAATTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGAGCTTAAGGTCAGAATGCAGAATCCGGCGCtcat                                                                                                                                         >  1:444009/1‑147 (MQ=255)
    aTTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCaggc                                                                                                                                    >  1:370057/1‑148 (MQ=255)
      ttGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCTCCATGCCTAATACCATCAGTGGCAGAGAAGAGTGCCCATACATTGAGGACCAGCGGGCGCTTAAGGGCAGAATGCAGAATCCGGCGATCATCaggac                                                                                                                                   >  6:193771/1‑145 (MQ=255)
       tGCTACATGGCAACCGGGTCGATAAACAGGTCCGCAAAGCCCATCACCGCCAGGCCTAATCCCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATCCAGAATCCGGCGCTCATCaggcca                                                                                                                                  <  7:154897/147‑1 (MQ=255)
         cGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCaggccaag                                                                                                                                >  2:80207/1‑147 (MQ=255)
         cGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCaggccaag                                                                                                                                <  1:80207/147‑1 (MQ=255)
              tGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGgccaaggccaag                                                                                                                          <  8:332348/148‑1 (MQ=255)
               ggCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCAGGCGCTCATCAggccaaggccaa                                                                                                                           >  7:397827/1‑146 (MQ=255)
                     cGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCCTCTGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAAt                                                                                                                   <  2:421348/148‑1 (MQ=255)
                         tCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               >  3:383324/1‑148 (MQ=255)
                         tCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               <  8:205623/148‑1 (MQ=255)
                                       cGCAAAGCCCATCACCGCCCGGCCTTTTACCATCAGTGGCAGAGTAGAGTGACCATAAATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTGATCAGTCCAAGGCCAAGAGCAAATTTCCCCCAGAGgcgc                                                                                                    <  6:387654/145‑1 (MQ=255)
                                             gCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGa                                                                                           <  8:174360/148‑1 (MQ=255)
                                                 aTCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTAc                                                                                       >  1:23382/1‑148 (MQ=255)
                                                  tCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTAcc                                                                                      >  5:47955/1‑148 (MQ=255)
                                                  tCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTAcc                                                                                      >  5:380187/1‑148 (MQ=255)
                                                  tCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTAcc                                                                                      >  4:15617/1‑148 (MQ=255)
                                                  tCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAACAGCAAATTTCCCCCAGATGCGCACGGTACGATTAcc                                                                                      >  7:133198/1‑148 (MQ=255)
                                                    aCCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCcg                                                                                    <  2:370057/148‑1 (MQ=255)
                                                     ccGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCcgc                                                                                   >  3:166279/1‑148 (MQ=255)
                                                     ccGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTCCGGACCAGCGGGCGCATAAGGTCAGAATGCACAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCAAGATGCGCACGCTACGATTACCcgc                                                                                   >  3:269209/1‑148 (MQ=255)
                                                          aGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGaca                                                                               <  1:134195/147‑1 (MQ=255)
                                                          aGGCCTAATACCATCAGAGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGaca                                                                               >  2:134195/1‑147 (MQ=255)
                                                              cTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCCAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACActctc                                                                          <  8:156439/148‑1 (MQ=255)
                                                                aaTACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               >  6:335404/1‑109 (MQ=255)
                                                                aaTACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               <  5:335404/109‑1 (MQ=255)
                                                                aaTACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCAGCGCtc                                                                                                                                           <  1:9036/81‑1 (MQ=255)
                                                                aaTACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCAGCGCtc                                                                                                                                           >  2:9036/1‑81 (MQ=255)
                                                                 aTACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACActctct                                                                         <  8:72350/146‑1 (MQ=255)
                                                                  tACCCTCAGTCGCAGAGAAGCGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCtttt                                                                      <  5:193771/148‑1 (MQ=255)
                                                                    ccATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTac                                                                    <  4:166279/148‑1 (MQ=255)
                                                                       tCAGTGGCAGAGAAGAGTGACCATTCATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCAGCGCTCATCaggcc                                                                                                                                   >  6:185807/1‑82 (MQ=255)
                                                                       tCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCAGCGCTCATCaggcc                                                                                                                                   <  5:185807/82‑1 (MQ=255)
                                                                        cAGTGGCAGAGAAGAGTGACCATACATTGCGGACCATCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTaccacc                                                                >  1:15524/1‑148 (MQ=255)
                                                                        cAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCcg                                                                                    <  8:110422/128‑1 (MQ=255)
                                                                        cAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCcg                                                                                    >  7:110422/1‑128 (MQ=255)
                                                                           tGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTAccacccac                                                             >  1:126958/1‑148 (MQ=255)
                                                                            ggCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCACGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTGTTACCACCCACg                                                            >  6:335353/1‑148 (MQ=255)
                                                                            ggCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAg                                                                                                          >  5:31502/1‑102 (MQ=255)
                                                                            ggCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAg                                                                                                          <  6:31502/102‑1 (MQ=255)
                                                                               agagAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCa                                                         <  2:23382/148‑1 (MQ=255)
                                                                                 agaagaGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCaaa                                                                                                                    >  4:244068/1‑87 (MQ=255)
                                                                                 agaagaGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCaaa                                                                                                                    <  3:244068/87‑1 (MQ=255)
                                                                                  gaagaGTGACCATACATTGCGGACCATCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGa                                                      <  2:15524/148‑1 (MQ=255)
                                                                                  gaagaGTGACCATACATTGCGGACCAGCGTGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGa                                                      <  3:247361/148‑1 (MQ=255)
                                                                                  gaagaGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCtt                                                                        >  7:162166/1‑130 (MQ=255)
                                                                                  gaagaGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCtt                                                                        <  8:162166/130‑1 (MQ=255)
                                                                                   aagaGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGaa                                                     >  8:182135/1‑148 (MQ=255)
                                                                                      agTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAacac                                                  >  7:61765/1‑148 (MQ=255)
                                                                                       gTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGTCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAacac                                                  <  2:126958/147‑1 (MQ=255)
                                                                                            cATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGca                                            <  4:383324/148‑1 (MQ=255)
                                                                                               aCATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGAAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGAAACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAg                                         >  8:422145/1‑148 (MQ=255)
                                                                                                cATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAg                                         >  7:238047/1‑147 (MQ=255)
                                                                                                   tGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGt                                                                                               >  2:230283/1‑90 (MQ=255)
                                                                                                   tGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGt                                                                                               <  1:230283/90‑1 (MQ=255)
                                                                                                    gCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATa                                     <  6:380187/147‑1 (MQ=255)
                                                                                                          cAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACTGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGa                              <  7:422145/148‑1 (MQ=255)
                                                                                                           aGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGa                              <  5:385946/147‑1 (MQ=255)
                                                                                                           aGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGa                              >  6:385946/1‑147 (MQ=255)
                                                                                                              gggCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCttt                                                                       >  3:254692/1‑103 (MQ=255)
                                                                                                              gggCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCttt                                                                       <  4:254692/103‑1 (MQ=255)
                                                                                                                     tAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCCCGCCCGGCACACACCGCACAGCATAACTGCGAAGGCCTTAATc                   >  8:379284/1‑148 (MQ=255)
                                                                                                                                tGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACAt        <  7:182135/148‑1 (MQ=255)
                                                                                                                                    gAATCCGGCGCTCATCAGGCCACGGCCAAGAGAAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCCACACTCGCTTTTACCACCCACGCCAGTAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATcgcg    >  2:167427/1‑148 (MQ=255)
                                                                                                                                    gAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATcgcg    >  8:250313/1‑148 (MQ=255)
                                                                                                                                      aTCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGa                              >  1:301375/1‑120 (MQ=255)
                                                                                                                                      aTCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGa                              <  2:301375/120‑1 (MQ=255)
                                                                                                                                      aTCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATCGCGGt  >  3:150985/1‑148 (MQ=255)
                                                                                                                                        ccGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               >  7:383910/1‑37 (MQ=255)
                                                                                                                                        ccGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               <  8:383910/37‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                               
CGTAATTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATCGCGGT  >  NC_000913/741348‑741629

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