Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,787,434 G→A pseudogene (1038/2253 nt) ygaQ → uncharacterized protein, putative enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,787,4340GA100.0% 237.8 / NA 73pseudogene (1038/2253 nt)ygaQuncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/34);  total (39/34)

ATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGGAA  >  NC_000913/2787304‑2787569
                                                                                                                                   |                                                                                                                                       
atatCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTc                                                                                                                                  >  1:437857/1‑139 (MQ=255)
             ttATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    <  2:630265/139‑1 (MQ=255)
             ttATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    <  1:287416/139‑1 (MQ=255)
                 aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGt                                                                                                                              >  2:96486/1‑125 (MQ=255)
                 aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGt                                                                                                                              <  1:96486/125‑1 (MQ=255)
                 aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattca                                                                                                                >  1:616318/1‑139 (MQ=255)
                 aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattca                                                                                                                >  1:689256/1‑139 (MQ=255)
                      tAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTaaa                                                                                                           >  2:26298/1‑139 (MQ=255)
                       aaCTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                                   <  2:96510/114‑1 (MQ=255)
                       aaCTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                                   >  1:96510/1‑114 (MQ=255)
                       aaCTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAg                                                                                                          >  2:201101/1‑139 (MQ=255)
                       aaCTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAg                                                                                                          >  2:614578/1‑139 (MQ=255)
                        aCTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGa                                                                                                         <  2:370455/139‑1 (MQ=255)
                             aTGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAg                                                                                                    >  2:139636/1‑139 (MQ=255)
                                         ggTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                        <  2:616318/139‑1 (MQ=255)
                                          gTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACg                                                                                       <  1:54921/139‑1 (MQ=255)
                                                   aaaaTGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                              <  1:139636/139‑1 (MQ=255)
                                                        gACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                         >  1:139020/1‑139 (MQ=255)
                                                         aCCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTGAACAAatccatcc                                                                        <  1:748388/139‑1 (MQ=255)
                                                         aCCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  1:607076/139‑1 (MQ=255)
                                                         aCCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  2:215513/139‑1 (MQ=255)
                                                          ccTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                       >  2:622727/1‑139 (MQ=255)
                                                            tGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                     <  2:139020/139‑1 (MQ=255)
                                                            tGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                     <  2:308452/139‑1 (MQ=255)
                                                            tGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAGCCATCCaga                                                                     <  2:575546/139‑1 (MQ=255)
                                                              gATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                                   <  2:108306/139‑1 (MQ=255)
                                                                   tCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                              <  2:465695/123‑1 (MQ=255)
                                                                   tCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                              >  1:465695/1‑123 (MQ=255)
                                                                        ctgct‑gATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                         >  1:406614/1‑139 (MQ=255)
                                                                        ctgct‑gATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                         >  1:793625/1‑139 (MQ=255)
                                                                        ctgct‑gATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                         >  1:720830/1‑139 (MQ=255)
                                                                        ctgct‑gATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                         >  1:686681/1‑139 (MQ=255)
                                                                           ct‑gATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCAc                                                      >  1:172022/1‑139 (MQ=255)
                                                                           ct‑gATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAAAAGATATACATCTCAc                                                      >  2:752497/1‑139 (MQ=255)
                                                                                tataCGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACaaa                                                                               >  2:644879/1‑110 (MQ=255)
                                                                                tataCGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACaaa                                                                               <  1:644879/110‑1 (MQ=255)
                                                                                   aCGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaaa                                               <  1:201101/139‑1 (MQ=255)
                                                                                      aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                                   <  1:606923/52‑1 (MQ=255)
                                                                                      aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                                   >  2:606923/1‑52 (MQ=255)
                                                                                      aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGAc                                            >  2:746476/1‑139 (MQ=255)
                                                                                      aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGAc                                            >  1:757126/1‑139 (MQ=255)
                                                                                      aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATACAGAAGATATACATCTCACTGCAAAAGAc                                            >  1:161272/1‑139 (MQ=255)
                                                                                           aGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGa                                       <  1:753747/139‑1 (MQ=255)
                                                                                           aGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGa                                       <  1:614578/139‑1 (MQ=255)
                                                                                            gATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGaa                                      <  1:218532/139‑1 (MQ=255)
                                                                                             aTCAGTAAATATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAAGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                                     <  2:161272/139‑1 (MQ=255)
                                                                                                    aaGATATAGAAGGTGATATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGCTATCCATCTCACTGCAAACGAATTCGAAGCACGTaa                              >  1:4456/1‑139 (MQ=255)
                                                                                                       atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGa                           >  2:500427/1‑139 (MQ=255)
                                                                                                        tataGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGAt                          <  2:720830/139‑1 (MQ=255)
                                                                                                        tataGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGAt                          <  2:689256/139‑1 (MQ=255)
                                                                                                           aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       <  1:500427/139‑1 (MQ=255)
                                                                                                           aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       <  1:130133/139‑1 (MQ=255)
                                                                                                           aGAAGGTGGTATCCCTAGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       <  2:757126/139‑1 (MQ=255)
                                                                                                            gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACa                      >  1:177612/1‑139 (MQ=255)
                                                                                                                gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       >  1:685418/1‑134 (MQ=255)
                                                                                                                gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       <  2:685418/134‑1 (MQ=255)
                                                                                                                gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACAtctc                  >  1:703064/1‑139 (MQ=255)
                                                                                                                gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACAtctc                  >  1:343498/1‑139 (MQ=255)
                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        >  1:538016/1‑84 (MQ=255)
                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  2:538016/84‑1 (MQ=255)
                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATctct                 >  1:246754/1‑139 (MQ=255)
                                                                                                                  ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTa                >  1:793123/1‑139 (MQ=255)
                                                                                                                  ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTa                >  1:661346/1‑139 (MQ=255)
                                                                                                                   gtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTaa               >  1:595915/1‑139 (MQ=255)
                                                                                                                    tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCa                                                       <  2:393030/98‑1 (MQ=255)
                                                                                                                    tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCa                                                       >  1:393030/1‑98 (MQ=255)
                                                                                                                           ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  2:312712/74‑1 (MQ=255)
                                                                                                                           ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        >  1:312712/1‑74 (MQ=255)
                                                                                                                           ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAt       >  2:177750/1‑139 (MQ=255)
                                                                                                                           ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAt       >  1:95940/1‑139 (MQ=255)
                                                                                                                              tATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                         >  2:476082/1‑86 (MQ=255)
                                                                                                                              tATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                         <  1:476082/86‑1 (MQ=255)
                                                                                                                                ttAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGGaa  <  2:343498/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
ATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCT‑GATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGGAA  >  NC_000913/2787304‑2787569

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