Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,425,289 C→T 6.3% G210S (GGT→AGT)  hisQ ← histidine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,425,2890CT 6.3% 162.6 / 2.7 63G210S (GGT→AGT) hisQhistidine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (35/24);  new base T (2/2);  total (37/26)
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.00e-01

AATAAGCAGCCACAGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATC  >  NC_000913/2425141‑2425437
                                                                                                                                                    |                                                                                                                                                    
aaTAAGCAGCCACAGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACAcc                                                                                                                                                      >  1:1139893/1‑149 (MQ=255)
 aTAAGCAGCCACAGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCa                                                                                                                                                     >  1:328487/1‑149 (MQ=255)
      cagcCACAGAGTGATCGCCACACCTGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGa                                                                                                                                                <  2:151913/149‑1 (MQ=255)
         ccACAGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCCCAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAAc                                                                                                                                             >  1:759252/1‑149 (MQ=255)
         ccACAGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCCCACCAGTGGAAAc                                                                                                                                             >  1:1180149/1‑149 (MQ=255)
         ccACAGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAAc                                                                                                                                             >  1:1183618/1‑149 (MQ=255)
           acaGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACCGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAg                                                                                                                                           >  2:324406/1‑149 (MQ=255)
            caGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCCGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAgt                                                                                                                                          <  2:124091/149‑1 (MQ=255)
             cgattGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGAGCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAgtg                                                                                                                                         <  2:93353/145‑1 (MQ=255)
             agagTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAgtg                                                                                                                                         <  1:908688/149‑1 (MQ=255)
              gagTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAgtgg                                                                                                                                        <  1:990624/149‑1 (MQ=255)
               agTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTCTTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAgtggt                                                                                                                                       <  1:526888/149‑1 (MQ=255)
                 tGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGa                                                                                                                                     >  1:365941/1‑149 (MQ=255)
                   aTCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACTATTGGAAACAGTGGTGaaa                                                                                                                                   >  2:1166501/1‑149 (MQ=255)
                   aTCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACTATTGGAAACAGTGGTGaaa                                                                                                                                   >  2:1161947/1‑149 (MQ=255)
                   aTCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGCAACAGCAGCACACCATTGGAAACAGTGGTGaaa                                                                                                                                   >  2:88091/1‑149 (MQ=255)
                    tCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACACACGGAGTAGCGGCGCTCAAGGAACCGCAGCACACCATTGGCAACAGTGGTGaaaa                                                                                                                                  >  1:872630/1‑149 (MQ=255)
                     cGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAAc                                                                                                                                 <  1:578629/149‑1 (MQ=255)
                       ccACACCAGTAATGCGATAACCGTCGGTCCACAGCAGCGGTATCCAGTCTTCATGTAAGAATTCGATAACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCa                                                                                                                               <  1:1092891/149‑1 (MQ=255)
                         acacCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGACCCGCAGCACACCATTGGAAACAGTGGTGAAAACCAgg                                                                                                                             >  1:383383/1‑149 (MQ=255)
                         acacCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGCAAACAGTGGTGAAAACCAgg                                                                                                                             >  1:267742/1‑149 (MQ=255)
                             cAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACGCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCAGTGGAAACAGTGGTGAAAACCAGGTaaa                                                                                                                         >  1:936458/1‑149 (MQ=255)
                              aGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAAt                                                                                                                        <  1:935047/149‑1 (MQ=255)
                                       aTAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACTCACACGGAGTAGCCGCGCTCAAGGAACAGCCGCACACCATTGGAAACAGTGGTGAAACCCAGGTAAATCACGCcaca                                                                                                               >  2:1157282/1‑149 (MQ=255)
                                       aTAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCCGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTCAAAACCAGGTAAATCACGCcaca                                                                                                               >  1:538337/1‑149 (MQ=255)
                                             gTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGAt                                                                                                         <  2:328487/149‑1 (MQ=255)
                                             gTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGAt                                                                                                         >  2:885241/1‑149 (MQ=255)
                                               cGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATgg                                                                                                       <  1:611652/149‑1 (MQ=255)
                                                 gTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATggag                                                                                                     >  2:431886/1‑147 (MQ=255)
                                                   ccACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACCGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCATCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGaa                                                                                                   >  1:733256/1‑149 (MQ=255)
                                                   ccACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCATCTTCACACCAACGGAGTAGCGGCGCTCAAGGAAAAGAAGCACACCATTGGAAACACTGGTGAACACCAGGTAAATAACGCCACAGCCGATGGCTaa                                                                                                   >  2:2185/1‑149 (MQ=255)
                                                    cacaGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGAAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACATACGATGGCGaaa                                                                                                  >  1:376209/1‑149 (MQ=255)
                                                     acaGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCTCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACAAGCTAAATCACGCCACAGACGATGGCGAAAt                                                                                                 >  2:283919/1‑149 (MQ=255)
                                                              gTTTCCAGTATTCATGTAAGATTTCGATCACAGATCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCCCAGACGATGGCGAAATAGAACGGtt                                                                                        >  1:284822/1‑149 (MQ=255)
                                                               tttCCAGTTTTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCCCACCCACGGAGTAGCGGCGCTCACGGAACAGCCGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTc                                                                                       >  2:966390/1‑149 (MQ=255)
                                                                 tCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGATccc                                                                                     <  1:91026/149‑1 (MQ=255)
                                                                   cAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAg                                                                                   <  2:867298/149‑1 (MQ=255)
                                                                       aTTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACTATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTAc                                                                               <  1:1166501/149‑1 (MQ=255)
                                                                       aTTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACTATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTAc                                                                               <  1:1161947/149‑1 (MQ=255)
                                                                            tGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCCCACCATTGGAAACAGTGGTGAAAACAAGGTAAATCACGCCCCAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTtg                                                                          >  2:1120844/1‑149 (MQ=255)
                                                                              tAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTtgcc                                                                        >  2:948372/1‑149 (MQ=255)
                                                                                      cGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGtt                                                                <  2:1139893/149‑1 (MQ=255)
                                                                                        aTCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTgg                                                              >  1:313498/1‑149 (MQ=255)
                                                                                         tCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCCGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTggg                                                             >  1:917540/1‑149 (MQ=255)
                                                                                               gTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACGAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCtt                                                       <  2:759252/149‑1 (MQ=255)
                                                                                                tCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGGGGTGAACACCAGGTAAATCACGCCACAGACGATGGCGCTATAGCACGGTTCCCAGGGACTTTTGCCTGCCCGTTGGGTGGCgtt                                                      >  1:875800/1‑149 (MQ=255)
                                                                                                 cAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTg                                                     <  2:325964/149‑1 (MQ=255)
                                                                                                   gCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGACAAGCAGCCCACCATTGGAAACAGTGGTGAAAACCAGGTAACTCACGCCACAGACGATGGCGCAATAGAACGGTTCCCAGGTACTTTTTCCTGCCAGTTGGGTGGCTTTGAc                                                   >  1:1072895/1‑149 (MQ=255)
                                                                                                      ctctTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACcac                                                >  1:839616/1‑149 (MQ=255)
                                                                                                      cacatCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACcac                                                <  2:365941/145‑1 (MQ=255)
                                                                                                              acCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACTCCATTGGAAAGAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCa                                        <  1:1657/149‑1 (MQ=255)
                                                                                                              acCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGTTTTGACCACATCTTCCa                                        <  1:885241/149‑1 (MQ=255)
                                                                                                              acCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCa                                        >  1:586345/1‑149 (MQ=255)
                                                                                                                  aCGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGcc                                    <  2:943505/149‑1 (MQ=255)
                                                                                                                    ggAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGa                                  <  1:324406/149‑1 (MQ=255)
                                                                                                                        ttgctgcgcTCAAGGAACAGCACCACACCAATGGAAACAGTGGTGAAAACCATGTAAATCACCCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTaa                              <  2:1072895/144‑1 (MQ=255)
                                                                                                                         aGCGGCGCTCCAGGAACAGCAGCTCACCTTTGGAAACATTGGTGAAAACAAGGTAAATCACGCCACAGACGATGGCGAACTAGAACGGTTCCCAGGTACTTTTGCCTGCGAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAAc                             <  2:538337/149‑1 (MQ=255)
                                                                                                                                  cAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGTTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGTCTTTGACCACATCATCCAGTCCGAGTAACGACACCTGt                    >  2:990190/1‑149 (MQ=255)
                                                                                                                                   aaGGAACAGCAGCACACCATTGGAAACAGTGGTGAATACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTg                   <  2:313498/149‑1 (MQ=255)
                                                                                                                                   aaGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGCACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGCCCACATCTTCCAGGCCGAGTACCGAACCCAGTg                   >  1:1086768/1‑149 (MQ=255)
                                                                                                                                    aGGAACAGCAGCACACCATTGGACACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGGACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGCGTACCGCAACCAGTGc                  >  1:808681/1‑149 (MQ=255)
                                                                                                                                        aCAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCATTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGTa              >  2:237796/1‑149 (MQ=255)
                                                                                                                                               cacaCCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGa       >  1:120307/1‑149 (MQ=255)
                                                                                                                                                 cacCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAgg     >  2:722900/1‑149 (MQ=255)
                                                                                                                                                   ccATTGGAACCAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCGCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGAGCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGGAGATTTGAGGAt   >  2:1165964/1‑149 (MQ=255)
                                                                                                                                                    cATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATc  >  2:900186/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AATAAGCAGCCACAGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACAGTGGTGAAAACCAGGTAAATCACGCCACAGACGATGGCGAAATAGAACGGTTCCCAGGTACTTTTGCCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATC  >  NC_000913/2425141‑2425437

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: ‑