Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 3,815,859 Δ82 bp 100% [rph][rph] [rph], [rph]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3815859 3815940 82 89 [0] [0] 86 [rph]–[rph] [rph], [rph]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 38158580 (0.000)69 (0.890) 48/158 NT 100% pseudogene (24/48 nt) rph ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
?NC_000913 3815941 = 0 (0.000)pseudogene (609/669 nt) rph ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

ATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTC  >  NC_000913/3815759‑3815955
                                                                                                   |                                                                                                 
atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                                   <  1:1161489‑M1/100‑1 (MQ=255)
atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                                   >  2:293878‑M1/1‑100 (MQ=255)
 tGGTTTCATGCCTTCGCTGCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                  <  1:582459‑M1/100‑2 (MQ=255)
  ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                                                                 <  1:106192‑M1/100‑3 (MQ=255)
   gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                                                >  2:44800‑M1/1‑97 (MQ=255)
     ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                              <  2:772256‑M1/100‑6 (MQ=255)
      tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                             >  2:2053388‑M1/1‑94 (MQ=255)
      ncaTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                              >  2:1904191‑M1/2‑94 (MQ=255)
        aTGCGTTCGCGCCTCTTCTTACTTTAAAAAAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGccgtcccc                                                                                            <  2:1945412/99‑8 (MQ=255)
               cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                                                                    <  1:903502‑M1/100‑16 (MQ=255)
               cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                                                                    <  1:181577‑M1/100‑16 (MQ=255)
                gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttca                                                                                   >  2:281020‑M1/1‑84 (MQ=255)
                gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttca                                                                                   >  2:281041‑M1/1‑84 (MQ=255)
                 ctcctcATCTTACTTTTATAAAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcact                                                                                      <  2:613129‑M1/96‑14 (MQ=255)
                  tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCTCCTTCTGCcgtcccctgcacttcaat                                                                                 <  1:1436632‑M1/100‑19 (MQ=255)
                  tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                 <  1:810994‑M1/100‑19 (MQ=255)
                   cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                >  1:606507‑M1/1‑81 (MQ=255)
                   cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                >  1:840775‑M1/1‑81 (MQ=255)
                   cctcATCTTACTTTAAAACAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                                                            <  2:1207009‑M1/88‑8 (MQ=255)
                    ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatga                                                                               >  1:1737360‑M1/1‑80 (MQ=255)
                     tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                              >  1:422924‑M1/1‑79 (MQ=255)
                      cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                             <  1:1749541‑M1/100‑23 (MQ=255)
                       aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgc                                                                            <  1:1704273‑M1/100‑24 (MQ=255)
                         cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                          >  1:2194443‑M1/1‑75 (MQ=255)
                          ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                         <  1:1372860‑M1/100‑27 (MQ=255)
                          ttACTTTTCTACAGACAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                         <  2:985788‑M1/100‑27 (MQ=255)
                             cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                                                                      >  1:855412‑M1/1‑71 (MQ=255)
                                  cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                 >  1:1009516‑M1/1‑66 (MQ=255)
                                    aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                                                               >  1:5990‑M1/1‑64 (MQ=255)
                                       gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                            >  1:1802396‑M1/1‑61 (MQ=255)
                                        aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                           <  1:122714‑M1/100‑41 (MQ=255)
                                        aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                           >  1:2015449‑M1/1‑60 (MQ=255)
                                        aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                           >  1:2015450‑M1/1‑60 (MQ=255)
                                         cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                          >  1:1187022‑M1/1‑59 (MQ=255)
                                         cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                          >  2:605829‑M1/1‑59 (MQ=255)
                                         cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                          >  2:323149‑M1/1‑59 (MQ=255)
                                          naaaaaaagGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                                                         >  1:278442‑M1/2‑58 (MQ=255)
                                            aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                                       >  2:882594‑M1/1‑56 (MQ=255)
                                             aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactac                                                      <  1:1090683‑M1/100‑46 (MQ=255)
                                               aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                    <  2:362259‑M1/100‑48 (MQ=255)
                                               aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                    <  1:1347577‑M1/100‑48 (MQ=255)
                                                aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                   <  2:883290‑M1/100‑49 (MQ=255)
                                                 aGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                  >  1:1984254‑M1/1‑51 (MQ=255)
                                                 aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                  <  1:186841‑M1/100‑50 (MQ=255)
                                                  ggCGACTCATCAGTCGCCTTAATAATCAGTTTGCCAGCGCCGCCTTCAGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttca                                                 <  2:1587416‑M1/100‑51 (MQ=255)
                                                  ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttga                                                 >  2:1431013‑M1/1‑50 (MQ=255)
                                                    cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                               >  2:1780911‑M1/1‑48 (MQ=255)
                                                    cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                               >  1:2215513‑M1/1‑48 (MQ=255)
                                                    cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                               >  2:371970‑M1/1‑48 (MQ=255)
                                                     gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt                                              <  2:897290‑M1/100‑54 (MQ=255)
                                                     gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt                                              >  2:2238949‑M1/1‑47 (MQ=255)
                                                      aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtc                                             <  1:845987‑M1/100‑55 (MQ=255)
                                                        tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                           <  2:1272706‑M1/100‑57 (MQ=255)
                                                        tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                           <  2:133095‑M1/100‑57 (MQ=255)
                                                        tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                           <  1:955537‑M1/100‑57 (MQ=255)
                                                        tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                           <  1:941331‑M1/100‑57 (MQ=255)
                                                        tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                           <  1:167241‑M1/100‑57 (MQ=255)
                                                         catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcactacgtccatgtcggt                                          >  2:577035‑M1/1‑43 (MQ=255)
                                                         catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt                                          >  1:668013‑M1/1‑43 (MQ=255)
                                                            caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                       >  1:2208456‑M1/1‑40 (MQ=255)
                                                            caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                       <  2:1053146‑M1/100‑61 (MQ=255)
                                                                cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcg                                   <  2:1918459‑M1/100‑65 (MQ=255)
                                                                 gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcg                                   >  2:1627196‑M1/1‑35 (MQ=255)
                                                                    ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcag                               >  2:658044‑M1/1‑32 (MQ=255)
                                                                    ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcag                               >  1:264582‑M1/1‑32 (MQ=255)
                                                                     tAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                              >  2:1617911‑M1/1‑31 (MQ=255)
                                                                         aaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                          <  1:1907542‑M1/100‑74 (MQ=255)
                                                                           tCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttc                        <  2:2222412‑M1/100‑76 (MQ=255)
                                                                           tCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttc                        <  1:934332‑M1/100‑76 (MQ=255)
                                                                           tCAGTTTGCCAGCGCCGCCTTCTGCcgtccccagcactttaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttc                        <  2:1649054‑M1/100‑76 (MQ=255)
                                                                            cAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttca                       >  2:1773924‑M1/1‑24 (MQ=255)
                                                                            cAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttca                       >  2:1332453‑M1/1‑24 (MQ=255)
                                                                              gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                     >  2:198236‑M1/1‑22 (MQ=255)
                                                                               tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                    <  1:1028674‑M1/100‑80 (MQ=255)
                                                                                  gCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtat                 <  1:436997‑M1/100‑83 (MQ=255)
                                                                                    cAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattc               <  1:2107748‑M1/100‑85 (MQ=255)
                                                                                    cAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattc               >  1:284019‑M1/1‑16 (MQ=255)
                                                                                        gccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccaga           >  2:169046‑M1/1‑12 (MQ=255)
                                                                                         ccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat          <  1:177637‑M1/100‑90 (MQ=255)
                                                                                         ccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat          >  1:313552‑M1/1‑11 (MQ=255)
                                                                                          cgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc         <  1:868447‑M1/100‑91 (MQ=255)
                                                                                             cTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtgtcttcggcagagtcttcaacggcttccagatcgca      >  1:354205‑M1/1‑7 (MQ=255)
                                                                                             cTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgca      >  1:1526567‑M1/1‑7 (MQ=255)
                                                                                             cTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtatgccagatcgca      >  2:1539672‑M1/1‑7 (MQ=255)
                                                                                               tCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaa    >  1:137200‑M1/1‑5 (MQ=255)
                                                                                                cTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaac   >  1:1115353‑M1/1‑4 (MQ=255)
                                                                                                cTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaac   <  2:797986‑M1/100‑97 (MQ=255)
                                                                                                cTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaac   >  2:803591‑M1/1‑4 (MQ=255)
                                                                                                 tGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaacc  >  1:1218100‑M1/1‑3 (MQ=255)
                                                                                                   |                                                                                                 
ATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTC  >  NC_000913/3815759‑3815955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: