Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,274,740 G→A 7.3% 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 7.3% 740.8 / 31.8 232intergenic (+79/+20)ssb/yjcBsingle‑stranded DNA‑binding protein/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (114/100);  new base A (9/8);  total (124/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

Only 100 of 239 total aligned reads displayed.

GCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACA  >  NC_000913/4274592‑4274888
                                                                                                                                                    |                                                                                                                                                     
gcagcaGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCg                                                                                                                                                       <  2:1230990/149‑1 (MQ=255)
       tCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGa                                                                                                                                               >  2:1000877/1‑149 (MQ=255)
       tCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGa                                                                                                                                               <  1:222327/149‑1 (MQ=255)
          gCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGc                                                                                                                                            <  3:620018/149‑1 (MQ=255)
          gCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAAC‑GGAGGc                                                                                                                                            <  2:516627/149‑1 (MQ=255)
            gccGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCtt                                                                                                                                          >  2:1490697/2‑149 (MQ=255)
                     gcCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATc                                                                                                                                 <  1:1723086/149‑1 (MQ=255)
                      ccGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAAC‑GGAGGCTTc                                                                                                                                         >  1:1098369/1‑140 (MQ=255)
                            aaCGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCaa                                                                                                                          >  1:474927/1‑149 (MQ=255)
                             aCGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCaaa                                                                                                                         >  3:1381487/1‑149 (MQ=255)
                                 gccgccGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCa                                                                                                                     <  2:840383/149‑1 (MQ=255)
                                   cgccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCAtt                                                                                                                   <  2:1070508/149‑1 (MQ=255)
                                   cgccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAAC‑GGAGGCTTCtttttt                                                                                                                                   <  1:1595476/133‑6 (MQ=255)
                                    gccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTACCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTg                                                                                                                  >  2:1289893/1‑149 (MQ=255)
                                     ccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTgg                                                                                                                 <  4:450441/149‑1 (MQ=255)
                                              tttGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAGGGTGGATGATTAAACAACCTGCCAGCCATAAAAAAAAAG‑CTCCGGTATGGGGGGCTTCTAGGATTAAGGCCAAAA‑CATTGGGCAtttttgg                                                                                                       >  1:961494/1‑144 (MQ=255)
                                              tttGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTgg                                                                                                        >  1:1337126/1‑149 (MQ=255)
                                                tgatgatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCg                                                                                                      <  1:1372227/149‑1 (MQ=255)
                                                    gatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCa                                                                                                                     >  4:324713/1‑130 (MQ=255)
                                                    gatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCa                                                                                                                     <  3:324713/130‑1 (MQ=255)
                                                      tgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCt                                                                                                >  2:1089658/1‑149 (MQ=255)
                                                       gaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTg                                                                                               >  3:1090531/1‑149 (MQ=255)
                                                        ataTTCCCTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCTTAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAAACATTGGCCATTGTGGCGGTGTCTGa                                                                                              <  1:1702391/147‑1 (MQ=255)
                                                        aCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                              >  3:1506122/1‑149 (MQ=255)
                                                          aTTCCGTTCTGATTTGTCATTAAAACCATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTc                                                                                                 >  3:776096/1‑144 (MQ=255)
                                                          aTTCCGTTCTGATTTGTCATTAAAACCATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTc                                                                                                 <  4:776096/144‑1 (MQ=255)
                                                                ttCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtctc                                                                                      >  4:905155/1‑149 (MQ=255)
                                                                 tCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCAtt                                                                                                                   <  2:1748848/119‑1 (MQ=255)
                                                                 tCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCAtt                                                                                                                   >  1:1748848/1‑119 (MQ=255)
                                                                  cTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtctt                                                                                    <  1:1463046/149‑1 (MQ=255)
                                                                   tGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCa                                                                                                                     <  3:467031/115‑1 (MQ=255)
                                                                   tGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtc                                                                                        >  4:725152/1‑144 (MQ=255)
                                                                    gATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                              <  1:1279269/137‑1 (MQ=255)
                                                                      tttGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTat                                                                                >  3:626432/1‑149 (MQ=255)
                                                                          tCATTAAAATAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATATCTTCTatcatc                                                                            >  2:340332/1‑149 (MQ=255)
                                                                          tCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTatcatc                                                                            <  2:1266471/149‑1 (MQ=255)
                                                                            aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                          <  1:1000877/149‑1 (MQ=255)
                                                                            aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                          >  3:191561/1‑149 (MQ=255)
                                                                            aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                          >  2:1623525/1‑149 (MQ=255)
                                                                            aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                          >  3:819586/1‑149 (MQ=255)
                                                                            aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                          >  2:45057/1‑149 (MQ=255)
                                                                             ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAcc                                                                         >  4:14684/1‑149 (MQ=255)
                                                                             ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAcc                                                                         <  4:819586/149‑1 (MQ=255)
                                                                               aaaaCAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGa                                                                       >  3:782646/1‑149 (MQ=255)
                                                                                 aaCAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACt                                                                     >  3:1450478/1‑149 (MQ=255)
                                                                                   cAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTatc                                                                               >  2:1138277/1‑137 (MQ=255)
                                                                                   cAATAGGTTATATTGTTTAAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTatc                                                                               <  1:1138277/137‑1 (MQ=255)
                                                                                    aaTAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                              >  3:957932/1‑121 (MQ=255)
                                                                                    aaTAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                              <  4:957932/121‑1 (MQ=255)
                                                                                        ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgtt                                                              <  4:108960/149‑1 (MQ=255)
                                                                                         gTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAg                                                                                                                               <  2:638585/83‑1 (MQ=255)
                                                                                         gTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttg                                                             >  2:1171114/1‑149 (MQ=255)
                                                                                              atTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCa                                                        >  2:927182/1‑149 (MQ=255)
                                                                                               ttGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGGTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTAACATCACCGACTAGATGTTGTGCCAg                                                       >  2:1570287/1‑149 (MQ=255)
                                                                                                  ttttAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATa                                                    >  4:1718479/1‑149 (MQ=255)
                                                                                                     tAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCg                                                 <  3:420126/149‑1 (MQ=255)
                                                                                                      aaGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCgg                                                <  4:442910/149‑1 (MQ=255)
                                                                                                      aaGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCgg                                                <  1:150594/149‑1 (MQ=255)
                                                                                                       aGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGt                                               <  2:1154086/149‑1 (MQ=255)
                                                                                                          tGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGAt                                            >  1:53460/1‑149 (MQ=255)
                                                                                                          tGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGAt                                            >  3:1607302/1‑149 (MQ=255)
                                                                                                            gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGt                                          >  4:1348244/1‑149 (MQ=255)
                                                                                                            gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGt                                          <  4:543384/149‑1 (MQ=255)
                                                                                                                  tAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTTTATCATCATCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGGCa                                    >  4:1396137/1‑149 (MQ=255)
                                                                                                                    aaGCATCTGCCAGCCATAAAAGAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttgt                                                            >  1:1220486/1‑123 (MQ=255)
                                                                                                                    aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAg                                                                                                                                              >  3:1181049/1‑41 (MQ=255)
                                                                                                                    aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAg                                                                                                                                              <  4:1181049/41‑1 (MQ=255)
                                                                                                                    aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCAtt                                  <  1:883217/149‑1 (MQ=255)
                                                                                                                    aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCAtt                                  >  1:592468/1‑149 (MQ=255)
                                                                                                                    aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCAtt                                  <  2:1723447/149‑1 (MQ=255)
                                                                                                                    aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCAtt                                  >  3:1298876/1‑149 (MQ=255)
                                                                                                                      gCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATtgt                                >  2:541363/1‑149 (MQ=255)
                                                                                                                      gCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATtgt                                <  4:1557354/149‑1 (MQ=255)
                                                                                                                        aTCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGa                              <  2:961494/149‑1 (MQ=255)
                                                                                                                         tCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                             <  2:474927/149‑1 (MQ=255)
                                                                                                                         tCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                             >  2:794377/1‑149 (MQ=255)
                                                                                                                         tCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                             >  2:1461672/1‑149 (MQ=255)
                                                                                                                         tCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                             <  1:460350/149‑1 (MQ=255)
                                                                                                                              cagccaTAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCaa                        <  2:691820/149‑1 (MQ=255)
                                                                                                                              cagccaTAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCaa                        <  4:626432/149‑1 (MQ=255)
                                                                                                                                gccaTAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGa                                                                  >  4:344613/1‑105 (MQ=255)
                                                                                                                                  caTAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAGCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCaa                    <  3:1494594/149‑1 (MQ=255)
                                                                                                                                  caTAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCaa                    >  4:516632/1‑149 (MQ=255)
                                                                                                                                    tAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGa                                             >  1:1391541/1‑122 (MQ=255)
                                                                                                                                    tAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGa                                             <  2:1391541/122‑1 (MQ=255)
                                                                                                                                        gtagaagCCTCCATAAC‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCAt              >  4:1610092/3‑149 (MQ=255)
                                                                                                                                        aaaGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCAt              <  1:927182/149‑1 (MQ=255)
                                                                                                                                         tagaagCCTCCATAAC‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgtt                                                              <  1:1265956/99‑1 (MQ=39)
                                                                                                                                         aagaagCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCAtt             <  4:1298876/149‑1 (MQ=255)
                                                                                                                                         aagaagCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCAt              <  1:263485/148‑1 (MQ=255)
                                                                                                                                         aagaagCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCAGTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCAtt             <  1:756614/149‑1 (MQ=255)
                                                                                                                                               ccTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAc       <  3:1530670/149‑1 (MQ=255)
                                                                                                                                                cTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTg                                              <  1:1167287/109‑1 (MQ=255)
                                                                                                                                                cTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAcc      <  4:1506786/149‑1 (MQ=255)
                                                                                                                                                cTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAcc      >  2:397343/1‑149 (MQ=255)
                                                                                                                                                cTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAcc      >  4:1097459/1‑149 (MQ=255)
                                                                                                                                                cTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAcc      >  4:1018082/1‑149 (MQ=255)
                                                                                                                                                cTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAcc      >  1:845743/1‑149 (MQ=255)
                                                                                                                                                 tCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCa     >  1:1496685/1‑149 (MQ=255)
                                                                                                                                                    gTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACa  <  3:1018082/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAT‑GGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACA  >  NC_000913/4274592‑4274888

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