Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,649,813 T→A 5.5% L825F (TTA→TTT 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,649,8130TA 5.5% 219.6 / 3.2 73L825F (TTA→TTTyfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (32/37);  new base A (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

TGGGACGCATTTTTGCTGGTGTTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGAT  >  NC_000913/2649674‑2649959
                                                                                                                                           |                                                                                                                                                  
tGGGACGCTTTTTTGCTGGTGTTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCTGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACa                                                                                                                                             <  4:270225/147‑1 (MQ=255)
  ggACGCATTTTTGCTGGTGTTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGAt                                                                                                                                          <  7:330994/148‑1 (MQ=255)
      gCATTTTTGCTGGTGTTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATg                                                                                                                                      <  4:170909/148‑1 (MQ=255)
       cATTTTTGCTGGTGTTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGa                                                                                                                                     >  5:197006/1‑148 (MQ=255)
          ttttGCTGGTGTTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGa                                                                                                                                  >  1:241769/1‑148 (MQ=255)
                 ggTGTTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAgg                                                                                                                           >  5:227188/1‑148 (MQ=255)
                   tgtTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGtt                                                                                                                         >  5:373255/1‑148 (MQ=255)
                    gtTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTAAAATACAGATCATGACGATTCCAGGttt                                                                                                                        >  5:124392/1‑148 (MQ=255)
                    gtTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTAAAATACAGATCATGACGATTCCAGGttt                                                                                                                        >  7:262247/1‑148 (MQ=255)
                      ttCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGGGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACa                                                                                                                                             >  1:148257/1‑125 (MQ=255)
                      ttCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACa                                                                                                                                             <  2:148257/125‑1 (MQ=255)
                        ccAGCGCCAGATCGAGGCGACGGTGTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCACGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTAt                                                                                                                     <  4:281447/147‑1 (MQ=255)
                                   tCGAGGCGACGATTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCg                                                                                                         <  2:158011/148‑1 (MQ=255)
                                    cGAGGCGACGGTTTTCATCGCCAATCGGCAGATGCAACACACCTACCGCTCGTTTTGGCTTCGCGGAGCGAGATTTATCGCATTGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCg                                                                                                         <  2:289414/147‑1 (MQ=255)
                                     gAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGt                                                                                                        <  4:8805/147‑1 (MQ=255)
                                        gCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGAACCACCAGCGTACTTAAATACAGATCATGACGAATCCAGGTTTTATCGACCGGAATCGTCAGa                                                                                                    <  3:88083/148‑1 (MQ=255)
                                         cGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGAt                                                                                                   <  7:48450/148‑1 (MQ=255)
                                           aCGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATc                                                                                                  >  8:286025/1‑147 (MQ=255)
                                                    aTCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGa                                                                                        <  6:375956/148‑1 (MQ=255)
                                                         cAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTGTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGccc                                                                                    >  1:261767/1‑147 (MQ=255)
                                                          aaGCGGCAGTTGGAACACACCAAGCGCGTGTTTTGGCGTCGCGGAGCGAGTTTTATCGTCAGGACGTTCCACCTGCGTACTTAAATACTGATAATGAAGATTCTAGGTTTTATCGAGGGGATTCGTGAGATCGAGCCGTTGAGCGCGa                                                                                  <  6:197006/148‑1 (MQ=255)
                                                          aaGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGa                                                                                  <  8:15168/148‑1 (MQ=255)
                                                           aGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATAc                                                                                                                                              <  6:369365/87‑1 (MQ=255)
                                                           aGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATAc                                                                                                                                              >  5:369365/1‑87 (MQ=255)
                                                            gCGGGAGATGAAACACACCAAGCGCTCGTTTTGGCGTCGCGGAGCGAGATTCATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAAc                                                                                <  6:375669/148‑1 (MQ=255)
                                                                cAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGt                                                                                                        <  2:309630/120‑1 (MQ=255)
                                                                cAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGt                                                                                                        >  1:309630/1‑120 (MQ=255)
                                                                 aGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCaa                                                                           <  6:350539/148‑1 (MQ=255)
                                                                  gATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAAt                                                                          >  6:155147/1‑148 (MQ=255)
                                                                  gATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAAt                                                                          >  8:66090/1‑148 (MQ=255)
                                                                   aTGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAAtc                                                                         >  3:257833/1‑148 (MQ=255)
                                                                   aTGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAAt                                                                          >  1:250451/1‑147 (MQ=255)
                                                                    tGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGcc                                                                                     >  5:201132/1‑135 (MQ=255)
                                                                    tGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGcc                                                                                     <  6:201132/135‑1 (MQ=255)
                                                                      cAACACACCAACAGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATctct                                                                      >  1:61831/1‑148 (MQ=255)
                                                                        acacacCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGa                                                                                                    <  4:151647/116‑1 (MQ=255)
                                                                        acacacCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGa                                                                                                    >  3:151647/1‑116 (MQ=255)
                                                                        acacacCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTg                                                                    >  2:91645/1‑148 (MQ=255)
                                                                        acacacCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTg                                                                    >  7:240865/1‑148 (MQ=255)
                                                                         cacacCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGc                                                                   >  6:274904/1‑148 (MQ=255)
                                                                          ccacCAACCCCCCCGTCTGGCCCCGCGGAGCGAGACTCACCGCCACCCCGCACCACCAGCGCACCTCAATCCAGATCACGACGATTCCAGGTTTCATCGACCGGAATCGTCAGATCCAGCCCGTGAGCCCCAACCTCACCCGCTTGcc                                                                  <  4:293746/147‑1 (MQ=255)
                                                                          acacCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTAAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGcc                                                                  <  8:262247/148‑1 (MQ=255)
                                                                          acacCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTAAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGcc                                                                  <  6:124392/148‑1 (MQ=255)
                                                                          acacCAACCGCGCGGTTTGGCGTCGCGGAGCGGGATTTCTCGCCAGGACGTACCACCAGCGTACTTAAATTCAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGcc                                                                  <  4:276856/148‑1 (MQ=255)
                                                                            acCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGccac                                                                >  1:348905/1‑148 (MQ=255)
                                                                               aaCCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGccacca                                                              >  6:248356/1‑147 (MQ=255)
                                                                                   gcgcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGa                                                                                                             <  8:326749/96‑1 (MQ=255)
                                                                                   gcgcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGa                                                                                                             >  7:326749/1‑96 (MQ=255)
                                                                                   gcgcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCg                                                         >  4:284611/1‑148 (MQ=255)
                                                                                    cgcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATcc                                                                                                 <  4:352430/107‑1 (MQ=255)
                                                                                    cgcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATcc                                                                                                 >  3:352430/1‑107 (MQ=255)
                                                                                      cgTGTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGcc                                                      <  6:140775/148‑1 (MQ=255)
                                                                                      cgTGTTCGCCTCGCGGAGCGAGACTTATCGCCAGCACGTACCACCACCGTACTTAAAAACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGcc                                                      <  8:96148/148‑1 (MQ=255)
                                                                                        ttttGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCt                                                    <  3:119227/148‑1 (MQ=255)
                                                                                          ttGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCtt                                                   >  1:84927/1‑147 (MQ=255)
                                                                                           tGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTc                                                  <  2:348905/147‑1 (MQ=255)
                                                                                               gTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTgg                                             <  2:241769/148‑1 (MQ=255)
                                                                                               gTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTgg                                             <  7:286025/148‑1 (MQ=255)
                                                                                                  gcgGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTCTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACt                                          <  6:227188/148‑1 (MQ=255)
                                                                                                         gagaTTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCAt                                   >  4:276879/1‑148 (MQ=255)
                                                                                                            aTTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGc                                <  5:155147/148‑1 (MQ=255)
                                                                                                              ttATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGCCCATCGCa                               >  5:10648/1‑147 (MQ=255)
                                                                                                                 tCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATaa                            <  5:274904/147‑1 (MQ=255)
                                                                                                                 tCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAAc                           <  2:250451/148‑1 (MQ=255)
                                                                                                                   gCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAAcc                          <  4:169557/147‑1 (MQ=255)
                                                                                                                      aGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACAt                                                                              >  3:196516/1‑92 (MQ=255)
                                                                                                                      aGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACAt                                                                              <  4:196516/92‑1 (MQ=255)
                                                                                                                              ccaccaGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGt              >  3:355056/1‑148 (MQ=255)
                                                                                                                               caccaGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGtt             >  2:189678/1‑148 (MQ=255)
                                                                                                                               caccaGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGtt             <  1:91645/148‑1 (MQ=255)
                                                                                                                                accaGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCGGCCCTTGAGCCCGAACCTCAATCGCTTGCCACCACAGCGGCCCTTCACTGGCCTCGACCATCGCATAACCTTTACCCGCCGTTg            >  4:247552/1‑148 (MQ=255)
                                                                                                                                  caGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGgc          >  7:254016/1‑148 (MQ=255)
                                                                                                                                         cTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCgcggcg    <  3:284611/147‑1 (MQ=255)
                                                                                                                                          ttAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGat  <  8:240865/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
TGGGACGCATTTTTGCTGGTGTTTCCAGCGCCAGATCGAGGCGACGGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGAT  >  NC_000913/2649674‑2649959

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