Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,274,744 T→C 6.4% intergenic (+83/+16) ssb → / ← yjcB single‑stranded DNA‑binding protein/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,274,7440TC 6.4% 551.6 / 18.0 169intergenic (+83/+16)ssb/yjcBsingle‑stranded DNA‑binding protein/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (83/75);  new base C (6/5);  total (89/80)
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.18e-01

Only 100 of 170 total aligned reads displayed.

CAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGC  >  NC_000913/4274596‑4274891
                                                                                                                                                    |                                                                                                                                                   
cagTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAt                                                                                                                                                     <  1:81487/149‑1 (MQ=255)
    ccGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAg                                                                                                                                                 <  3:351408/149‑1 (MQ=255)
        tCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCtt                                                                                                                                             <  4:500569/149‑1 (MQ=255)
         ccGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTc                                                                                                                                            <  4:1148085/149‑1 (MQ=255)
         ccGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCTTc                                                                                                                                            <  4:655117/149‑1 (MQ=255)
           ggCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCTTCtt                                                                                                                                          <  1:658922/149‑2 (MQ=255)
            gCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTAc                                                                                                                                         >  3:139224/1‑149 (MQ=255)
            gCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCTTCttt                                                                                                                                         <  2:246877/149‑3 (MQ=255)
              aGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGt                                                                                                                                       <  4:616972/149‑1 (MQ=255)
                       tAACGAGCCGCGGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCa                                                                                                                              >  3:1145505/1‑149 (MQ=255)
                        aaCGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCaa                                                                                                                             >  1:922419/1‑149 (MQ=255)
                             gccgccGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCa                                                                                                                        >  4:707619/1‑149 (MQ=255)
                             gccgccGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCa                                                                                                                        <  1:791859/149‑1 (MQ=255)
                               cgccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCAtt                                                                                                                      <  2:608896/149‑1 (MQ=255)
                               cgccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCAtt                                                                                                                      >  4:745964/1‑149 (MQ=255)
                                 ccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTgg                                                                                                                    >  2:40166/1‑149 (MQ=255)
                                       gACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATtg                                                                                                              >  3:549683/1‑149 (MQ=255)
                                          tttGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTgg                                                                                                           >  3:685343/1‑149 (MQ=255)
                                          tttGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTgg                                                                                                           <  4:1076205/149‑1 (MQ=255)
                                                gatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCaaaa                                                                                                                           >  3:443330/1‑127 (MQ=255)
                                                gatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCaaaa                                                                                                                           <  4:443330/127‑1 (MQ=255)
                                                   gaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCAt                                                                                                                       <  1:1104101/128‑1 (MQ=255)
                                                   gaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCAt                                                                                                                       >  2:1104101/1‑128 (MQ=255)
                                                   gaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTg                                                                                                  >  1:209299/1‑149 (MQ=255)
                                                     cATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAt                                                                                                >  4:1211426/1‑149 (MQ=255)
                                                     cATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAt                                                                                                <  4:629534/149‑1 (MQ=255)
                                                      aTTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAtt                                                                                               <  1:310068/149‑1 (MQ=255)
                                                       ttCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTc                                                                                              <  3:1351398/149‑1 (MQ=255)
                                                            ttCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtctc                                                                                         >  2:144795/1‑149 (MQ=255)
                                                               tGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttc                                                                                      <  2:1157321/149‑1 (MQ=255)
                                                               tGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttc                                                                                      >  2:703628/1‑149 (MQ=255)
                                                                gATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttct                                                                                     <  3:1414668/149‑1 (MQ=255)
                                                                  tttGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCGTCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTat                                                                                   >  4:188377/1‑149 (MQ=255)
                                                                  tttGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTat                                                                                   >  4:460460/1‑149 (MQ=255)
                                                                      tCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTatcatc                                                                               >  2:66767/1‑149 (MQ=255)
                                                                      tCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTatcatc                                                                               <  1:1143366/149‑1 (MQ=255)
                                                                      tCATTAAAACAATAAGCTCTATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTatcatc                                                                               >  4:408331/1‑149 (MQ=255)
                                                                       cATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTAtcatca                                                                              <  2:797192/149‑1 (MQ=255)
                                                                        aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                             >  3:921671/1‑149 (MQ=255)
                                                                        aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                             >  2:369600/1‑149 (MQ=255)
                                                                        aTTAAAACAATAGGTTATATTGTTTTAAGGAGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTAATGGGGCGTTCTTGGTATCGGGTCAAAATCTTTGGCCTTTGTGGCGGTTTTTGATTCATTTCTTCTATCATCAc                                                                             >  3:1113360/1‑149 (MQ=255)
                                                                         ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtctc                                                                                         >  1:612799/1‑136 (MQ=255)
                                                                         ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAcc                                                                            >  1:581251/1‑149 (MQ=255)
                                                                          tAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCg                                                                           >  4:389787/1‑149 (MQ=255)
                                                                                aaTAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGa                                                                     <  2:922419/149‑1 (MQ=255)
                                                                                  tAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtg                                                                   >  1:1321894/1‑149 (MQ=255)
                                                                                    ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgtt                                                                 <  3:968570/149‑1 (MQ=255)
                                                                                    ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgtt                                                                 >  3:1292028/1‑149 (MQ=255)
                                                                                     gTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAtt                                                                                               >  2:316597/1‑118 (MQ=255)
                                                                                      ttATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGTGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttgt                                                               >  4:840456/1‑149 (MQ=255)
                                                                                           ttGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                                 >  4:547890/1‑110 (MQ=255)
                                                                                           ttGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                                 <  3:547890/110‑1 (MQ=255)
                                                                                              ttttAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATa                                                       >  3:315384/1‑149 (MQ=255)
                                                                                                 tAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCg                                                    <  1:1011623/149‑1 (MQ=255)
                                                                                                  aaGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCgg                                                   <  1:270523/149‑1 (MQ=255)
                                                                                                     gTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGa                                                <  4:1292028/149‑1 (MQ=255)
                                                                                                     gTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGa                                                <  3:707619/149‑1 (MQ=255)
                                                                                                        gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGt                                             >  3:817169/1‑149 (MQ=255)
                                                                                                         atgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTc                                            >  4:1383272/1‑149 (MQ=255)
                                                                                                           gatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCg                                          >  2:1365658/1‑149 (MQ=255)
                                                                                                             ttAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGcc                                        >  1:997672/1‑149 (MQ=255)
                                                                                                              tAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGc                                         >  3:1220060/1‑147 (MQ=255)
                                                                                                              tAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGc                                         <  4:1220060/147‑1 (MQ=255)
                                                                                                              tAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCa                                       <  3:52567/149‑1 (MQ=255)
                                                                                                              tAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCa                                       >  1:6824/1‑149 (MQ=255)
                                                                                                               aaaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGACGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGAt                                               <  1:3305/140‑1 (MQ=255)
                                                                                                                  gCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATtgt                                   <  3:620942/149‑1 (MQ=255)
                                                                                                                     tCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATg                                              >  3:1037968/1‑135 (MQ=255)
                                                                                                                     tCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                                >  1:1449111/1‑149 (MQ=255)
                                                                                                                     tCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                                <  3:745964/149‑1 (MQ=255)
                                                                                                                       tGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGc                              <  1:1263172/149‑1 (MQ=255)
                                                                                                                        gccagccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGcc                             <  1:369600/149‑1 (MQ=255)
                                                                                                                          cagccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCaa                           <  4:921671/149‑1 (MQ=255)
                                                                                                                          cagccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCaa                           <  1:1215275/149‑1 (MQ=255)
                                                                                                                              caTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCaa                       <  3:285921/149‑1 (MQ=255)
                                                                                                                               aTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGc                                         >  1:318019/1‑130 (MQ=255)
                                                                                                                               aTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGc                                         <  2:318019/130‑1 (MQ=255)
                                                                                                                               aTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAc                      <  3:262762/149‑1 (MQ=255)
                                                                                                                               aTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAc                      <  1:474742/149‑1 (MQ=255)
                                                                                                                                 tacgtagaagCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGa                                                                     >  2:478305/6‑100 (MQ=38)
                                                                                                                                 tacgtagaagCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGa                                                                     <  1:478305/95‑1 (MQ=38)
                                                                                                                                  aaaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAcccc                   >  3:396996/1‑149 (MQ=255)
                                                                                                                                   aaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCa                  >  1:459309/1‑149 (MQ=255)
                                                                                                                                     aagaagCCTACTTTATGGAGACTTCTACATATCAAGTCAAAATCATTGGCCATTGTGCCGGTGTCTGATTCGTCTCTTCTATCATCACCGGCTCGATGTTGTGCCAGATAGCGGTGATGTTGGCCATTGTGATGGCAAGCAACCCCAtt                >  4:276201/1‑149 (MQ=255)
                                                                                                                                      agaagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTg               >  4:444865/1‑149 (MQ=255)
                                                                                                                                      agaagCCGCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATAGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTg               <  3:1255964/149‑1 (MQ=255)
                                                                                                                                        aagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCg             >  2:99290/1‑149 (MQ=255)
                                                                                                                                        aagCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCCTTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCg             >  1:957407/1‑149 (MQ=255)
                                                                                                                                          gCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGaa           <  4:817169/149‑1 (MQ=255)
                                                                                                                                          gCCTCCGTAATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGaa           <  2:1211474/149‑1 (MQ=255)
                                                                                                                                           ccTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAc          <  4:1038071/149‑1 (MQ=255)
                                                                                                                                           ccTCCGTTATGGAGGCATCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAc          <  3:444865/149‑1 (MQ=255)
                                                                                                                                            cTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAcc         >  2:1309095/1‑149 (MQ=255)
                                                                                                                                            cTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAcc         >  2:1439325/1‑149 (MQ=255)
                                                                                                                                             tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCa        >  4:190342/1‑149 (MQ=255)
                                                                                                                                             tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCa        >  3:790785/1‑149 (MQ=255)
                                                                                                                                               cGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAAc      >  4:427865/1‑149 (MQ=255)
                                                                                                                                               cGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAAc      <  4:549683/149‑1 (MQ=255)
                                                                                                                                                gTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACa     <  2:997672/149‑1 (MQ=255)
                                                                                                                                                  tATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAgg   <  3:427865/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGC  >  NC_000913/4274596‑4274891

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