Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,274,740 G→A 6.2% intergenic (+79/+20) 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,7400GA 6.2% 613.7 / 19.6 187intergenic (+79/+20)ssb/yjcBsingle‑stranded DNA‑binding protein/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (91/83);  new base A (7/5);  total (99/88)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 189 total aligned reads displayed.

GGTGGCGCTCCGGCAGGTGGCAATATCGGTGGTGGTCAGCCGCAGGGCGGTTGGGGTCAGCCTCAGCAGCCGCAGGGTGGCAATCAGTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACG  >  NC_000913/4274479‑4274988
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
ggtggCGCTCCGGCAGGTGGCAATATCGGTGGTGGTCAGCCGCAGGGCGGTTGGGGTCAGCCTCAGCAGCCGCAGGGTGGCAATCAGTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTACAACAATAGGTTATCTTGTTTTAAGGGGGATGATTAAGCCAATCACCAGCCATAAAAAGACAGCCTCCATTCAGGAGGCTTc                                                                                                                                                                                                                                             >  4:253821/1‑275 (MQ=255)
                                                       gTCAGCCTCAGCAGCCGCAGGGTGGCAATCAGTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGGCTGATTCATCtcttct                                                                                                                                                                                      >  2:341643/1‑275 (MQ=255)
                                                               gagcagcCGAAGGGGGGAAGTCAGTTCGGCGGCGGCGGGCAGTCTGGCCGGCAGCAGTCCGATCCGGCATCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                                                                                                                              <  4:97034/273‑1 (MQ=255)
                                                                                    cggTTCGGCGGCGGCGCGCAGGCTCGCCGTCAGCAGTCCGCTCCGGCGGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGa                                                                                                                                                          <  4:86654/272‑1 (MQ=255)
                                                                                      gTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAAAAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCAt                                                                                                                                                                                                                        >  1:52583/1‑210 (MQ=255)
                                                                                                                  cagcagTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGt                                                                                                                                                                                                                                                         <  6:898746/149‑1 (MQ=255)
                                                                                                                     cagTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAt                                                                                                                                                                                                                                                      <  5:81487/149‑1 (MQ=255)
                                                                                                                             tCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCtt                                                                                                                                                                                                                                              <  8:500569/149‑1 (MQ=255)
                                                                                                                                ggCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCTTCtt                                                                                                                                                                                                                                           <  5:658922/149‑2 (MQ=255)
                                                                                                                                 gCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCTTCttt                                                                                                                                                                                                                                          <  6:1364897/149‑3 (MQ=255)
                                                                                                                                   aGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGt                                                                                                                                                                                                                                        <  8:616972/149‑1 (MQ=255)
                                                                                                                                          tCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGt                                                                                                                                                                                                                                 >  7:727859/1‑149 (MQ=255)
                                                                                                                                               cGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCaaaa                                                                                                                                                                                                                            <  5:296115/149‑1 (MQ=255)
                                                                                                                                                    cgccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGCATCAGGTCAAAATCAtt                                                                                                                                                                                                                       <  5:476877/149‑1 (MQ=255)
                                                                                                                                                      ccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTgg                                                                                                                                                                                                                     >  6:40166/1‑149 (MQ=255)
                                                                                                                                                         aTGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCa                                                                                                                                                                                                                  >  6:957936/1‑149 (MQ=255)
                                                                                                                                                               tttGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTgg                                                                                                                                                                                                            >  7:685343/1‑149 (MQ=255)
                                                                                                                                                                 tgatgatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGgtgt                                                                                                                                                                                                      >  4:184388/1‑153 (MQ=255)
                                                                                                                                                                     gatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                                                                                                                                  <  2:184388/153‑1 (MQ=255)
                                                                                                                                                                       tgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCt                                                                                                                                                                                                    <  8:139224/149‑1 (MQ=255)
                                                                                                                                                                        gaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCAt                                                                                                                                                                                                                        <  5:1104101/128‑1 (MQ=255)
                                                                                                                                                                        gaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCAt                                                                                                                                                                                                                        >  6:1104101/1‑128 (MQ=255)
                                                                                                                                                                        gaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTg                                                                                                                                                                                                   >  5:209299/1‑149 (MQ=255)
                                                                                                                                                                          cATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAt                                                                                                                                                                                                 <  8:629534/149‑1 (MQ=255)
                                                                                                                                                                           aTTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAtt                                                                                                                                                                                                <  5:778643/149‑1 (MQ=255)
                                                                                                                                                                            ttCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTc                                                                                                                                                                                               >  8:421542/1‑149 (MQ=255)
                                                                                                                                                                             tCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCa                                                                                                                                                                                              >  5:1305215/1‑149 (MQ=255)
                                                                                                                                                                                   cTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtctt                                                                                                                                                                                        <  5:1358815/149‑1 (MQ=255)
                                                                                                                                                                                   cTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGCATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtctt                                                                                                                                                                                        <  8:1100279/149‑1 (MQ=255)
                                                                                                                                                                                    tGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttc                                                                                                                                                                                       <  6:1157321/149‑1 (MQ=255)
                                                                                                                                                                                       tttGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTat                                                                                                                                                                                    >  8:460460/1‑149 (MQ=255)
                                                                                                                                                                                           tCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTatcatc                                                                                                                                                                                >  6:66767/1‑149 (MQ=255)
                                                                                                                                                                                           tCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTatcatc                                                                                                                                                                                <  5:1143366/149‑1 (MQ=255)
                                                                                                                                                                                            cATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTAtcatca                                                                                                                                                                               <  6:797192/149‑1 (MQ=255)
                                                                                                                                                                                             aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                                                                                                                              >  6:369600/1‑149 (MQ=255)
                                                                                                                                                                                             aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                                                                                                                              >  7:921671/1‑149 (MQ=255)
                                                                                                                                                                                             aTTAAAACAATAGGTTATATTGTTTTAAGGAGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTAATGGGGCGTTCTTGGTATCGGGTCAAAATCTTTGGCCTTTGTGGCGGTTTTTGATTCATTTCTTCTATCATCAc                                                                                                                                                                              >  7:1113360/1‑149 (MQ=255)
                                                                                                                                                                                              ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtctc                                                                                                                                                                                          >  5:612799/1‑136 (MQ=255)
                                                                                                                                                                                              ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtctc                                                                                                                                                                                          <  6:612799/136‑1 (MQ=255)
                                                                                                                                                                                               tAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCg                                                                                                                                                                            >  8:389787/1‑149 (MQ=255)
                                                                                                                                                                                                         ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgtt                                                                                                                                                                  <  7:968570/149‑1 (MQ=255)
                                                                                                                                                                                                         ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgtt                                                                                                                                                                  >  7:1292028/1‑149 (MQ=255)
                                                                                                                                                                                                          gTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAtt                                                                                                                                                                                                <  5:316597/118‑1 (MQ=255)
                                                                                                                                                                                                           ttATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGTGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttgt                                                                                                                                                                >  8:840456/1‑149 (MQ=255)
                                                                                                                                                                                                           ttATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttgt                                                                                                                                                                <  8:685343/149‑1 (MQ=255)
                                                                                                                                                                                                                ttGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                                                                                                                                  <  7:547890/110‑1 (MQ=255)
                                                                                                                                                                                                                ttGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAg                                                                                                                                                           >  7:407037/1‑149 (MQ=255)
                                                                                                                                                                                                                   ttttAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATa                                                                                                                                                        >  7:315384/1‑149 (MQ=255)
                                                                                                                                                                                                                       aaGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCgg                                                                                                                                                    <  5:270523/149‑1 (MQ=255)
                                                                                                                                                                                                                        aGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGt                                                                                                                                                   <  7:1211426/149‑1 (MQ=255)
                                                                                                                                                                                                                          gTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGa                                                                                                                                                 <  8:1292028/149‑1 (MQ=255)
                                                                                                                                                                                                                           tGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGAt                                                                                                                                                >  5:929609/1‑149 (MQ=255)
                                                                                                                                                                                                                             gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAt                                                                                                                                                                                                 >  4:12766/1‑98 (MQ=255)
                                                                                                                                                                                                                             gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGt                                                                                                                                              >  8:533630/1‑149 (MQ=255)
                                                                                                                                                                                                                             gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGt                                                                                                                                              >  7:817169/1‑149 (MQ=255)
                                                                                                                                                                                                                              atgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTc                                                                                                                                             >  8:1383272/1‑149 (MQ=255)
                                                                                                                                                                                                                                gatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCg                                                                                                                                           >  6:1365658/1‑149 (MQ=255)
                                                                                                                                                                                                                                 atTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAt                                                                                                                                                                                             <  2:12766/98‑1 (MQ=255)
                                                                                                                                                                                                                                 atTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGc                                                                                                                                          >  7:581766/1‑149 (MQ=255)
                                                                                                                                                                                                                                  ttAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGcc                                                                                                                                         >  5:997672/1‑149 (MQ=255)
                                                                                                                                                                                                                                   tAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGc                                                                                                                                          <  8:1220060/147‑1 (MQ=255)
                                                                                                                                                                                                                                   tAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCa                                                                                                                                        <  7:52567/149‑1 (MQ=255)
                                                                                                                                                                                                                                    aaaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTCTCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCAt                                                                                                                                       >  6:286508/1‑149 (MQ=255)
                                                                                                                                                                                                                                    aaaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGACGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGAt                                                                                                                                                <  5:3305/140‑1 (MQ=255)
                                                                                                                                                                                                                                       gCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATtgt                                                                                                                                    <  7:389787/149‑1 (MQ=255)
                                                                                                                                                                                                                                       gCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATtgt                                                                                                                                    <  7:620942/149‑1 (MQ=255)
                                                                                                                                                                                                                                       gCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATtgt                                                                                                                                    <  6:581251/149‑1 (MQ=255)
                                                                                                                                                                                                                                         aTCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGa                                                                                                                                  >  6:1298631/1‑149 (MQ=255)
                                                                                                                                                                                                                                          tCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATg                                                                                                                                               >  7:1037968/1‑135 (MQ=255)
                                                                                                                                                                                                                                           cTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGtgtg                       <  4:100480/254‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTATGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGtgtg                       >  2:100480/1‑254 (MQ=255)
                                                                                                                                                                                                                                            tGCCCGCCATAAAAAAGAAGCCTCCGTTTAGGGGGCTTCTACGGTTCAGGTCAAAATCATTGGTCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCTGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACg  <  4:341643/274‑1 (MQ=255)
                                                                                                                                                                                                                                             gccagccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGcc                                                                                                                              <  5:369600/149‑1 (MQ=255)
                                                                                                                                                                                                                                              ccagccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAt                                                                                                                                                                                             <  7:442332/85‑1 (MQ=255)
                                                                                                                                                                                                                                              ccagccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCAATTGTGGCGGTGTCTGATTCAt                                                                                                                                                                                             >  8:442332/1‑85 (MQ=255)
                                                                                                                                                                                                                                               cagccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCaa                                                                                                                            <  8:921671/149‑1 (MQ=255)
                                                                                                                                                                                                                                               cagccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCaa                                                                                                                            <  5:1215275/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   caTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCaa                                                                                                                        <  7:285921/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   caTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCaa                                                                                                                        <  6:1449111/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    aTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGc                                                                                                                                          <  6:318019/130‑1 (MQ=255)
                                                                                                                                                                                                                                                    aTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAc                                                                                                                       <  7:262762/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     tAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAcc                                                                                                                      >  7:1396882/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        aaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCa                                                                                                                   >  5:459309/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         aaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCAt                                                                                                                  <  8:727859/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          aagaagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCAtt                                                                                                                 >  8:1173614/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          aagaagCCTACTTTATGGAGACTTCTACATATCAAGTCAAAATCATTGGCCATTGTGCCGGTGTCTGATTCGTCTCTTCTATCATCACCGGCTCGATGTTGTGCCAGATAGCGGTGATGTTGGCCATTGTGATGGCAAGCAACCCCAtt                                                                                                                 >  8:276201/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           agaagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTg                                                                                                                >  8:444865/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           agaagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTg                                                                                                                <  5:66767/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           agaagCCGCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATAGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTg                                                                                                                <  7:1255964/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            gaagCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGc                                                                                                               >  5:522139/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               gCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGaa                                                                                                            <  8:817169/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                ccTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAc                                                                                                           <  8:1038071/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                ccTCCGTTATGGAGGCATCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAc                                                                                                           <  7:444865/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAcc                                                                                                          >  6:1439325/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 cTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAcc                                                                                                          <  7:1383272/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCa                                                                                                         >  7:1246886/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCa                                                                                                         >  8:190342/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCa                                                                                                         >  7:790785/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    cGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAAc                                                                                                       <  8:549683/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACa                                                                                                      <  6:997672/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
GGTGGCGCTCCGGCAGGTGGCAATATCGGTGGTGGTCAGCCGCAGGGCGGTTGGGGTCAGCCTCAGCAGCCGCAGGGTGGCAATCAGTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACG  >  NC_000913/4274479‑4274988

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