Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 745,001 T→C 5.6% G253G (GGT→GGC ybgK → putative allophanate hydrolase, subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913745,0010TC 5.6% 209.4 / 3.0 73G253G (GGT→GGCybgKputative allophanate hydrolase, subunit 2
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base C (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

CGCATGGGCTATCGCTTACAGGGGCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAAT  >  NC_000913/744858‑745146
                                                                                                                                               |                                                                                                                                                 
cGCATGGGCTATCGCTTACAGGGGCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTa                                                                                                                                              >  5:4017/1‑149 (MQ=255)
      ggCTATCGCTTACAGGGGCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCcgcg                                                                                                                                        >  8:106836/1‑149 (MQ=255)
      ggCTATCGCTTACAGGGGCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCcgc                                                                                                                                         <  6:286867/148‑1 (MQ=255)
       gCTATCGCTTACAGGGGCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGt                                                                                                                                       <  4:326143/149‑1 (MQ=255)
        cTATCGCTTACAGGGGCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTa                                                                                                                                      >  4:245132/1‑149 (MQ=255)
          aTCGCTTACAGCGGCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGAGTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGAAtt                                                                                                                                    >  6:264571/1‑149 (MQ=255)
              cTTACAGGGGCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGGTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCt                                                                                                                                <  7:153476/149‑1 (MQ=255)
                     gggCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGGGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCAt                                                                                                                         <  6:328510/149‑1 (MQ=255)
                       gCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCAtt                                                                                                                        >  2:329980/1‑148 (MQ=255)
                        cacaTTTTAAAACGCACCACCGATCGCGATCTGTTATCTCACGGCTTGTTACCGGGCGTGGTGCAGGTGCCACAAAACGGTCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGTTGACCCGCGTATTGCCTGTATCAttta                                                                                                                      <  8:137256/146‑3 (MQ=255)
                         aaaTTTTAAAACGCACCAACGATCGCGAACTGTGAACTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACAGAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTGACCAGCGGATTGCCTGTATCaattag                                                                                                                     >  2:247289/1‑144 (MQ=255)
                          aaatttAAAACACACAACCGATCGCGAACTGATATCTCACGGTTTGTTACCGGGAGTGGTCCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGGGGTTACCCGCGAATTGCCTGTATCATTGAgg                                                                                                                    >  8:42940/4‑149 (MQ=255)
                          aaTTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAgg                                                                                                                    >  4:307792/1‑149 (MQ=255)
                            ttttAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCt                                                                                                                  <  7:97541/149‑1 (MQ=255)
                              ttAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCt                                                                                                                  >  1:323035/1‑147 (MQ=255)
                              ttAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCt                                                                                                                  <  2:323035/147‑1 (MQ=255)
                               tAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGATGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGa                                                                                                                <  8:19582/148‑1 (MQ=255)
                                 aaaCGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGata                                                                                                              <  1:312250/148‑1 (MQ=255)
                                   aCGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGt                                                                                                           <  1:329980/149‑1 (MQ=255)
                                    cGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTa                                                                                                          <  7:106836/149‑1 (MQ=255)
                                     gCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTAc                                                                                                         <  4:218728/149‑1 (MQ=255)
                                           ccGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTg                                                                                                   <  2:188879/149‑1 (MQ=255)
                                            cGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGa                                                                                                                      >  2:237450/1‑129 (MQ=255)
                                            cGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGa                                                                                                                      <  1:237450/129‑1 (MQ=255)
                                              aTCGCGAACTGTTATCTCAAGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCATCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGgcg                                                                                                >  2:196617/1‑149 (MQ=255)
                                                  cgAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTCACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTa                                                                                                          <  8:337043/135‑1 (MQ=255)
                                                  cgAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTCACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTa                                                                                                          >  7:337043/1‑135 (MQ=255)
                                                    aaCTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGTGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCTGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAAtt                                                                                          <  7:42940/149‑1 (MQ=255)
                                                    aaCTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAAtt                                                                                          <  3:307792/149‑1 (MQ=255)
                                                     aCTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTTTTGATGAACTACACACATACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAAtt                                                                                          <  5:256049/148‑1 (MQ=255)
                                                         ttATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCt                                                                                     >  6:140292/1‑149 (MQ=255)
                                                            tctcACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCg                                                                                   <  3:256076/148‑1 (MQ=255)
                                                            tctcACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACAACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCgg                                                                                  >  3:4816/1‑149 (MQ=255)
                                                              tcACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGt                                                                                 >  3:314237/1‑148 (MQ=255)
                                                                aCGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGACCGACGCACAGACCACAGGTGGTTAGCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAg                                                                              >  4:19222/1‑149 (MQ=255)
                                                                            cGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCAtttt                                                                  <  5:140292/149‑1 (MQ=255)
                                                                               gCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTc                                                               <  3:184120/149‑1 (MQ=255)
                                                                                     tgCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAgtgt                                                         >  8:186706/1‑149 (MQ=255)
                                                                                        aGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCa                                                      <  8:202533/149‑1 (MQ=255)
                                                                                        aGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCa                                                      >  5:228803/1‑149 (MQ=255)
                                                                                          gTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCa                                                      >  1:129858/1‑147 (MQ=255)
                                                                                          gTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCa                                                      <  2:129858/147‑1 (MQ=255)
                                                                                              cacaTAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTgg                                                  >  5:321578/1‑147 (MQ=255)
                                                                                              cacaTAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTgg                                                  >  6:137796/1‑147 (MQ=255)
                                                                                                caTCACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAaga                                              <  2:277666/149‑1 (MQ=255)
                                                                                                 aTAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAagag                                             <  7:212132/149‑1 (MQ=255)
                                                                                                   aaCGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAgg                                            >  6:105709/1‑148 (MQ=255)
                                                                                                   aaCGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACAGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGc                                           <  6:122558/149‑1 (MQ=255)
                                                                                                     cGGGCAGCCGATTGTGGTGATGTACGACGCACAGAAGAACGGGGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCCTCTGGCGCAAATTCCGCTCGGTCAGCCGATTCTTTTTGTCCAGTGTTCACTGGAAGAGGCAc                                         <  4:4816/149‑1 (MQ=255)
                                                                                                      gggCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGCGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTAATTTTGTCCAGTGTTCACTGGAAGAGGCAc                                         >  7:15923/1‑148 (MQ=255)
                                                                                                         cAGCCGATTGTGTTGATGAACGACGCAGAGACCAGCTGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTaaa                                     <  5:264571/149‑1 (MQ=255)
                                                                                                              gATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCAc                                >  8:256677/1‑149 (MQ=255)
                                                                                                              gATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCAc                                <  3:245132/149‑1 (MQ=255)
                                                                                                              gATTGTGTTGATGAACGACGAACATACCACCGGTGGTTACCCGCTTATTGCCTGTATCATTTAGGCTGATATGTAACATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCAc                                <  8:15923/149‑1 (MQ=255)
                                                                                                                    gtTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGc                             >  1:299224/1‑146 (MQ=255)
                                                                                                                       gatgaACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGAt                        >  2:340600/1‑148 (MQ=255)
                                                                                                                       gatgaACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCAGCTGGCGCAAATTACGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGATGCACTAAAAGCACGGCAAgct                        >  6:325350/1‑146 (MQ=255)
                                                                                                                       gatgaACGACGCACAGACCACCGGCGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCAtt                                                                    >  4:130665/1‑104 (MQ=255)
                                                                                                                       gatgaACGACGCACAGACCACCGGCGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCAtt                                                                    >  4:125889/1‑104 (MQ=255)
                                                                                                                       gatgaACGACGCACAGACCACCGGCGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCAtt                                                                    <  3:130665/104‑1 (MQ=255)
                                                                                                                       gatgaACGACGCACAGACCACCGGCGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCAtt                                                                    <  3:125889/104‑1 (MQ=255)
                                                                                                                         tgaACGACGCACAGAACACCGTTGGTTACACGAGTAATGCCTCTAGCATTGAGGCGGATATGTACCATCTGGCGCATATTACGCTCGGTCAGCCGATTCATTTTTTACAGTGTTCACTGGAAGAGGCACTAAAAGCACCGCAAGATCAg                     <  1:247289/149‑1 (MQ=255)
                                                                                                                          gaACGACGCACAGACCACCGGTGGTTACACGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAg                     >  5:102096/1‑148 (MQ=255)
                                                                                                                               acgCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGt               >  7:162347/1‑149 (MQ=255)
                                                                                                                                  cacaGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAt            >  8:302775/1‑149 (MQ=255)
                                                                                                                                   acaGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGTCTGATATGTACCATCTGGCGCAAATTCCGCTCGCTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAAccaaatt           >  8:291337/1‑143 (MQ=255)
                                                                                                                                      gACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtttt         >  5:135773/1‑148 (MQ=255)
                                                                                                                                       accaccGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAgg                                            >  8:110716/1‑112 (MQ=255)
                                                                                                                                       accaccGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAgg                                            <  7:110716/112‑1 (MQ=255)
                                                                                                                                        ccaccGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGaa      <  2:299224/149‑1 (MQ=255)
                                                                                                                                         caccGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGCTGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGaa      >  5:272117/1‑148 (MQ=255)
                                                                                                                                         caccGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAt            >  5:30359/1‑142 (MQ=255)
                                                                                                                                         caccGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAt            <  6:30359/142‑1 (MQ=255)
                                                                                                                                            ctgtggtTACCCGCGTATTGCCTGTATCATTGATGCTGATATTTACCATCTGGCGCAATTTCCGCTCGGTCAGACGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAAt  <  1:196617/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CGCATGGGCTATCGCTTACAGGGGCAAATTTTAAAACGCACCACCGATCGCGAACTGTTATCTCACGGTTTGTTACCGGGCGTGGTGCAGGTGCCACATAACGGGCAGCCGATTGTGTTGATGAACGACGCACAGACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAAT  >  NC_000913/744858‑745146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: