Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,535,441 C→T pseudogene (439/1029 nt) yjhR → pseudogene, helicase family, putative frameshift suppressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,535,4410CT100.0% 317.4 / NA 97pseudogene (439/1029 nt)yjhRpseudogene, helicase family, putative frameshift suppressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (47/50);  total (47/50)

TCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGCCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTGTGG  >  NC_000913/4535314‑4535573
                                                                                                                               |                                                                                                                                    
tcCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCAtttt                                                                                                                             <  1:784777/137‑1 (MQ=255)
tcCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCAtttt                                                                                                                             >  2:784777/1‑137 (MQ=255)
   tGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTc                                                                                                                            >  2:163569/1‑135 (MQ=255)
   tGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTc                                                                                                                            <  1:163569/135‑1 (MQ=255)
       tACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAg                                                                                                                    >  1:1081362/1‑139 (MQ=255)
              aaCATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTAtt                                                                                                             >  2:138185/1‑139 (MQ=255)
               aCATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTg                                                                                                            >  1:547919/1‑139 (MQ=255)
                   gAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                                      >  2:1301471/1‑125 (MQ=255)
                   gAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                                      <  1:1301471/125‑1 (MQ=255)
                   gAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCaa                                                                                                        <  2:1081362/139‑1 (MQ=255)
                    aaGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCaaa                                                                                                       <  1:1057412/139‑1 (MQ=255)
                       aTGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCGGTCTCACGCGCTAAAGATAGTTTCTGGTATTTGGCGATATGGCCCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                                   <  1:1358670/139‑1 (MQ=255)
                        tGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGa                                                                                                                     >  1:594522/1‑121 (MQ=255)
                        tGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGa                                                                                                                     <  2:594522/121‑1 (MQ=255)
                        tGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                                   >  1:554173/1‑139 (MQ=255)
                        tGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                                   <  2:267670/139‑1 (MQ=255)
                        tGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                                   >  2:26650/1‑139 (MQ=255)
                         ggCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATAtc                                                                                                  >  1:1007557/1‑139 (MQ=255)
                               tttCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                            <  2:343107/139‑1 (MQ=255)
                               tttCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                            <  2:554173/139‑1 (MQ=255)
                                 tCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttct                                                                                          <  1:138185/139‑1 (MQ=255)
                                 tCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttct                                                                                          <  2:496599/139‑1 (MQ=255)
                                   ttGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                                      <  1:1502865/109‑1 (MQ=255)
                                   ttGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                                      >  2:1502865/1‑109 (MQ=255)
                                    tGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCa                                                                                       <  1:544931/139‑1 (MQ=255)
                                          cAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTcgccg                                                                                                                        <  2:1314208/100‑1 (MQ=255)
                                          cAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTcgccg                                                                                                                        >  1:1314208/1‑100 (MQ=255)
                                            aCAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGAtaataa                                                                               >  2:1557728/1‑139 (MQ=255)
                                                cACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttcttc                                                                                        >  1:1445216/1‑126 (MQ=255)
                                                cACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttcttc                                                                                        <  2:1445216/126‑1 (MQ=255)
                                                   tATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTg                                                                        <  2:42230/139‑1 (MQ=255)
                                                          aaTGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                            >  1:811551/1‑112 (MQ=255)
                                                          aaTGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                            <  2:811551/112‑1 (MQ=255)
                                                          aaTGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGataata                                                                                >  2:949243/1‑124 (MQ=255)
                                                          aaTGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGataata                                                                                <  1:949243/124‑1 (MQ=255)
                                                          aaTGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTg                                                                 >  1:1339307/1‑139 (MQ=255)
                                                               tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                                      >  2:15730/1‑81 (MQ=255)
                                                               tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                                      <  1:15730/81‑1 (MQ=255)
                                                               tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGataata                                                                                <  1:279019/119‑1 (MQ=255)
                                                               tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGataata                                                                                >  2:279019/1‑119 (MQ=255)
                                                               tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGttt                                                            >  1:1239168/1‑139 (MQ=255)
                                                               tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGttt                                                            >  1:856280/1‑139 (MQ=255)
                                                                 cTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCa                                                          >  2:967813/1‑139 (MQ=255)
                                                                    tctcACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtt                                                                                              <  2:169606/100‑1 (MQ=255)
                                                                    tctcACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtt                                                                                              >  1:169606/1‑100 (MQ=255)
                                                                        aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGa                                                                                                                     >  1:1194038/1‑73 (MQ=255)
                                                                        aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGa                                                                                                                     <  2:1194038/73‑1 (MQ=255)
                                                                        aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGa                                                   <  1:1557728/139‑1 (MQ=255)
                                                                              tAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGa                                                                                                                     <  1:818652/67‑1 (MQ=255)
                                                                              tAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGa                                                                                                                     >  2:818652/1‑67 (MQ=255)
                                                                                      gTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTAt                                                                                                              <  1:719498/66‑1 (MQ=255)
                                                                                      gTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTAt                                                                                                              >  2:719498/1‑66 (MQ=255)
                                                                                      gTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTg                                                                        <  1:118000/104‑1 (MQ=255)
                                                                                      gTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTg                                                                        >  2:118000/1‑104 (MQ=255)
                                                                                       ttttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttcttc                                                                                        >  2:66946/1‑87 (MQ=255)
                                                                                       ttttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttcttc                                                                                        <  1:66946/87‑1 (MQ=255)
                                                                                        tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCAtttt                                                                                                                             <  1:1539069/49‑1 (MQ=255)
                                                                                        tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCAtttt                                                                                                                             >  2:1539069/1‑49 (MQ=255)
                                                                                        tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCg                                   >  2:1228403/1‑139 (MQ=255)
                                                                                            tGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGcaca                               <  1:1414847/139‑1 (MQ=255)
                                                                                            tGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGcaca                               <  2:1583114/139‑1 (MQ=255)
                                                                                             ggTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACAt                              >  2:883709/1‑139 (MQ=255)
                                                                                                aTTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATaca                           <  1:26650/139‑1 (MQ=255)
                                                                                                 tttGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGa                                                                                                                     <  2:89984/48‑1 (MQ=255)
                                                                                                 tttGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGa                                                                                                                     >  1:89984/1‑48 (MQ=255)
                                                                                                   tGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGa                                                                                                                     <  1:1120626/46‑1 (MQ=255)
                                                                                                   tGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGa                                                                                                                     >  2:1120626/1‑46 (MQ=255)
                                                                                                   tGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACaa                        >  1:1344381/1‑139 (MQ=255)
                                                                                                    ggCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                                      <  1:51984/44‑1 (MQ=255)
                                                                                                    ggCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                                      >  2:51984/1‑44 (MQ=255)
                                                                                                    ggCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                            >  1:71129/1‑70 (MQ=255)
                                                                                                    ggCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                            <  2:71129/70‑1 (MQ=255)
                                                                                                      cGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAAtt                     >  2:621435/1‑139 (MQ=255)
                                                                                                        atatGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATAcaca                         <  1:1277287/133‑1 (MQ=255)
                                                                                                        atatGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATAcaca                         >  2:1277287/1‑133 (MQ=255)
                                                                                                        atatGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTc                   >  2:1169977/1‑139 (MQ=255)
                                                                                                         tatGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCg                                   >  2:196048/1‑122 (MQ=255)
                                                                                                         tatGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCg                                   <  1:196048/122‑1 (MQ=255)
                                                                                                          atGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCaa                 <  2:547919/139‑1 (MQ=255)
                                                                                                            ggACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACg               <  1:1169977/139‑1 (MQ=255)
                                                                                                            ggACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACg               <  2:1344381/139‑1 (MQ=255)
                                                                                                               ccTTATCGAGATGGAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGa                                                   >  1:1077670/1‑100 (MQ=255)
                                                                                                               ccTTATCGAGATGCAGTCCGCATTTACGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATCATGCGTTGCAATTTGAGTTTCAGAAACGa                                                   <  2:1077670/100‑1 (MQ=255)
                                                                                                                 ttATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa          <  1:967813/139‑1 (MQ=255)
                                                                                                                 ttATCGAGATGCAGTCCGCATTCTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa          <  2:1158782/139‑1 (MQ=255)
                                                                                                                  tATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCAt         >  1:1414949/1‑139 (MQ=255)
                                                                                                                   aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTgg                                                                                                           >  2:621593/1‑40 (MQ=255)
                                                                                                                   aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTgg                                                                                                           <  1:621593/40‑1 (MQ=255)
                                                                                                                   aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa          >  2:1090543/1‑137 (MQ=255)
                                                                                                                   aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa          <  1:1090543/137‑1 (MQ=255)
                                                                                                                   aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATg        <  2:856280/139‑1 (MQ=255)
                                                                                                                      gagaTGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGgtg     <  1:621435/139‑1 (MQ=255)
                                                                                                                        gaTGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                                   >  1:1504435/1‑43 (MQ=255)
                                                                                                                        gaTGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                                   <  2:1504435/43‑1 (MQ=255)
                                                                                                                         aTGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTGTgg  >  2:1587834/1‑139 (MQ=255)
                                                                                                                          tGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATAcaca                         <  2:61118/115‑1 (MQ=255)
                                                                                                                          tGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATAcaca                         >  1:61118/1‑115 (MQ=255)
                                                                                                                               |                                                                                                                                    
TCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGCCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTGTGG  >  NC_000913/4535314‑4535573

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