Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,353,128 C→A V167V (GTG→GTT sapD ← antimicrobial peptide ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,353,1280CA98.7% 277.4 / NA 74V167V (GTG→GTTsapDantimicrobial peptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (39/34);  total (40/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 = 9.97e-01

GGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGCCTTTATAGGTCCA  >  NC_000913/1352980‑1353274
                                                                                                                                                    |                                                                                                                                                  
ggCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAAc                                                                                                                                                  <  1:470538/151‑1 (MQ=255)
 gCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACt                                                                                                                                                 >  1:157236/1‑151 (MQ=255)
    gATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACttt                                                                                                                                               <  2:718971/150‑1 (MQ=255)
       aaGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTg                                                                                                                                            <  2:1140035/150‑1 (MQ=255)
        aGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGa                                                                                                                                           >  1:707024/1‑150 (MQ=255)
        aGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGAc                                                                                                                                          >  1:504266/1‑151 (MQ=255)
         gCAAAATAGTGGCATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGAc                                                                                                                                          >  2:521926/1‑150 (MQ=255)
          cAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACAt                                                                                                                                        >  1:62431/1‑151 (MQ=255)
          cAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACAt                                                                                                                                        >  1:189902/1‑151 (MQ=255)
          cAAAATAGTGGTATTCCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGATGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGAc                                                                                                                                          >  2:243740/1‑149 (MQ=255)
             aaTAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTc                                                                                                                                      >  1:156103/1‑150 (MQ=255)
             aaTAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCa                                                                                                                                     <  1:142116/151‑1 (MQ=255)
             aaTAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCa                                                                                                                                     <  2:810055/151‑1 (MQ=255)
             aaTAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCa                                                                                                                                     <  2:281193/151‑1 (MQ=255)
             aaTAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCa                                                                                                                                     <  1:683113/151‑1 (MQ=255)
              aTAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCAc                                                                                                                                    >  1:900547/1‑151 (MQ=255)
                aGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCt                                                                                                                                  <  1:1132873/151‑1 (MQ=255)
                aGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCt                                                                                                                                  <  2:284308/151‑1 (MQ=255)
                 gTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCGccgt                                                                                                                                 <  2:289429/151‑3 (MQ=255)
                      aTTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGt                                                                                                                             <  2:690264/150‑1 (MQ=255)
                      aTTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGt                                                                                                                             <  1:831143/150‑1 (MQ=255)
                         gCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAAc                                                                                                                         >  1:553515/1‑151 (MQ=255)
                           tGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTc                                                                                                                       >  1:775769/1‑151 (MQ=255)
                            gTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACt                                                                                                                        >  1:629227/1‑149 (MQ=255)
                             ttttGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCAt                                                                                                                     >  2:644242/1‑151 (MQ=255)
                             ttttGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTACGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCAt                                                                                                                     >  2:457828/1‑151 (MQ=255)
                             ttgtGGTTGAGACGCGTCAGCAGGCAAAAGATTTGAGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCa                                                                                                                      <  2:10359/147‑1 (MQ=255)
                              tttGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCAt                                                                                                                     >  1:622361/1‑150 (MQ=255)
                                tGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAg                                                                                                                   <  1:264747/150‑1 (MQ=255)
                                 ggTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCCTAgg                                                                                                                  <  2:745931/150‑1 (MQ=255)
                                 ggTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAgg                                                                                                                  <  1:433340/150‑1 (MQ=255)
                                   ttGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGa                                                                                                                >  1:918336/1‑150 (MQ=255)
                                     gagaCGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCGCCTTCGGTCAACTCATAGGGaaa                                                                                                              <  2:770033/150‑1 (MQ=255)
                                      agaCGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGaaaa                                                                                                             >  1:566430/1‑150 (MQ=255)
                                      agaCGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAAc                                                                                                            >  2:1123871/1‑151 (MQ=255)
                                       gaCGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACt                                                                                                           >  2:558176/1‑151 (MQ=255)
                                        aCGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTg                                                                                                          <  1:26688/151‑1 (MQ=255)
                                           cgTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTgcgc                                                                                                       <  2:826109/151‑1 (MQ=255)
                                              cagcagGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgcatc                                                                                                    >  1:150197/1‑151 (MQ=255)
                                               agcagGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAAATGcgcatcg                                                                                                   >  2:84054/1‑151 (MQ=255)
                                                 cagGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgcatcgc                                                                                                  >  1:610125/1‑150 (MQ=255)
                                                  agGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgcatcgcat                                                                                                <  2:601696/151‑1 (MQ=255)
                                                  agGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgcatcgcat                                                                                                >  2:425169/1‑151 (MQ=255)
                                                       aaaGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTtg                                                                                           <  2:703797/151‑1 (MQ=255)
                                                          gATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTtgtg                                                                                         >  1:83317/1‑150 (MQ=255)
                                                            tttGCGCCTGGGCTGTTGGCTCCATTGAGTTGGCCGGTTCGTCAGCACTCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACACTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGAt                                                                                       >  2:598264/1‑150 (MQ=255)
                                                             ttGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCt                                                                                     >  1:423617/1‑151 (MQ=255)
                                                                    tGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATcc                                                                              >  1:1145104/1‑151 (MQ=255)
                                                                          gttgGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCAcg                                                                         <  2:85008/150‑1 (MQ=255)
                                                                          gttgGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCAcg                                                                         <  1:571688/150‑1 (MQ=255)
                                                                                 ccATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTgca                                                                  <  2:773695/150‑1 (MQ=255)
                                                                                           ggTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATg                                                       >  2:598608/1‑151 (MQ=255)
                                                                                             tCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGc                                                     >  2:450296/1‑151 (MQ=255)
                                                                                                gTTCATCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGGGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCAc                                                   >  2:1209456/1‑150 (MQ=255)
                                                                                                  tCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGc                                                >  1:231638/1‑151 (MQ=255)
                                                                                                   cGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCg                                               <  2:71862/151‑1 (MQ=255)
                                                                                                       aGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTg                                           <  2:984268/151‑1 (MQ=255)
                                                                                                       aGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTg                                           <  1:795877/151‑1 (MQ=255)
                                                                                                          aaTCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTgcgc                                        >  2:950288/1‑151 (MQ=255)
                                                                                                               gcagACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCgccagcc                                   <  2:821491/151‑1 (MQ=255)
                                                                                                                  gACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCaaa                                <  1:1047495/151‑1 (MQ=255)
                                                                                                                   aCGCGGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCaaaa                               <  1:324429/151‑1 (MQ=255)
                                                                                                                    cgcgGTTGATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAAc                              >  1:806320/1‑151 (MQ=255)
                                                                                                                           gATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGcca                       <  2:900160/151‑1 (MQ=255)
                                                                                                                           gATTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGcc                        >  1:635119/1‑150 (MQ=255)
                                                                                                                            aTTCGCCAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGccac                      >  2:375975/1‑151 (MQ=255)
                                                                                                                                 ccAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACg                 <  2:838144/151‑1 (MQ=255)
                                                                                                                                  cAGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACgg                >  2:301997/1‑151 (MQ=255)
                                                                                                                                   aGTGCAATGGCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGc               >  1:279297/1‑151 (MQ=255)
                                                                                                                                        aaTGGATATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGCCTTta          <  1:1020497/151‑1 (MQ=255)
                                                                                                                                           ggCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGCCTTtat         >  1:381264/1‑149 (MQ=255)
                                                                                                                                           ggCTATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGCCTTTATAg       <  1:56822/151‑1 (MQ=255)
                                                                                                                                            gCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCATCCAAAACGCTGCCACCAACGGCCTTTATAgg      >  1:1166590/1‑151 (MQ=255)
                                                                                                                                              tATCATAACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGCCTTTATAGGTc    >  1:976309/1‑151 (MQ=255)
                                                                                                                                               atcataACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGCCTTTATAGGt     <  2:294871/149‑1 (MQ=255)
                                                                                                                                               atcataACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGCCTTTATAGGTcc   <  1:633654/151‑1 (MQ=255)
                                                                                                                                                 cataaCTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGCCTTTATAGGTcc   >  2:52521/1‑149 (MQ=255)
                                                                                                                                                 cataaCTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGCCTTTATAGGTCCa  <  2:137625/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCACGCGGTGCAGCAGTTCAATGGCACGGCGTTTGCGCCAGCCAAAACGCTGCCACCAACGGCCTTTATAGGTCCA  >  NC_000913/1352980‑1353274

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