Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,419,721 (G)7→8 intergenic (‑265/+44) ydaG ← / ← racR Rac prophage, uncharacterized protein/Rac prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,419,7141.G98.6% 243.5 / ‑3.6 69intergenic (‑258/+51)ydaG/racRRac prophage, uncharacterized protein/Rac prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (38/30);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGATACTTACCCACTTTCATTTCCATCG  >  NC_000913/1419579‑1419847
                                                                                                                                        |                                                                                                                                     
gAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATggg                                                                                                                                    <  1:816442/139‑1 (MQ=255)
gAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATggg                                                                                                                                    <  2:334803/139‑1 (MQ=255)
     aTTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATgggggggg                                                                                                                               <  2:1011423/139‑2 (MQ=255)
                                        tGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATggg                                                                                                                                    <  1:1606338/99‑1 (MQ=255)
                                        tGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATggg                                                                                                                                    >  2:1606338/1‑99 (MQ=255)
                                           tGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCata                                                                                          >  1:1521275/1‑139 (MQ=255)
                                           tGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCat                                                                                           <  1:394483/138‑1 (MQ=255)
                                           tGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCat                                                                                           >  2:394483/1‑138 (MQ=255)
                                                  ataGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATAtaagttt                                                                                   >  2:1487474/1‑139 (MQ=255)
                                                   taGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   <  1:841072/122‑1 (MQ=255)
                                                   taGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   >  2:841072/1‑122 (MQ=255)
                                                     ggTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATAtaagtttaag                                                                                >  2:1421109/1‑139 (MQ=255)
                                                        aaTGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   >  1:1546646/1‑117 (MQ=255)
                                                        aaTGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   <  2:1546646/117‑1 (MQ=255)
                                                          tGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   >  1:1563742/1‑115 (MQ=255)
                                                          tGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   <  2:1563742/115‑1 (MQ=255)
                                                                attgattgTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATggg                                                                                                                                    <  1:778125/75‑1 (MQ=255)
                                                                attgattgTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATggg                                                                                                                                    >  2:778125/1‑75 (MQ=255)
                                                                      tgTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACg                                                                                                                      >  2:994662/1‑84 (MQ=255)
                                                                      tgTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   >  2:1503540/1‑103 (MQ=255)
                                                                      tgTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   <  1:1503540/103‑1 (MQ=255)
                                                                      tgTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGa                                                                                               >  1:229859/1‑107 (MQ=255)
                                                                      tgTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGa                                                                                               <  2:229859/107‑1 (MQ=255)
                                                                      tgTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATg                                                               >  1:1327724/1‑139 (MQ=255)
                                                                      tgTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAACAGAAAGTGATGGGGGGGGTTAAATAACg                                                                                                                      <  1:994662/84‑1 (MQ=255)
                                                                          aaTAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCatat                                                                                         >  1:583906/1‑109 (MQ=255)
                                                                          aaTAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCatat                                                                                         <  2:583906/109‑1 (MQ=255)
                                                                            tAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTgg                                                         <  2:1521275/139‑1 (MQ=255)
                                                                               gCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCtt                                                      <  1:1487474/139‑1 (MQ=255)
                                                                                  aTGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   >  2:773441/1‑91 (MQ=255)
                                                                                  aTGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   <  1:773441/91‑1 (MQ=255)
                                                                                  aTGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTcc                                                                            >  1:545854/1‑114 (MQ=255)
                                                                                  aTGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTcc                                                                            <  2:545854/114‑1 (MQ=255)
                                                                                   tGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   >  1:1459960/1‑90 (MQ=255)
                                                                                   tGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                                                   <  2:1459960/90‑1 (MQ=255)
                                                                                    gCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCg                                                 >  2:160670/1‑139 (MQ=255)
                                                                                    gCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATgggggggg                                                                                                                               >  1:1241043/1‑59 (MQ=255)
                                                                                    gCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATgggggggg                                                                                                                               <  2:1241043/60‑2 (MQ=255)
                                                                                     ccTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGt                                                <  2:301883/139‑1 (MQ=255)
                                                                                         aTGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGGGTTTTCCGTCAGACCatat                                                                                         >  1:90826/1‑94 (MQ=255)
                                                                                         aTGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCatat                                                                                         >  1:1538006/1‑94 (MQ=255)
                                                                                         aTGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCatat                                                                                         <  2:1538006/94‑1 (MQ=255)
                                                                                         aTGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCatat                                                                                         <  2:90826/94‑1 (MQ=255)
                                                                                              ttGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCa                                       >  1:1165563/1‑139 (MQ=255)
                                                                                                 aGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAaga                                    >  1:368817/1‑139 (MQ=255)
                                                                                                  gCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAagag                                   <  2:1327724/139‑1 (MQ=255)
                                                                                                     tAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGtttt                                                                                                       >  1:368565/1‑68 (MQ=255)
                                                                                                     tAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGtttt                                                                                                       <  2:368565/68‑1 (MQ=255)
                                                                                                     tAACCTAATAGGCGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAAGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGGTTTTTTGATGTTCTTGCTTTTCCGTTCAGATTCAAGAGCttt                               >  1:85680/1‑139 (MQ=255)
                                                                                                                  gATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCa                                                                                                 >  2:426618/1‑61 (MQ=255)
                                                                                                                  gATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCa                                                                                                 <  1:426618/61‑1 (MQ=255)
                                                                                                                  gATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGa                           >  1:473439/1‑131 (MQ=255)
                                                                                                                  gATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGa                           <  2:473439/131‑1 (MQ=255)
                                                                                                                  gATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGATACTTAcc                   >  1:703421/1‑139 (MQ=255)
                                                                                                                  gATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGATACTTAcc                   >  1:1025151/1‑139 (MQ=255)
                                                                                                                  gATGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGATACTTAcc                   >  2:208412/1‑139 (MQ=255)
                                                                                                                   aTGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGa                                                                                               >  1:456451/1‑62 (MQ=255)
                                                                                                                   aTGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGa                                                                                               <  2:456451/62‑1 (MQ=255)
                                                                                                                   aTGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGa                                                                                               >  2:315433/1‑62 (MQ=255)
                                                                                                                   aTGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGa                                                                                               <  2:246105/62‑1 (MQ=255)
                                                                                                                   aTGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGa                                                                                               <  1:315433/62‑1 (MQ=255)
                                                                                                                   aTGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGa                                                                                               >  1:246105/1‑62 (MQ=255)
                                                                                                                   aTGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCtt                                                      <  2:161202/103‑1 (MQ=255)
                                                                                                                   aTGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCtt                                                      >  1:161202/1‑103 (MQ=255)
                                                                                                                   aTGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGATACTTAccc                  >  2:47654/1‑139 (MQ=255)
                                                                                                                   aTGGCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGATACTTAccc                  >  1:228346/1‑139 (MQ=255)
                                                                                                                     ggCGACAGCAGAAAGTGATGGGGGGGGTTAAATaa                                                                                                                        >  2:1271228/1‑35 (MQ=255)
                                                                                                                     ggCGACAGCAGAAAGTGATGGGGGGGGTTAAATaa                                                                                                                        <  1:1271228/35‑1 (MQ=255)
                                                                                                                      gCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTTCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGGTCTGGCTTTTTCGTTCAGATTCAAGAGCTTTCAGATACTTACCCACt               <  2:1165563/139‑1 (MQ=255)
                                                                                                                      gCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTc                                                                             <  2:639219/77‑1 (MQ=255)
                                                                                                                      gCGACAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTc                                                                             >  1:639219/1‑77 (MQ=255)
                                                                                                                         aCAGCAGAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGATACTTACCCACTTTc            >  1:841156/1‑139 (MQ=255)
                                                                                                                          cagcagAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGATACTTACCCACTTTCa           >  1:436607/1‑139 (MQ=255)
                                                                                                                           agcagAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGt                                                                         >  1:1476776/1‑76 (MQ=255)
                                                                                                                           agcagAAAGTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGt                                                                         <  2:1476776/76‑1 (MQ=255)
                                                                                                                                  aGTGATGGGGGGGGTTAAATAACGGAATCCAGGAg                                                                                                           <  2:373494/35‑1 (MQ=255)
                                                                                                                                  aGTGATGGGGGGGGTTAAATAACGGAATCCAGGAg                                                                                                           >  1:373494/1‑35 (MQ=255)
                                                                                                                                   gTGATGGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGATACTTACCCACTTTCATTTCCATCg  <  1:47654/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
GAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAGCTTTCAGATACTTACCCACTTTCATTTCCATCG  >  NC_000913/1419579‑1419847

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