Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,501,225 C→T 7.7% intergenic (‑434/‑35) yjgX ← / → yjgZ pseudogene fragment/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,501,2250CT 7.7% 239.5 / 6.4 76intergenic (‑434/‑35)yjgX/yjgZpseudogene fragment/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (29/41);  new base T (2/4);  total (31/45)
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.65e-01

CCCCGGTGGGTATGACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATC  >  NC_000913/4501077‑4501368
                                                                                                                                                    |                                                                                                                                               
ccccGGTGGGTATGACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTcc                                                                                                                                                 >  4:275242/1‑149 (MQ=255)
 cccGGTGGGTATGACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCt                                                                                                                                                <  4:293741/149‑1 (MQ=255)
 cccGGTGGGTATGACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCt                                                                                                                                                <  3:275242/149‑1 (MQ=255)
    ggTGGGTATGACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTgaa                                                                                                                                             <  4:110544/149‑1 (MQ=255)
          tATGACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCTTCAGGAAGTc                                                                                                                                       <  2:261130/149‑1 (MQ=255)
          tATGACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCTTCAGGAAGTc                                                                                                                                       <  1:83651/149‑1 (MQ=255)
            tGACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGt                                                                                                                                     >  3:351372/1‑149 (MQ=255)
             gACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCTTCAGGAAGTcaac                                                                                                                                    <  4:185142/149‑4 (MQ=255)
             gACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGtt                                                                                                                                    <  3:205590/149‑1 (MQ=255)
                 gATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAg                                                                                                                                <  3:443759/149‑1 (MQ=255)
                  aTGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAgg                                                                                                                               <  2:429576/149‑1 (MQ=255)
                  aTGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAgg                                                                                                                               <  4:351372/149‑1 (MQ=255)
                     ccTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGgag                                                                                                                            >  3:183905/1‑149 (MQ=255)
                            cACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCa                                                                                                                     <  4:456882/149‑1 (MQ=255)
                             aCTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGTGAGATTCAg                                                                                                                    >  2:224825/1‑149 (MQ=255)
                              ctctctTCCCGGTCATCGTCGGCCTTCACCGTTATTGTGGCGTGTCAGCACCCTCGAAGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTAAGc                                                                                                                   >  2:381161/1‑149 (MQ=255)
                               tctctTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGcc                                                                                                                  >  2:386802/1‑149 (MQ=255)
                                 tctTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAg                                                                                                                                <  4:144404/133‑1 (MQ=255)
                                 tctTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAg                                                                                                                                >  3:144404/1‑133 (MQ=255)
                                   ttCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCttt                                                                                                              <  1:93718/149‑1 (MQ=255)
                                    tCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCtttt                                                                                                             <  3:248336/149‑1 (MQ=255)
                                       cGGTCATCGTCTGCCTTCCCCTGTATTGTGGCGTGTCAGCACACTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGtt                                                                                                          <  4:142084/149‑1 (MQ=255)
                                          tCATCGTCGGCCGTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTAcc                                                                                                       <  2:379569/149‑1 (MQ=255)
                                             tcgtcgGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTAcctcc                                                                                                    <  4:151475/149‑1 (MQ=255)
                                                      ttCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGct                                                                                           >  1:319691/1‑149 (MQ=255)
                                                        ccccGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGcttc                                                                                         <  1:386802/149‑1 (MQ=255)
                                                        ccccGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGcttc                                                                                         <  4:22175/149‑1 (MQ=255)
                                                         cccGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGcttct                                                                                        >  3:134071/1‑149 (MQ=255)
                                                          ccGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGgagaga                                                                                                                         >  2:406183/1‑115 (MQ=255)
                                                          ccGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGgagaga                                                                                                                         <  1:406183/115‑1 (MQ=255)
                                                              tATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTcc                                                                                   <  2:64984/149‑1 (MQ=255)
                                                                 tgtgGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTc                                                                                <  3:143373/149‑1 (MQ=255)
                                                                  gtgGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCAGAACCCGTTGACTTCCTGAGAAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTACTTcc                                                                               <  4:302612/149‑1 (MQ=255)
                                                                       gtgtCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGt                                                                                                                                     <  1:312270/90‑1 (MQ=255)
                                                                       gtgtCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGt                                                                                                                                     >  2:312270/1‑90 (MQ=255)
                                                                       gtgtCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCg                                                                          >  3:188198/1‑149 (MQ=255)
                                                                             gCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCga                                                                    <  4:188198/149‑1 (MQ=255)
                                                                               aCCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCgaga                                                                  <  3:252113/149‑1 (MQ=255)
                                                                                 ccTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAg                                                                                                                    >  2:331738/1‑97 (MQ=255)
                                                                                 ccTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAg                                                                                                                    <  1:331738/97‑1 (MQ=255)
                                                                                       cGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTgaagaa                                                                                                                                          >  4:277854/1‑69 (MQ=255)
                                                                                       cGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTgaagaa                                                                                                                                          <  3:277854/69‑1 (MQ=255)
                                                                                       cGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCg                                                          <  3:116937/149‑1 (MQ=255)
                                                                                          cATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCgggg                                                       <  3:339222/149‑1 (MQ=255)
                                                                                           aTTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGAtt                                                                                                                       <  2:176427/84‑1 (MQ=255)
                                                                                           aTTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGAtt                                                                                                                       >  1:176427/1‑84 (MQ=255)
                                                                                           aTTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGc                                                      <  2:268988/149‑1 (MQ=255)
                                                                                           aTTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGc                                                      >  1:343966/1‑149 (MQ=255)
                                                                                               cGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGAc                                                  <  4:171554/149‑1 (MQ=255)
                                                                                                      tgctCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCAcc                                           >  3:215675/1‑149 (MQ=255)
                                                                                                      tgctCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCAcc                                           >  3:348731/1‑149 (MQ=255)
                                                                                                           aGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACgg                                      >  2:54086/1‑149 (MQ=255)
                                                                                                              gtTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGaaa                                   >  1:454559/1‑149 (MQ=255)
                                                                                                              gtTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGATTCAGGGACAGCGTCAGCCCTCTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGGTAGTTAGTCGGGGCGGACTTTCACTGACGGaaa                                   >  4:56053/1‑149 (MQ=255)
                                                                                                               ttCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAg                                  <  4:215675/149‑1 (MQ=255)
                                                                                                               ttCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTACTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAg                                  >  3:399434/1‑149 (MQ=255)
                                                                                                                 cactcacGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGAt                                <  2:319691/149‑1 (MQ=255)
                                                                                                                  actcacGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATc                               >  3:182191/1‑149 (MQ=255)
                                                                                                                    tcacGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGct                                                                                           >  1:233756/1‑87 (MQ=255)
                                                                                                                    tcacGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGt                             <  1:54086/149‑1 (MQ=255)
                                                                                                                    tcacGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTGACCTCCTGGACCGct                                                                                           <  2:233756/87‑1 (MQ=255)
                                                                                                                          tACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCAttt                       <  3:39292/149‑1 (MQ=255)
                                                                                                                              gCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAg                   <  2:454559/149‑1 (MQ=255)
                                                                                                                              gCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAg                   <  4:20181/149‑1 (MQ=255)
                                                                                                                              gCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAg                   <  1:359973/149‑1 (MQ=255)
                                                                                                                               cccGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAgg                  >  4:423668/1‑149 (MQ=255)
                                                                                                                                     tAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGcacaca            >  4:180042/1‑149 (MQ=255)
                                                                                                                                        cccGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTc                                              >  3:458353/1‑112 (MQ=255)
                                                                                                                                        cccGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTc                                              <  4:458353/112‑1 (MQ=255)
                                                                                                                                           aacgACTTCTTCAGGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACAcaacaa      >  3:383957/4‑149 (MQ=255)
                                                                                                                                              gACTTCTTCAGGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCg                                                     <  4:297532/99‑1 (MQ=255)
                                                                                                                                              gACTTCTTCAGGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCg                                                     >  3:297532/1‑99 (MQ=255)
                                                                                                                                               aCTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACacaa         >  2:286204/1‑142 (MQ=255)
                                                                                                                                               aCTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACacaa         <  1:286204/142‑1 (MQ=255)
                                                                                                                                               aCTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATc  >  4:382096/1‑149 (MQ=255)
                                                                                                                                                 ttCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTcc                                                                                   >  2:424519/1‑66 (MQ=255)
                                                                                                                                                 ttCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTcc                                                                                   <  1:424519/66‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCCCGGTGGGTATGACCGATGCCTTTCTCACTCTCTTCCCGGTCATCGTCGGCCTTCCCCGTTATTGTGGCGTGTCAGCACCCTCGACGGCATTTCGATGGCTGCTCAGTGTTCACTCACGTTACGGCCCGCATAACCCGTTGACTTCCTGAAGAAGTCGTTTCAGGGAGAGATTCAGCCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATC  >  NC_000913/4501077‑4501368

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