Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,101 T→A 9.0% L666Q (CTG→CAG)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,1010TA 9.0% 230.1 / 11.8 67L666Q (CTG→CAG) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base A (3/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

TGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCC  >  NC_000913/163955‑164244
                                                                                                                                                  |                                                                                                                                               
tGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGctgctg                                                                                                                                                <  3:642899/148‑1 (MQ=255)
              gCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTgg                                                                                                                                  <  1:437783/148‑1 (MQ=255)
                   cTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTAc                                                                                                                             >  2:385063/1‑148 (MQ=255)
                     cAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACAt                                                                                                                           <  4:923448/148‑1 (MQ=255)
                        cTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCa                                                                                                                        <  1:385063/148‑1 (MQ=255)
                            aCGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTAc                                                                                                                    >  2:451443/1‑148 (MQ=255)
                                  gaagcggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGcc                                                                                                              <  1:1004908/148‑1 (MQ=255)
                                          aCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCa                                                                                                      <  1:451443/148‑1 (MQ=255)
                                                aCGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGAc                                                                                                >  3:947425/1‑148 (MQ=255)
                                                     tGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTa                                                                                           <  4:947425/148‑1 (MQ=255)
                                                     tGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTa                                                                                           >  4:303389/1‑148 (MQ=255)
                                                      gCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTAt                                                                                          <  3:303389/148‑1 (MQ=255)
                                                           tGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGc                                                                                     <  1:914785/148‑1 (MQ=255)
                                                                            tGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTa                                                                        >  4:627201/1‑144 (MQ=255)
                                                                            tGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTa                                                                        <  3:627201/144‑1 (MQ=255)
                                                                                   cGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATc                                                                                         >  3:805803/1‑120 (MQ=255)
                                                                                   cGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATc                                                                                         <  4:805803/120‑1 (MQ=255)
                                                                                    ggAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGaa                                                                                                          <  2:753194/102‑1 (MQ=255)
                                                                                    ggAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGaa                                                                                                          >  1:753194/1‑102 (MQ=255)
                                                                                                 ccAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACAtt                                                                                             >  2:426689/1‑102 (MQ=255)
                                                                                                 ccAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACAtt                                                                                             <  1:426689/102‑1 (MQ=255)
                                                                                                  cAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACata                                                                                                                                         <  4:810085/57‑1 (MQ=255)
                                                                                                  cAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACata                                                                                                                                         >  3:810085/1‑57 (MQ=255)
                                                                                                   aGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTg                                             >  3:46088/1‑148 (MQ=255)
                                                                                                               tgaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAg                                        >  3:579758/1‑141 (MQ=255)
                                                                                                               tgaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAg                                        <  4:579758/141‑1 (MQ=255)
                                                                                                               tgaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTg                                 <  4:46088/148‑1 (MQ=255)
                                                                                                                gaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTgg                                            <  2:514376/136‑1 (MQ=255)
                                                                                                                gaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTgg                                            >  1:514376/1‑136 (MQ=255)
                                                                                                                 aaaGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAg                                        >  3:791841/1‑139 (MQ=255)
                                                                                                                 aaaGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAg                                        <  4:791841/139‑1 (MQ=255)
                                                                                                                     tttATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                           >  1:561402/1‑148 (MQ=255)
                                                                                                                      ttATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACt                                                                      <  2:704701/104‑1 (MQ=255)
                                                                                                                      ttATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACt                                                                      >  1:704701/1‑104 (MQ=255)
                                                                                                                      ttATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgc                                                        <  2:672711/118‑1 (MQ=255)
                                                                                                                      ttATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgc                                                        >  1:672711/1‑118 (MQ=255)
                                                                                                                       tATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTTGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTAt                                                                                          >  4:994920/1‑83 (MQ=255)
                                                                                                                       tATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTAt                                                                                          <  3:994920/83‑1 (MQ=255)
                                                                                                                        aTTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGt                                       >  2:1059775/1‑133 (MQ=255)
                                                                                                                        aTTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGt                                       <  1:1059775/133‑1 (MQ=255)
                                                                                                                          tGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATg                                                         >  3:525649/1‑113 (MQ=255)
                                                                                                                          tGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATg                                                         <  4:525649/113‑1 (MQ=255)
                                                                                                                            gCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                       >  4:754717/1‑113 (MQ=255)
                                                                                                                            gCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                       <  3:754717/113‑1 (MQ=255)
                                                                                                                            gCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTaca                    <  2:561402/148‑1 (MQ=255)
                                                                                                                                cgcTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGa                                                                   >  2:162661/1‑97 (MQ=255)
                                                                                                                                cgcTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGa                                                                   <  1:162661/97‑1 (MQ=255)
                                                                                                                                    ggAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTgg                                                               >  3:597639/1‑97 (MQ=255)
                                                                                                                                    ggAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTgg                                                               <  4:597639/97‑1 (MQ=255)
                                                                                                                                    ggAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGa                                     <  4:411961/123‑1 (MQ=255)
                                                                                                                                    ggAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGa                                     >  3:411961/1‑123 (MQ=255)
                                                                                                                                      aaaCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTAcgcg                                                                             >  4:71101/1‑81 (MQ=255)
                                                                                                                                      aaaCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTAcgcg                                                                             <  3:71101/81‑1 (MQ=255)
                                                                                                                                      aaaCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgc                                                        <  3:40565/102‑1 (MQ=255)
                                                                                                                                      aaaCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgc                                                        >  4:40565/1‑102 (MQ=255)
                                                                                                                                      aaaCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGAt                                          <  4:143020/116‑1 (MQ=255)
                                                                                                                                      aaaCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGAt                                          >  3:143020/1‑116 (MQ=255)
                                                                                                                                       tatgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTgg                                                               >  1:99970/4‑94 (MQ=255)
                                                                                                                                       tatgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTgg                                                               <  2:99970/91‑1 (MQ=255)
                                                                                                                                        atgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgc                                                        <  3:355176/98‑1 (MQ=255)
                                                                                                                                        atgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgc                                                        >  4:355176/3‑100 (MQ=255)
                                                                                                                                        aCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                       <  2:460019/101‑1 (MQ=255)
                                                                                                                                        aCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                       >  1:460019/1‑101 (MQ=255)
                                                                                                                                         tgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTgg                                                               <  1:490659/91‑1 (MQ=255)
                                                                                                                                         tgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTgg                                                               >  2:490659/2‑92 (MQ=255)
                                                                                                                                           tGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                       <  3:243422/98‑1 (MQ=255)
                                                                                                                                           tGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                       >  4:243422/1‑98 (MQ=255)
                                                                                                                                           tGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGaa     >  4:218628/1‑148 (MQ=255)
                                                                                                                                             gctgctTCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGc   >  4:340302/1‑148 (MQ=255)
                                                                                                                                             gctgctGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGc   >  4:905753/1‑148 (MQ=255)
                                                                                                                                              ctgctgCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGcc  >  2:85883/1‑148 (MQ=255)
                                                                                                                                              ctgctgCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGcc  >  4:283766/1‑148 (MQ=255)
                                                                                                                                              ctgctgCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTAACAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAAt                                   <  1:1035812/115‑1 (MQ=255)
                                                                                                                                              ctgctgCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTAACAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAAt                                   >  2:1035812/1‑115 (MQ=255)
                                                                                                                                               tgctgcCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGa                                                                                                 <  3:21230/52‑1 (MQ=255)
                                                                                                                                               tgctgcCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGa                                                                                                 >  4:21230/1‑52 (MQ=255)
                                                                                                                                                gctgcCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGAt                                                                  <  3:1035556/82‑1 (MQ=255)
                                                                                                                                                gctgcCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGAt                                                                  >  4:1035556/1‑82 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCC  >  NC_000913/163955‑164244

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