Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,984,411 C→T *394* (TAG→TAA yqeF ← short chain acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4110CT100.0% 243.8 / NA 76*394* (TAG→TAAyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (39/37);  total (39/37)

GTGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTA  >  NC_000913/2984276‑2984544
                                                                                                                                       |                                                                                                                                       
gTGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTAt                                                                                                                                      <  1:790302/139‑1 (MQ=255)
                gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATgg                                                                                                                      >  2:476654/1‑139 (MQ=255)
                     tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAAt                                                                                                                 >  2:1338232/1‑139 (MQ=255)
                          tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCcac                                                                                                            >  2:945636/1‑139 (MQ=255)
                          tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCcac                                                                                                            >  1:1449552/1‑139 (MQ=255)
                              aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑ccc                                                                                                      >  2:58673/1‑139 (MQ=255)
                                tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTg                                                                                                    >  2:1022216/1‑139 (MQ=255)
                                tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATACCACA‑‑CCCTg                                                                                                    >  1:480121/1‑139 (MQ=255)
                                   aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTCTCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGcc                                                                                                 >  1:1043516/1‑139 (MQ=255)
                                   aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGcc                                                                                                 <  1:1338232/139‑1 (MQ=255)
                                       aTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATACCACA‑‑CCCTGGCCCccg                                                                                             <  2:480121/139‑1 (MQ=255)
                                        tGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCCCTGGccccc                                                                                              >  1:1153650/1‑139 (MQ=255)
                                            ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                                     >  1:455308/1‑112 (MQ=255)
                                            ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                                     <  2:455308/112‑1 (MQ=255)
                                                 ccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACaaa                                                                                   <  1:945636/139‑1 (MQ=255)
                                                   aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGc                                                                                 <  2:1449552/139‑1 (MQ=255)
                                                               tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                              <  1:1385417/100‑1 (MQ=255)
                                                               tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                              >  2:1385417/1‑100 (MQ=255)
                                                               tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGCAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                     >  1:135665/1‑139 (MQ=255)
                                                                 gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                   >  1:1140213/1‑139 (MQ=255)
                                                                   ctggccagccaTAATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGTTCAATGCCAGA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                 <  2:1043516/136‑1 (MQ=255)
                                                                   cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                 <  1:58673/139‑1 (MQ=255)
                                                                    ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGcccc                                                                                               >  2:23826/1‑108 (MQ=255)
                                                                    ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGcccc                                                                                               <  1:23826/108‑1 (MQ=255)
                                                                    ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                <  1:1667960/139‑1 (MQ=255)
                                                                     tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCca                                                                                                             <  1:1102710/95‑1 (MQ=255)
                                                                     tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCca                                                                                                             >  2:1102710/1‑95 (MQ=255)
                                                                                 ttattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                   <  1:476654/139‑1 (MQ=255)
                                                                                         cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                           >  1:470272/1‑139 (MQ=255)
                                                                                            aTTAATGACGGCAGTTATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGCTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGc                                                                                 <  2:437524/98‑1 (MQ=255)
                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                              >  2:614625/1‑71 (MQ=255)
                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                              <  1:614625/71‑1 (MQ=255)
                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTg                                                                             <  1:173256/102‑1 (MQ=255)
                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTg                                                                             >  2:173256/1‑102 (MQ=255)
                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                        >  2:760467/1‑139 (MQ=255)
                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                        >  2:187080/1‑139 (MQ=255)
                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                        >  2:458909/1‑139 (MQ=255)
                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGCTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGc                                                                                 >  1:437524/1‑98 (MQ=255)
                                                                                             ttAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                       <  1:760467/139‑1 (MQ=255)
                                                                                              tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCa                                                                          <  1:1284911/103‑1 (MQ=255)
                                                                                              tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCa                                                                          >  2:1284911/1‑103 (MQ=255)
                                                                                              tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                  <  1:1130235/111‑1 (MQ=255)
                                                                                              tAATGACGGCAGATAGTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                  >  2:1130235/1‑111 (MQ=255)
                                                                                               aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                     >  1:114201/1‑139 (MQ=255)
                                                                                                   aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                 <  2:470272/139‑1 (MQ=255)
                                                                                                    cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                <  2:1140213/139‑1 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCaa                                                                                                                  >  1:1586369/1‑57 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCaa                                                                                                                  <  2:1586369/57‑1 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                              >  2:1658577/1‑61 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                              <  1:1658577/61‑1 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCTCA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                        <  2:1600216/129‑1 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGc                                                                                 >  1:1567475/1‑88 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGc                                                                                 <  2:1567475/88‑1 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                        >  1:1600216/1‑129 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                               <  1:938852/138‑1 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                               >  2:938852/1‑138 (MQ=255)
                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                              >  1:9377/1‑139 (MQ=255)
                                                                                                       cAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa                             >  1:1206767/1‑139 (MQ=255)
                                                                                                        aGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                            <  2:9377/139‑1 (MQ=255)
                                                                                                        aGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                            <  2:114201/139‑1 (MQ=255)
                                                                                                          atatTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                            <  2:1153650/139‑1 (MQ=255)
                                                                                                                 tCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                                   >  2:1348098/1‑45 (MQ=255)
                                                                                                                 tCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                                   <  1:1348098/45‑1 (MQ=255)
                                                                                                                 tCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                   <  1:1022216/139‑1 (MQ=255)
                                                                                                                     tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCca                                                                                                             >  2:1608175/1‑47 (MQ=255)
                                                                                                                     tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCca                                                                                                             <  1:1608175/47‑1 (MQ=255)
                                                                                                                      aaGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc              >  2:527644/1‑139 (MQ=255)
                                                                                                                          aggaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCTCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa          >  1:1610021/1‑139 (MQ=255)
                                                                                                                          aggaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa          <  1:187080/139‑1 (MQ=255)
                                                                                                                            gaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg        <  2:135665/139‑1 (MQ=255)
                                                                                                                            gaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg        <  2:1206767/139‑1 (MQ=255)
                                                                                                                               aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc     >  1:808179/1‑139 (MQ=255)
                                                                                                                               aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc     >  2:1560572/1‑139 (MQ=255)
                                                                                                                                  cAGAGTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTa  <  2:808179/139‑1 (MQ=255)
                                                                                                                                    gagTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                        >  2:1238210/1‑99 (MQ=255)
                                                                                                                                    gagTTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                        <  1:1238210/99‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
GTGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACA‑‑CCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTA  >  NC_000913/2984276‑2984544

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