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% 324.5 / NA 98pseudogene (1038/2253 nt)ygaQuncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (50/48);  total (50/48)

ATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGGAA  >  NC_000913/2787304‑2787569
                                                                                                                                  |                                                                                                                                       
atatCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCa                                                                                                                                 <  2:505410/139‑1 (MQ=255)
     aTTCCTGATCATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                            <  1:951663/139‑1 (MQ=255)
       tCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTg                                                                                                                          >  2:645517/1‑139 (MQ=255)
       tCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTg                                                                                                                          >  2:1542349/1‑139 (MQ=255)
             ttATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    <  2:1182191/139‑1 (MQ=255)
                taCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattc                                                                                                                 >  1:392203/1‑139 (MQ=255)
                 aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattca                                                                                                                >  1:1313227/1‑139 (MQ=255)
                      tAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTa                                                                                                             >  1:465879/1‑137 (MQ=255)
                      tAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTa                                                                                                             <  2:465879/137‑1 (MQ=255)
                             aTGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                            >  2:929502/1‑115 (MQ=255)
                             aTGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                            <  1:929502/115‑1 (MQ=255)
                             aTGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                      >  2:1315358/1‑137 (MQ=255)
                             aTGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                      <  1:1315358/137‑1 (MQ=255)
                             aTGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAg                                                                                                    >  1:1225779/1‑139 (MQ=255)
                                    ttGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                             >  2:1136695/1‑139 (MQ=255)
                                      gggggTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAg                                                                                           >  2:33604/1‑139 (MQ=255)
                                         ggTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                        <  2:1313227/139‑1 (MQ=255)
                                                   aaaaTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTaaa                                                                                   <  1:149161/134‑1 (MQ=255)
                                                   aaaaTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTaaa                                                                                   >  2:149161/1‑134 (MQ=255)
                                                   aaaaTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                              <  1:1542349/139‑1 (MQ=255)
                                                        gACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    >  2:674689/1‑96 (MQ=255)
                                                        gACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    <  1:674689/96‑1 (MQ=255)
                                                        gACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                      >  1:743792/1‑110 (MQ=255)
                                                        gACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                      <  2:743792/110‑1 (MQ=255)
                                                         aCCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  2:392203/139‑1 (MQ=255)
                                                          ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                       >  1:1346183/1‑139 (MQ=255)
                                                           cttgATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                                                                      <  1:1136695/136‑1 (MQ=255)
                                                           cTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACATATCCATCCag                                                                      <  1:645517/139‑1 (MQ=255)
                                                           cTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                                                                      >  2:582433/1‑139 (MQ=255)
                                                           cTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                                                                      <  2:1225779/139‑1 (MQ=255)
                                                              gATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATACCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACCAAAGAACGTAAACAAATCCATCCagaag                                                                   <  1:33604/139‑1 (MQ=255)
                                                                ttATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGat                                                                 <  1:993776/139‑1 (MQ=255)
                                                                ttATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGat                                                                 <  2:1551895/139‑1 (MQ=255)
                                                                   tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                                   <  1:1008776/70‑1 (MQ=255)
                                                                   tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                                   >  2:1008776/1‑70 (MQ=255)
                                                                   tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    <  2:925742/85‑1 (MQ=255)
                                                                   tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    >  1:925742/1‑85 (MQ=255)
                                                                   tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        >  2:1542867/1‑129 (MQ=255)
                                                                   tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  1:1542867/129‑1 (MQ=255)
                                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGt                                                                                                                              >  2:322714/1‑70 (MQ=255)
                                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGt                                                                                                                              <  1:322714/70‑1 (MQ=255)
                                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    <  1:266229/80‑1 (MQ=255)
                                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    >  2:266229/1‑80 (MQ=255)
                                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                      >  1:1149482/1‑94 (MQ=255)
                                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                      <  2:1149482/94‑1 (MQ=255)
                                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                             >  1:786818/1‑103 (MQ=255)
                                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                             <  2:786818/103‑1 (MQ=255)
                                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  2:385382/124‑1 (MQ=255)
                                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        >  1:385382/1‑124 (MQ=255)
                                                                         tgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        >  1:296038/1‑123 (MQ=255)
                                                                         tgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  2:296038/123‑1 (MQ=255)
                                                                           ctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCAc                                                      >  1:358077/1‑139 (MQ=255)
                                                                                ataCGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaa                                                 >  2:1244299/1‑139 (MQ=255)
                                                                                     aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                             <  2:1418378/90‑1 (MQ=255)
                                                                                     aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                             >  1:1418378/1‑90 (MQ=255)
                                                                                     aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGAc                                            >  1:1227255/1‑139 (MQ=255)
                                                                                        gAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTc                                         >  1:1003491/1‑139 (MQ=255)
                                                                                           gATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGaa                                      <  2:1003491/139‑1 (MQ=255)
                                                                                               aGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                                     >  1:427811/1‑136 (MQ=255)
                                                                                               aGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                                     <  2:427811/136‑1 (MQ=255)
                                                                                                gTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        >  1:978783/1‑100 (MQ=255)
                                                                                                gTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  2:978783/100‑1 (MQ=255)
                                                                                                      atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGAtt                                                                                                                       <  2:899586/47‑1 (MQ=255)
                                                                                                      atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGAtt                                                                                                                       >  1:899586/1‑47 (MQ=255)
                                                                                                      atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcat                                                                                                                   >  1:1200540/1‑51 (MQ=255)
                                                                                                      atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcat                                                                                                                   <  2:1200540/51‑1 (MQ=255)
                                                                                                      atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  2:997178/94‑1 (MQ=255)
                                                                                                      atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        >  1:997178/1‑94 (MQ=255)
                                                                                                          aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       <  1:1244299/139‑1 (MQ=255)
                                                                                                          aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       <  1:441386/139‑1 (MQ=255)
                                                                                                          aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       <  1:358954/139‑1 (MQ=255)
                                                                                                              ggtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACAtct                   >  2:782474/1‑139 (MQ=255)
                                                                                                               gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                                   >  1:1015197/1‑90 (MQ=255)
                                                                                                               gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                                   <  2:1015197/90‑1 (MQ=255)
                                                                                                               gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                     <  1:782078/88‑1 (MQ=255)
                                                                                                               gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                     >  2:782078/1‑88 (MQ=255)
                                                                                                               gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaa                                                >  2:1462534/1‑109 (MQ=255)
                                                                                                               gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaa                                                <  1:1462534/109‑1 (MQ=255)
                                                                                                               gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACAtctc                  >  1:1133037/1‑139 (MQ=255)
                                                                                                                tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaa                                                <  1:1201133/108‑1 (MQ=255)
                                                                                                                tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaa                                                >  2:1201133/1‑108 (MQ=255)
                                                                                                                tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGAt                          >  2:1441914/1‑130 (MQ=255)
                                                                                                                tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGAt                          <  1:1441914/130‑1 (MQ=255)
                                                                                                                tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATctct                 >  1:1449329/1‑139 (MQ=255)
                                                                                                                 ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTc                                                                                                                      <  1:123551/37‑1 (MQ=255)
                                                                                                                 ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTc                                                                                                                      >  2:123551/1‑37 (MQ=255)
                                                                                                                 ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                        >  2:441461/1‑67 (MQ=255)
                                                                                                                 ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                        <  1:441461/67‑1 (MQ=255)
                                                                                                                 ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTa                >  2:518700/1‑139 (MQ=255)
                                                                                                                     tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATa                         >  2:71279/1‑126 (MQ=255)
                                                                                                                     tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATa                         <  1:71279/126‑1 (MQ=255)
                                                                                                                          ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                              >  1:1427361/1‑68 (MQ=255)
                                                                                                                          ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                              <  2:1427361/68‑1 (MQ=255)
                                                                                                                          ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAt       >  2:566189/1‑139 (MQ=255)
                                                                                                                          ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAt       >  1:1440691/1‑139 (MQ=255)
                                                                                                                               ttAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                       <  1:1512012/70‑1 (MQ=255)
                                                                                                                               ttAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                       >  2:1512012/1‑70 (MQ=255)
                                                                                                                               ttAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGGaa  <  2:1133037/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
ATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGGAA  >  NC_000913/2787304‑2787569

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