Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,090,229 G→A 6.1% G159S (GGT→AGT)  yggI → Zn‑dependent metalloprotease‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,090,2290GA 6.1% 481.7 / 15.1 147G159S (GGT→AGT) yggIZn‑dependent metalloprotease‑related protein
Reads supporting (aligned to +/- strand):  ref base G (51/87);  new base A (4/5);  total (55/92)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

Only 100 of 154 total aligned reads displayed.

CCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGG  >  NC_000913/3090082‑3090377
                                                                                                                                                   |                                                                                                                                                    
cccgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                     <  2:610100/149‑1 (MQ=255)
cccgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                     <  2:770957/149‑1 (MQ=255)
   gcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTg                                                                                                                                                   >  3:211724/1‑148 (MQ=255)
        tcgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAg                                                                                                                                              <  4:760852/148‑1 (MQ=255)
        tcgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACCGCAGTGAACaca                                                                                                                                             <  1:856228/149‑3 (MQ=255)
        tcgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACCGCAGTGAACac                                                                                                                                              <  4:359764/148‑2 (MQ=255)
         cgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACCGCAGTGAACacag                                                                                                                                            <  1:446046/149‑4 (MQ=255)
               cATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTg                                                                                                                                       <  3:453036/148‑1 (MQ=255)
               cATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTg                                                                                                                                       >  4:477143/1‑148 (MQ=255)
               cATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTg                                                                                                                                       >  3:812063/1‑148 (MQ=255)
                 tCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTTGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGa                                                                                                                                    >  1:280847/1‑149 (MQ=255)
                     ttCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                     <  2:993480/128‑1 (MQ=255)
                     ttCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                     >  1:993480/1‑128 (MQ=255)
                       cGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAAc                                                                                                                              >  1:795121/1‑149 (MQ=255)
                          aCTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCa                                                                                                                            <  4:869347/148‑1 (MQ=255)
                              cAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGc                                                                                                                                             <  3:741655/127‑1 (MQ=255)
                               aaTCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGa                                                                                                                       <  4:825101/148‑1 (MQ=255)
                               aaTCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCAGTGAACAGCTGGTTGCGAAATAACCATCTGaa                                                                                                                      <  2:53272/149‑1 (MQ=255)
                               aaTCCGTGCGTCGCAACACCTTCCCCTACCACTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGa                                                                                                                       <  3:916493/148‑1 (MQ=255)
                                 tCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACt                                                                                                                    >  1:662014/1‑149 (MQ=255)
                                    gTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTCCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATc                                                                                                                 <  2:1055155/149‑1 (MQ=255)
                                      gCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCa                                                                                                                <  4:58842/148‑1 (MQ=255)
                                         tCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGaa                                                                                                            <  1:589801/149‑1 (MQ=255)
                                           gCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACt                                                                                                          <  2:409845/149‑1 (MQ=255)
                                             aaCACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAAt                                                                                                                                 <  1:506081/124‑1 (MQ=255)
                                             aaCACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACttt                                                                                                        <  1:656823/149‑1 (MQ=255)
                                               ccccTCCCCCTACCGCTGCAAGTGCCAGGAGCAGCAGCTTACCGTGCGCTGCCATAAGCGGGGAGTTGGTGGCGAGGAAATCTATCGCTGTATTAACTGCGGTGATCAGATGGTTGCGGAATAACCATCTGAACTATCAGGAACTTTc                                                                                                       <  4:180192/146‑1 (MQ=255)
                                               cacCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTcc                                                                                                      <  1:237110/149‑1 (MQ=255)
                                               cacCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTcc                                                                                                      <  2:598228/149‑1 (MQ=255)
                                                acCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTcc                                                                                                      <  3:541259/148‑1 (MQ=255)
                                                acCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTcc                                                                                                      <  4:877150/148‑1 (MQ=255)
                                                  cTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGa                                                                                                   >  2:866470/1‑149 (MQ=255)
                                                   ttCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGa                                                                                                   >  3:813650/1‑148 (MQ=255)
                                                   ttCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGa                                                                                                   >  3:88276/1‑148 (MQ=255)
                                                   ttCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGAt                                                                                                  <  2:579541/149‑1 (MQ=255)
                                                   ttCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCATACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGa                                                                                                   >  3:391388/1‑148 (MQ=255)
                                                    tCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGAt                                                                                                  <  4:597889/148‑1 (MQ=255)
                                                      cccTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTg                                                                                               <  2:1012466/149‑1 (MQ=255)
                                                       ccTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACttt                                                                                                        <  1:89018/139‑1 (MQ=255)
                                                       ccTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACttt                                                                                                        >  2:89018/1‑139 (MQ=255)
                                                       ccTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTgg                                                                                              <  2:836355/149‑1 (MQ=255)
                                                       ccTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAAAGTTCCTGAtagtt                                                                                              <  1:51065/149‑5 (MQ=255)
                                                            cGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGa                                                                                          <  3:885631/148‑1 (MQ=255)
                                                             gCTGCACTTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTCTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAAATTTCCTGATCTGGCTGAtt                                                                                        <  2:698929/149‑1 (MQ=255)
                                                             gCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGAtt                                                                                        <  1:159803/149‑1 (MQ=255)
                                                             gCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGAtt                                                                                        <  1:953135/149‑1 (MQ=255)
                                                                  aaGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATa                                                                                   <  1:406366/149‑1 (MQ=255)
                                                                        cAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaaaa                                                                             <  2:306625/149‑1 (MQ=255)
                                                                        cAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaaa                                                                              <  3:546645/148‑1 (MQ=255)
                                                                        cAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaaa                                                                              <  4:211724/148‑1 (MQ=255)
                                                                         aGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaaaa                                                                             <  4:545447/148‑1 (MQ=255)
                                                                         aGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaaaa                                                                             >  3:746460/1‑148 (MQ=255)
                                                                          ggAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAAc                                                                            >  4:25446/1‑148 (MQ=255)
                                                                               aTCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCttt                                                                      >  1:381111/1‑149 (MQ=255)
                                                                                tCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTc                                                                     >  1:102646/1‑149 (MQ=255)
                                                                                 cAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACGATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTc                                                                     >  3:665148/1‑148 (MQ=255)
                                                                                   gCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCt                                                                  >  2:272291/1‑149 (MQ=255)
                                                                                              gccgccATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc                                                       >  1:783042/1‑149 (MQ=255)
                                                                                                  ccATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt                                                   <  1:883610/149‑1 (MQ=255)
                                                                                                   cATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGtt                                                  <  2:102646/149‑1 (MQ=255)
                                                                                                   cATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt                                                   <  4:608984/148‑1 (MQ=255)
                                                                                                     tAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAACCAGCTTTCGCTCCGTTGCTGGCTCGttt                                                 >  3:565956/1‑148 (MQ=255)
                                                                                                       aTCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTaa                                              <  2:216744/149‑1 (MQ=255)
                                                                                                        tCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTaa                                              >  4:296951/1‑148 (MQ=255)
                                                                                                             tAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACgg                                         >  3:153162/1‑148 (MQ=255)
                                                                                                               gTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGt                                      <  2:783042/149‑1 (MQ=255)
                                                                                                                 tCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTa                                               <  4:916524/138‑1 (MQ=255)
                                                                                                                 tCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTa                                               >  3:916524/1‑138 (MQ=255)
                                                                                                                 tCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTa                                     <  4:565956/148‑1 (MQ=255)
                                                                                                                     ggCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGGTTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTg                                 <  3:25446/148‑1 (MQ=255)
                                                                                                                     ggCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTg                                 <  4:812063/148‑1 (MQ=255)
                                                                                                                     ggCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTg                                 <  4:111772/148‑1 (MQ=255)
                                                                                                                     ggCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGa                                <  2:568518/149‑1 (MQ=255)
                                                                                                                     ggCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGa                                <  1:949085/149‑1 (MQ=255)
                                                                                                                      gCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGAt                               <  1:554213/149‑1 (MQ=255)
                                                                                                                            ccGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTAcc                          >  4:706336/1‑148 (MQ=255)
                                                                                                                              gTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGt                        >  3:1029950/1‑148 (MQ=255)
                                                                                                                              gTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGt                        >  3:52759/1‑148 (MQ=255)
                                                                                                                                cTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGtt                                                  <  3:398099/120‑1 (MQ=255)
                                                                                                                                  aTCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATATCCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTAtt                    <  4:744747/148‑1 (MQ=255)
                                                                                                                                  aTCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTAttt                   <  1:1035049/149‑1 (MQ=255)
                                                                                                                                  aTCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTAttt                   <  1:343149/149‑1 (MQ=255)
                                                                                                                                  aTCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTAtt                    <  3:557937/148‑1 (MQ=255)
                                                                                                                                  aTCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTAtt                    <  4:664668/148‑1 (MQ=255)
                                                                                                                                     gCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACgg                                         <  4:632640/124‑1 (MQ=255)
                                                                                                                                     gCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACgg                                         >  3:632640/1‑124 (MQ=255)
                                                                                                                                     gCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGt                 >  2:669608/1‑148 (MQ=255)
                                                                                                                                      cTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTc                >  4:741367/1‑148 (MQ=255)
                                                                                                                                      cTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCt               >  2:1070506/1‑149 (MQ=255)
                                                                                                                                      agctgtTCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCt               >  1:221753/4‑149 (MQ=255)
                                                                                                                                         tgtTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTAtt            <  2:889181/149‑1 (MQ=255)
                                                                                                                                         tgtTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTAtt            >  2:387649/1‑149 (MQ=255)
                                                                                                                                         tgtTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTAtt            >  1:590085/1‑149 (MQ=255)
                                                                                                                                         tgtTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTAt             <  4:665148/148‑1 (MQ=255)
                                                                                                                                           ttCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATtgc          >  1:858871/1‑149 (MQ=255)
                                                                                                                                            tCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATtgc          >  4:128509/1‑148 (MQ=255)
                                                                                                                                            tCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATtgc          >  3:378682/1‑148 (MQ=255)
                                                                                                                                                 gCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCgg     <  4:391388/148‑1 (MQ=255)
                                                                                                                                                   ggTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCggtgg  >  2:237452/1‑149 (MQ=255)
                                                                                                                                                   ggTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCggtg   <  4:139744/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGG  >  NC_000913/3090082‑3090377

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