Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,098 T→A 9.0% L665Q (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,0980TA 9.0% 227.5 / 11.9 67L665Q (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.01e-01

CCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCC  >  NC_000913/163952‑164244
                                                                                                                                                  |                                                                                                                                                  
ccGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGgctg                                                                                                                                                   >  2:914785/1‑148 (MQ=255)
ccGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGgctg                                                                                                                                                   >  4:642899/1‑148 (MQ=255)
   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                                                                                         <  4:805803/120‑1 (MQ=255)
                                                                                      cGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATc                                                                                         >  3:805803/1‑120 (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                                            >  1:514376/1‑136 (MQ=255)
                                                                                                                   gaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTgg                                            <  2:514376/136‑1 (MQ=255)
                                                                                                                    aaaGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAg                                        <  4:791841/139‑1 (MQ=255)
                                                                                                                    aaaGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAg                                        >  3:791841/1‑139 (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                                                        >  4:40565/1‑102 (MQ=255)
                                                                                                                                         aaaCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgc                                                        <  3:40565/102‑1 (MQ=255)
                                                                                                                                         aaaCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGAt                                          >  3:143020/1‑116 (MQ=255)
                                                                                                                                         aaaCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGAt                                          <  4:143020/116‑1 (MQ=255)
                                                                                                                                          tatgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTgg                                                               >  1:99970/4‑94 (MQ=255)
                                                                                                                                          tatgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTgg                                                               <  2:99970/91‑1 (MQ=255)
                                                                                                                                           atgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgc                                                        >  4:355176/3‑100 (MQ=255)
                                                                                                                                           atgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgc                                                        <  3:355176/98‑1 (MQ=255)
                                                                                                                                           aCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                       >  1:460019/1‑101 (MQ=255)
                                                                                                                                           aCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                       <  2:460019/101‑1 (MQ=255)
                                                                                                                                            tgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTgg                                                               >  2:490659/2‑92 (MQ=255)
                                                                                                                                            tgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTgg                                                               <  1:490659/91‑1 (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  >  4:283766/1‑148 (MQ=255)
                                                                                                                                                 ctgctgCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGcc  >  2:85883/1‑148 (MQ=255)
                                                                                                                                                 ctgctgCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTAACAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAAt                                   >  2:1035812/1‑115 (MQ=255)
                                                                                                                                                 ctgctgCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTAACAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAAt                                   <  1:1035812/115‑1 (MQ=255)
                                                                                                                                                  tgctgcCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGa                                                                                                 <  3:21230/52‑1 (MQ=255)
                                                                                                                                                  tgctgcCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGa                                                                                                 >  4:21230/1‑52 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCC  >  NC_000913/163952‑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: