Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,651,365 T→C Q308R (CAA→CGA)  yfhM ← bacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,651,3650TC98.7% 272.2 / ‑4.2 75Q308R (CAA→CGA) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (35/39);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTTGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCCTGCTGAAGCGGTT  >  NC_000913/2651218‑2651509
                                                                                                                                                   |                                                                                                                                                
ttACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCgg                                                                                                                                                <  2:1153699/150‑1 (MQ=255)
  aCGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGt                                                                                                                                              <  2:49959/150‑1 (MQ=255)
    gcgcCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTaaa                                                                                                                                           <  1:383413/151‑1 (MQ=255)
     cgcCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAg                                                                                                                                          <  2:282470/151‑1 (MQ=255)
      gcCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTTGGGTAAAGa                                                                                                                                         <  2:7574/151‑1 (MQ=255)
      gcCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTaa                                                                                                                                            >  2:766106/1‑148 (MQ=255)
        cAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGat                                                                                                                                        >  2:845960/1‑150 (MQ=255)
        cAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGat                                                                                                                                        >  2:943747/1‑150 (MQ=255)
         aaCAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGatat                                                                                                                                      <  1:1128654/151‑1 (MQ=255)
          aCAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATc                                                                                                                                     <  2:1098663/151‑1 (MQ=255)
          aCAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATc                                                                                                                                     >  1:1094749/1‑151 (MQ=255)
            aGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATcc                                                                                                                                    >  2:285109/1‑150 (MQ=255)
               aaTGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGa                                                                                                                                 <  2:1037592/150‑1 (MQ=255)
                       gTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGa                                                                                                                        <  1:379618/151‑1 (MQ=255)
                          tttATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATa                                                                                                                      <  2:871422/150‑1 (MQ=255)
                            tATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACg                                                                                                                   >  2:20687/1‑151 (MQ=255)
                             aTCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACgg                                                                                                                  >  2:431936/1‑151 (MQ=255)
                               cATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTg                                                                                                                >  1:28741/1‑151 (MQ=255)
                                  tttCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGc                                                                                                             >  2:1241311/1‑151 (MQ=255)
                                   ttCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCt                                                                                                            >  2:326097/1‑151 (MQ=255)
                                       aGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGaa                                                                                                         >  1:1104438/1‑150 (MQ=255)
                                        gCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAAc                                                                                                       <  1:600192/151‑1 (MQ=255)
                                         cTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACg                                                                                                      >  1:499732/1‑151 (MQ=255)
                                             acaTGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCg                                                                                                   >  1:1236839/1‑150 (MQ=255)
                                                    ccTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACt                                                                                           <  2:60573/151‑1 (MQ=255)
                                                       gAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACtt                                                                                          <  1:845922/149‑1 (MQ=255)
                                                       gAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTa                                                                                         <  2:827065/150‑1 (MQ=255)
                                                       gAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTa                                                                                         <  2:425584/150‑1 (MQ=255)
                                                       gAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTa                                                                                         >  2:1141257/1‑150 (MQ=255)
                                                            tCACTGGTTGCCTGAGTCAGAGTCTGCTCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTa                                                                                    <  2:785387/150‑1 (MQ=255)
                                                             cACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTaa                                                                                   <  2:644092/150‑1 (MQ=255)
                                                                 ggTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTACGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACa                                                                                <  1:766070/149‑1 (MQ=255)
                                                                  gTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCg                                                                             <  2:388281/151‑1 (MQ=255)
                                                                     gCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTcg                                                                          <  2:1145006/151‑1 (MQ=255)
                                                                         gAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCggg                                                                      <  1:374008/151‑1 (MQ=255)
                                                                          aGTCAGAGTCTGCCCTTTCTCATTTAATAACGAGACTTCAATTCCTTGCTGGGACGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCggg                                                                      >  1:340460/1‑150 (MQ=255)
                                                                             cAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATAGCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGAtt                                                                   <  1:662665/150‑1 (MQ=255)
                                                                               gagTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACt                                                                >  2:544804/1‑151 (MQ=255)
                                                                                     gCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAAt                                                           <  1:719216/150‑1 (MQ=255)
                                                                                     gCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATc                                                          <  1:105672/151‑1 (MQ=255)
                                                                                       ccTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGt                                                        <  1:1166328/151‑1 (MQ=255)
                                                                                          ttCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTaa                                                      <  2:1122934/150‑1 (MQ=255)
                                                                                              cATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGt                                                   <  2:905183/149‑1 (MQ=255)
                                                                                              cATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGt                                                   <  1:824145/149‑1 (MQ=255)
                                                                                              cATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTcc                                                 >  2:1022844/1‑151 (MQ=255)
                                                                                              cATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTc                                                  <  2:758360/150‑1 (MQ=255)
                                                                                                tttAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTcc                                                 >  2:275567/1‑149 (MQ=255)
                                                                                                tttAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTcc                                                 >  2:64221/1‑149 (MQ=255)
                                                                                                  taataaAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTCACGTcc                                                 >  1:980532/1‑147 (MQ=255)
                                                                                                   aataaAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGcc                                             <  1:273675/150‑1 (MQ=255)
                                                                                                    ataaAGCGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTg                                           >  1:912058/1‑151 (MQ=255)
                                                                                                        agagACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATtc                                       >  2:951468/1‑151 (MQ=255)
                                                                                                          agaCTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATtcat                                     >  2:395446/1‑151 (MQ=255)
                                                                                                               tCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATcaca                                 <  2:705922/150‑1 (MQ=255)
                                                                                                               tCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATcac                                  >  2:1091577/1‑149 (MQ=255)
                                                                                                                cAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGc                               >  1:690429/1‑151 (MQ=255)
                                                                                                                 aaTTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGcc                              <  2:306326/151‑1 (MQ=255)
                                                                                                                  aTTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCa                             <  2:373825/151‑1 (MQ=255)
                                                                                                                     ccTTGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGa                           <  2:940575/150‑1 (MQ=255)
                                                                                                                       ttGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATac                        >  2:288486/1‑151 (MQ=255)
                                                                                                                       ttGCTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATac                        >  2:733740/1‑151 (MQ=255)
                                                                                                                          cTGGGCCGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACAcg                     >  1:550190/1‑151 (MQ=255)
                                                                                                                               ccGCGCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAAGGTTGTGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACAcgcccgc                <  1:1009075/151‑1 (MQ=255)
                                                                                                                                cgcgCCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACAcgcccgcc               <  2:348366/151‑1 (MQ=255)
                                                                                                                                 gcgcCGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCct              >  2:434162/1‑151 (MQ=255)
                                                                                                                                    ccGTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCctgct           >  1:248092/1‑151 (MQ=255)
                                                                                                                                      gTTTTCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCctgctg          >  1:1206810/1‑150 (MQ=255)
                                                                                                                                       ttttCCAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGCATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCctgct           >  2:1253449/1‑148 (MQ=255)
                                                                                                                                           ccAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCCTGCTGAAGCgg    <  1:1022806/151‑1 (MQ=255)
                                                                                                                                           ccAGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCCTGCTGAAGCg     <  2:624240/150‑1 (MQ=255)
                                                                                                                                             aGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCCTGCTGAAGCGGtt  >  1:545520/1‑151 (MQ=255)
                                                                                                                                             aGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCCTGCTGAAGCGGtt  >  1:187543/1‑151 (MQ=255)
                                                                                                                                             aGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCCTGCTGAAGCGGtt  >  1:1268300/1‑151 (MQ=255)
                                                                                                                                             aGACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCCTGCTGAAGCGGtt  <  1:1141221/151‑1 (MQ=255)
                                                                                                                                              gACTTCGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCCTGCTGAAGCGGtt  <  2:968708/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTTTTCCAGACTTTGGGTAAAGATATCCAGACGATTGTGATAACGGTGAGCTGAAACGCCGATATCACTTAACGTAAACAGCGTCGCGGGATTACTGTAATCGTAACGTCCAGCCTGATTCATCACAGCCAGATACACGCCCGCCTGCTGAAGCGGTT  >  NC_000913/2651218‑2651509

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