Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,907,390 G→C 6.0% R184R (CGC→CGG eno ← enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,907,3900GC 6.0% 233.8 / 6.2 67R184R (CGC→CGGenoenolase
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base C (2/2);  total (33/34)
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.25e-01

GCTTTAACAGCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTC  >  NC_000913/2907247‑2907536
                                                                                                                                               |                                                                                                                                                  
gcTTTAACAGCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGAt                                                                                                                                               <  2:20676/149‑1 (MQ=255)
  tttAACAGCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATgg                                                                                                                                             >  1:114851/1‑149 (MQ=255)
         gCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCttt                                                                                                                                      <  4:179024/149‑1 (MQ=255)
               gCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGt                                                                                                                                <  1:187833/149‑1 (MQ=255)
                cGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGtt                                                                                                                               >  1:410851/1‑149 (MQ=255)
                cGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGtt                                                                                                                               <  3:306515/149‑1 (MQ=255)
                       aGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGcttctt                                                                                                                                       >  4:378182/1‑134 (MQ=255)
                       aGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGcttctt                                                                                                                                       <  3:378182/134‑1 (MQ=255)
                        gCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGcc                                                                                                                       >  3:267500/1‑149 (MQ=255)
                                ttCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCt                                                                                                               >  4:439529/1‑149 (MQ=255)
                                   aGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGACCTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGaa                                                                                                            <  4:109096/149‑1 (MQ=255)
                                            aCCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAAtt                                                                                                   >  3:284635/1‑149 (MQ=255)
                                               cAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCt                                                                                                >  4:282076/1‑149 (MQ=255)
                                                aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTca                                                                                                                                    <  1:3984/112‑1 (MQ=255)
                                                aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTca                                                                                                                                    >  2:3984/1‑112 (MQ=255)
                                                aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTg                                                                                               <  4:245561/149‑1 (MQ=255)
                                                aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTg                                                                                               <  3:277401/149‑1 (MQ=255)
                                                aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTg                                                                                               <  3:15901/149‑1 (MQ=255)
                                                  gTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGa                                                                                             >  4:81401/1‑149 (MQ=255)
                                                       gcgcATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCa                                                                                        >  2:465210/1‑149 (MQ=255)
                                                         gcATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAAc                                                                                      >  4:192538/1‑149 (MQ=255)
                                                          cATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACg                                                                                     <  1:94880/149‑1 (MQ=255)
                                                              gCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACgttgt                                                                                 >  2:58273/1‑149 (MQ=255)
                                                               ccACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACgttgtt                                                                                >  3:168294/1‑149 (MQ=255)
                                                                 aCCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGaa                                                                                                     >  1:61393/1‑126 (MQ=255)
                                                                 aCCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGaa                                                                                                     <  2:61393/126‑1 (MQ=255)
                                                                 aCCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTtgttgt                                                                              >  3:341801/1‑149 (MQ=255)
                                                                 aCCTTCGGCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCGGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACTGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTtgttgt                                                                              >  4:402583/1‑149 (MQ=255)
                                                                              aCAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTcacc                                                                 >  1:76413/1‑149 (MQ=255)
                                                                                   agTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGt                                                            <  3:81401/149‑1 (MQ=255)
                                                                                    gtgtTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGtt                                                           >  1:22795/1‑149 (MQ=255)
                                                                                      gtTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTtga                                                         <  4:326435/149‑1 (MQ=255)
                                                                                         cATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTtgatga                                                      <  4:267500/149‑1 (MQ=255)
                                                                                          aTGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTgatgat                                                     >  2:270523/1‑149 (MQ=255)
                                                                                           tGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGatgatg                                                    <  4:284635/149‑1 (MQ=255)
                                                                                            gCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGata                                                                                           >  1:255376/1‑109 (MQ=255)
                                                                                            gCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGata                                                                                           <  2:255376/109‑1 (MQ=255)
                                                                                                ttCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcat                                               <  2:438341/149‑1 (MQ=255)
                                                                                                ttCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcat                                               <  3:282076/149‑1 (MQ=255)
                                                                                                  cGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcatca                                             >  3:186803/1‑149 (MQ=255)
                                                                                                          gAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAAcc                                     >  3:427152/1‑149 (MQ=255)
                                                                                                             cTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGc                                  >  1:91801/1‑149 (MQ=255)
                                                                                                                ttGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATa                               <  3:439529/149‑1 (MQ=255)
                                                                                                                  gCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                             <  1:58273/149‑1 (MQ=255)
                                                                                                                  gCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                             <  3:250259/149‑1 (MQ=255)
                                                                                                                  gCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                             <  2:76413/149‑1 (MQ=255)
                                                                                                                  gCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                             >  4:439442/1‑149 (MQ=255)
                                                                                                                     aGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACCTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTa                          >  2:191660/1‑149 (MQ=255)
                                                                                                                      ggTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAt                         <  2:114851/149‑1 (MQ=255)
                                                                                                                        tGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAttt                       <  2:22795/149‑1 (MQ=255)
                                                                                                                           tGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGaa                                                                                                            >  2:272483/1‑61 (MQ=255)
                                                                                                                           tGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGaa                                                                                                            <  1:272483/61‑1 (MQ=255)
                                                                                                                             gAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCg                  <  2:410851/149‑1 (MQ=255)
                                                                                                                               aaaCTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGa                >  1:292971/1‑149 (MQ=255)
                                                                                                                               aaaCTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGa                <  1:270523/149‑1 (MQ=255)
                                                                                                                               aaaCTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGa                <  3:402583/149‑1 (MQ=255)
                                                                                                                                aaCTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAg               <  2:292971/149‑1 (MQ=255)
                                                                                                                                  cTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTa             <  4:168294/149‑1 (MQ=255)
                                                                                                                                    tCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTAcc           <  1:465210/149‑1 (MQ=255)
                                                                                                                                    aaaGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACgttgtt                                                                                >  1:273792/3‑80 (MQ=37)
                                                                                                                                    aaaGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACgttgtt                                                                                <  2:273792/78‑1 (MQ=37)
                                                                                                                                      aGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                             <  3:386668/129‑1 (MQ=255)
                                                                                                                                      aGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                             >  4:386668/1‑129 (MQ=255)
                                                                                                                                       gAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgtt        >  4:11296/1‑149 (MQ=255)
                                                                                                                                           ccATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAACTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgttcagt    >  1:225716/1‑149 (MQ=255)
                                                                                                                                             aTGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACGGGCTGAATCTTGAATTCCTGGATATCAACGTTGTTGCCGGCTTGCTCACCACCGTGGATGATGTCCACCATCGGAACCGGCATAGAGTTTTTGCCCGGAGTACCgttcagttc  >  2:405608/1‑149 (MQ=255)
                                                                                                                                             aTGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgttcagttc  >  3:128501/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCTTTAACAGCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTC  >  NC_000913/2907247‑2907536

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