Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,917,910 T→C 12.3% intergenic (+97/+135) barA → / ← gudD hybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,917,9100TC12.3% 207.7 / 18.4 75intergenic (+97/+135)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (5/4);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCT  >  NC_000913/2917762‑2918045
                                                                                                                                                    |                                                                                                                                       
gTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                                                         <  4:1334275/149‑1 (MQ=255)
 ttGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                                                        <  2:1474996/149‑1 (MQ=255)
    gACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                                                     >  2:611379/1‑149 (MQ=255)
    gACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                                                     >  3:447830/1‑149 (MQ=255)
    gACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                                                     <  4:1298081/149‑1 (MQ=255)
            ggATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                                                         <  3:357166/137‑1 (MQ=255)
            ggATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                                                         >  4:357166/1‑137 (MQ=255)
            ggATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTCGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                                                             >  4:1511129/1‑149 (MQ=255)
               tAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                                                                          <  1:955529/149‑1 (MQ=255)
                            gAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCa                                                                                                             >  4:353204/1‑149 (MQ=255)
                                 cAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGaa                                                                                                        <  1:611379/149‑1 (MQ=255)
                                    cAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTaaa                                                                                                                 <  4:1470098/137‑1 (MQ=255)
                                    cAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTaaa                                                                                                                 >  3:1470098/1‑137 (MQ=255)
                                       aaTTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTaaaa                                                                                                  <  4:447830/149‑1 (MQ=255)
                                       aaTTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTaaaa                                                                                                  >  3:1335353/1‑149 (MQ=255)
                                        aTTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTaaaaa                                                                                                 >  1:1327072/1‑149 (MQ=255)
                                             cGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatgat                                                                                            <  2:1327072/149‑1 (MQ=255)
                                              gggTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatg                                                                                              >  4:688903/1‑146 (MQ=255)
                                              gggTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatg                                                                                              <  3:688903/146‑1 (MQ=255)
                                                gTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                                               <  3:747731/111‑1 (MQ=255)
                                                gTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                                               >  2:1048348/1‑111 (MQ=255)
                                                gTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                                               >  4:747731/1‑111 (MQ=255)
                                                gTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                                               <  1:1048348/111‑1 (MQ=255)
                                                gTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTa                                                                                         >  3:1518954/1‑149 (MQ=255)
                                                gTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTa                                                                                         <  3:1511129/149‑1 (MQ=255)
                                                   aTGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGAc                                                                                      <  4:1335353/149‑1 (MQ=255)
                                                            ttGTCTAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAAt                                                                                                                     >  1:1341843/1‑109 (MQ=255)
                                                            ttGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAAt                                                                                                                     <  2:1341843/109‑1 (MQ=255)
                                                             tGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATa                                                                            <  3:353204/149‑1 (MQ=255)
                                                                 gAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCaa                                                                                                            <  4:1861988/113‑1 (MQ=255)
                                                                 gAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCaa                                                                                                            >  3:1861988/1‑113 (MQ=255)
                                                                     aGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACAt                                                                             <  1:1187149/140‑1 (MQ=255)
                                                                     aGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACAt                                                                             >  2:1187149/1‑140 (MQ=255)
                                                                       cAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAgg                                                                  <  4:1518954/149‑1 (MQ=255)
                                                                              cTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGa                                                                                  <  4:458191/126‑1 (MQ=255)
                                                                              cTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGa                                                                                  >  3:458191/1‑126 (MQ=255)
                                                                               tttttATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCaa                                                                                                            >  1:1759187/1‑99 (MQ=255)
                                                                               tttttATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCaa                                                                                                            <  2:1759187/99‑1 (MQ=255)
                                                                               tttttATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTg                                                          >  1:85389/1‑149 (MQ=255)
                                                                                       aCGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATa                                                                            >  4:209221/1‑123 (MQ=255)
                                                                                       aCGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATa                                                                            <  3:129808/123‑1 (MQ=255)
                                                                                       aCGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATa                                                                            >  4:129808/1‑123 (MQ=255)
                                                                                       aCGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATa                                                                            <  3:209221/123‑1 (MQ=255)
                                                                                                   gATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCaaa                                      <  2:85389/149‑1 (MQ=255)
                                                                                                          gACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCacaaca                                                                    >  4:1360702/1‑112 (MQ=255)
                                                                                                          gACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCacaaca                                                                    <  3:1360702/112‑1 (MQ=255)
                                                                                                             aaaGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATa                                                                            >  3:1001084/1‑101 (MQ=255)
                                                                                                             aaaGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATa                                                                            <  4:1001084/101‑1 (MQ=255)
                                                                                                             aaaGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGtct                                          >  2:1334342/1‑135 (MQ=255)
                                                                                                             aaaGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGtct                                          <  1:1334342/135‑1 (MQ=255)
                                                                                                                    cTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAGCCAATTAAATTCaa                                                                                                            <  1:271747/62‑1 (MQ=255)
                                                                                                                    cTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCaa                                                                                                            >  2:271747/1‑62 (MQ=255)
                                                                                                                     tgtgTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGa                                                                                       >  1:1716513/1‑82 (MQ=255)
                                                                                                                     tgtgTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGa                                                                                       <  2:1716513/82‑1 (MQ=255)
                                                                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCa                                                 <  3:200953/116‑1 (MQ=255)
                                                                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCa                                                 >  2:859220/1‑116 (MQ=255)
                                                                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCa                                                 >  4:200953/1‑116 (MQ=255)
                                                                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCa                                                 <  1:859220/116‑1 (MQ=255)
                                                                                                                              gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCg                                                        <  4:1866329/104‑1 (MQ=38)
                                                                                                                              gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCg                                                        >  3:1866329/1‑104 (MQ=38)
                                                                                                                                  gataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTTAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGtt                                                      >  1:1874127/6‑102 (MQ=39)
                                                                                                                                  gataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGtt                                                      >  1:490609/6‑102 (MQ=39)
                                                                                                                                  gataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGtt                                                      <  2:490609/97‑1 (MQ=39)
                                                                                                                                  gataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGtt                                                      <  2:1874127/97‑1 (MQ=39)
                                                                                                                                  gataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc       >  3:797844/6‑149 (MQ=255)
                                                                                                                                   ataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTACATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACATCAGGCTATTTTGCGTTTAGCATCAGTCTCAAAACGGCTCCAGATAGAGCCgg                  <  4:1605401/133‑1 (MQ=255)
                                                                                                                                   ataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCgg                  >  3:1605401/5‑137 (MQ=255)
                                                                                                                                       ggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCt  >  1:1327623/1‑149 (MQ=255)
                                                                                                                                       ggcgTTAACGCCTAATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGAt                                                                                 >  2:1757381/1‑70 (MQ=255)
                                                                                                                                       ggcgTTAACGCCTAATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGAt                                                                                 <  1:1757381/70‑1 (MQ=255)
                                                                                                                                        gcgTTCGCGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTa                                                    >  2:720083/1‑98 (MQ=255)
                                                                                                                                        gcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTa                                                    <  1:720083/98‑1 (MQ=255)
                                                                                                                                             aaCGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTAttt                                                            <  4:818444/85‑1 (MQ=255)
                                                                                                                                             aaCGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTAttt                                                            >  3:818444/1‑85 (MQ=255)
                                                                                                                                                gCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttt               >  1:362333/1‑127 (MQ=255)
                                                                                                                                                gCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttt               <  2:362333/127‑1 (MQ=255)
                                                                                                                                                 ccTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt        >  4:215709/1‑133 (MQ=255)
                                                                                                                                                 ccTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt        <  3:215709/133‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                       
GTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCT  >  NC_000913/2917762‑2918045

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