Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,517,059 G→C 5.0% intergenic (+57/+40) patD → / ← yncL gamma‑aminobutyraldehyde dehydrogenase/stress‑induced small inner membrane enterobacterial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,517,0590GC 5.0% 385.4 / 9.0 119intergenic (+57/+40)patD/yncLgamma‑aminobutyraldehyde dehydrogenase/stress‑induced small inner membrane enterobacterial protein
Reads supporting (aligned to +/- strand):  ref base G (57/55);  new base C (2/4);  total (59/60)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 121 total aligned reads displayed.

GGTGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTGCCACCTTTT  >  NC_000913/1516910‑1517208
                                                                                                                                                     |                                                                                                                                                     
ggTGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTg                                                                                                                                                       <  2:1188747/150‑1 (MQ=255)
       aGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTtctc                                                                                                                                                >  1:664676/1‑150 (MQ=255)
            cTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACa                                                                                                                                           >  1:1159775/1‑150 (MQ=255)
            cTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGAATGAGATGATCTTGCTTACTGTCAGAGCACTGATTTCTCTGCCa                                                                                                                                           >  1:1356505/1‑150 (MQ=255)
              ttCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGt                                                                                                                                         <  1:1031063/150‑1 (MQ=255)
              ttCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAGt                                                                                                                                         <  1:426069/150‑1 (MQ=255)
              ttCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAGt                                                                                                                                         <  1:1861232/150‑1 (MQ=255)
               tCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTa                                                                                                                                        <  1:1015063/150‑1 (MQ=255)
               tCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTa                                                                                                                                        <  1:866015/150‑1 (MQ=255)
                cTGGTTTCGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTAc                                                                                                                                       <  2:1126544/150‑1 (MQ=255)
                  ggTTACGGCAAGGATATGTCACTGTATGGGCTGGAGGATTACACCGTAGTCCGACACGTAATTGTTAAACATTAAGGATAATATTGCAGATCGGAAGAGTAGGAGAGGAGCTTGCGTACTGTCAGAGCACTGATTTCTATGAGAGTACGt                                                                                                                                     >  2:309841/1‑150 (MQ=255)
                   gTTACGGCAAGTATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAGTACgca                                                                                                                                    <  1:844426/150‑3 (MQ=39)
                    ttccGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACTCCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGCGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGtgc                                                                                                                                   <  1:1309580/147‑1 (MQ=255)
                    ttACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGtgc                                                                                                                                   <  2:914864/150‑1 (MQ=255)
                     tACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCCTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGtgct                                                                                                                                  <  1:1061555/150‑1 (MQ=255)
                       cGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGtgctgc                                                                                                                                >  2:714392/1‑150 (MQ=255)
                             ggATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCaa                                                                                                                          >  1:1263867/1‑150 (MQ=255)
                             ggATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCaa                                                                                                                          >  1:558250/1‑150 (MQ=255)
                             ggATATGTCACTTTATGGGCTGGAGGATTACACCGTCGGCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCaa                                                                                                                          <  1:1025867/150‑1 (MQ=255)
                              gATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCaaa                                                                                                                         <  1:139286/150‑1 (MQ=255)
                                tatGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAAt                                                                                                                       >  1:636636/1‑150 (MQ=255)
                                   gTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCAc                                                                                                                    >  1:450052/1‑150 (MQ=255)
                                   gTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCAc                                                                                                                    <  1:245199/150‑1 (MQ=255)
                                   gTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCATGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCAc                                                                                                                    <  1:538174/150‑1 (MQ=255)
                                     cACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTa                                                                                                                  >  2:649827/1‑150 (MQ=255)
                                      aCTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGCTTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTAc                                                                                                                 >  1:260978/1‑150 (MQ=255)
                                       cTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACa                                                                                                                <  1:462489/150‑1 (MQ=255)
                                       atttATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACa                                                                                                                <  1:1705767/149‑1 (MQ=255)
                                         ttATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATc                                                                                                              <  1:1117917/150‑1 (MQ=255)
                                          tATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCa                                                                                                             <  2:1792777/150‑1 (MQ=255)
                                            ttggCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCATAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCaaa                                                                                                           <  1:81156/148‑1 (MQ=255)
                                              ggCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATa                                                                                                         <  1:451258/150‑1 (MQ=255)
                                                cTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATaaa                                                                                                       >  1:1847444/1‑150 (MQ=255)
                                                 tGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCTAAATCACTACATCAAATAAAc                                                                                                      >  2:1239846/1‑150 (MQ=255)
                                                  ggaggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAAcc                                                                                                     >  1:1830826/1‑150 (MQ=255)
                                                   gaggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTATTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCa                                                                                                    <  1:854152/150‑1 (MQ=255)
                                                   gaggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCATAATCACTACATCAAATAAAcgc                                                                                                    >  1:1334118/1‑148 (MQ=255)
                                                   gaggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAATCCa                                                                                                    >  2:978206/1‑150 (MQ=255)
                                                   gaggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCa                                                                                                    >  1:219186/1‑150 (MQ=255)
                                                    aggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCCGCGTGCAAAATCACTACATCAAATAAACCAg                                                                                                   >  2:1399988/1‑150 (MQ=255)
                                                     ggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTTAGGATAATATTGCAGATCGTAAGAGTAGGAGATGATCTTGCGTACTGTCAGAGCACGGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATACACCAGc                                                                                                  >  2:1364424/1‑150 (MQ=255)
                                                     ggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGc                                                                                                  >  2:53049/1‑150 (MQ=255)
                                                      gaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGcc                                                                                                 <  1:526522/150‑1 (MQ=255)
                                                      gaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGcc                                                                                                 >  1:1164131/1‑150 (MQ=255)
                                                        ttACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCaa                                                                                               <  1:1367676/150‑1 (MQ=255)
                                                            acCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAAtct                                                                                           <  2:834559/150‑1 (MQ=255)
                                                             ccGTCGTCCGCCACGGCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAAtctc                                                                                          <  1:1370558/150‑1 (MQ=255)
                                                               gtcgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCAt                                                                                        <  1:697769/150‑1 (MQ=255)
                                                                 cgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATag                                                                                      <  1:1380772/150‑1 (MQ=255)
                                                                 cgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATag                                                                                      <  1:1799524/150‑1 (MQ=255)
                                                                 cgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATag                                                                                      <  2:1486279/150‑1 (MQ=255)
                                                                   tcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATagag                                                                                    <  2:506694/150‑1 (MQ=255)
                                                                   tcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATagag                                                                                    <  2:1741449/150‑1 (MQ=255)
                                                                    ccGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCTTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAgaga                                                                                   <  1:164140/150‑1 (MQ=255)
                                                                     cGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGAt                                                                                  <  1:970452/150‑1 (MQ=255)
                                                                        cACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAag                                                                               <  2:669515/150‑1 (MQ=255)
                                                                          cGTCATGGGTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAagag                                                                             <  2:1663973/150‑1 (MQ=255)
                                                                            tCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTa                                                                           <  1:110758/150‑1 (MQ=255)
                                                                              aTGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaa                                                                         >  1:606437/1‑150 (MQ=255)
                                                                                ggTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaaca                                                                       >  2:1078134/1‑150 (MQ=255)
                                                                                  ttAAACATTAAGGATAATATTGCAGATAGTAAGAGTATGAGATGATCTTGCGTACTTTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaacaaa                                                                     <  1:918046/150‑1 (MQ=255)
                                                                                   tAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGAGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaacaaac                                                                    <  1:269985/150‑1 (MQ=255)
                                                                                      aCATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCATATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCaa                                                                 >  1:1020679/1‑150 (MQ=255)
                                                                                       cATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAc                                                                >  1:1662372/1‑150 (MQ=255)
                                                                                        aTTAAGGATAATATTGCAGCTCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTTCAGCGTGCAAAATCACTTACTAAAAAAAACCAGCCAAATCTCATAGAGATAAGAGAAaacaaacaaaca                                                               >  2:512493/1‑150 (MQ=255)
                                                                                        aTTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAca                                                               <  2:1841595/150‑1 (MQ=255)
                                                                                        aTTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGAATTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAca                                                               >  1:213859/1‑150 (MQ=255)
                                                                                            aGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAcagcag                                                           <  2:885127/150‑1 (MQ=255)
                                                                                             ggATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACagcagc                                                          >  2:179152/1‑150 (MQ=255)
                                                                                              gATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTTCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTTAAAAAACAAACAgcagca                                                         >  2:319727/1‑150 (MQ=255)
                                                                                                   tatTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGaa                                                    >  2:344609/1‑150 (MQ=255)
                                                                                                   tatTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACAAACAGCAGCAAAGaa                                                    >  1:1383628/1‑150 (MQ=255)
                                                                                                       gCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAAttt                                                >  1:970657/1‑150 (MQ=255)
                                                                                                       gCAGATCGTAAGAGTATGAGAAGATCTTGCGTACTGTAAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCTAAATAACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAAATCAGCAAAGAAAttt                                                >  2:126756/1‑150 (MQ=255)
                                                                                                        cAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTa                                               >  1:1438182/1‑150 (MQ=255)
                                                                                                          gATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATc                                             >  2:1511933/1‑150 (MQ=255)
                                                                                                          gATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATc                                             >  1:1332172/1‑150 (MQ=255)
                                                                                                          gATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATc                                             <  2:1159775/150‑1 (MQ=255)
                                                                                                            tCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAg                                           >  2:643077/1‑150 (MQ=255)
                                                                                                              gTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGta                                         >  2:381077/1‑150 (MQ=255)
                                                                                                                cagagTATGAGCTGATCTTGCGTACTGTCAGAGCTCTTATTTCTCTGACAGTACGTGCGGCGTGCAAAATCACTACATGTAATAAACCAGCCAACTCTCAAAGAGATAAGAGTAAACAATCCAACAGCATCAAAGAAATTTAACAGtact                                       <  2:299587/149‑1 (MQ=255)
                                                                                                                        gagaTGCTCTTGCGTACTGTCAGAGCCCTGATTTCTCTGCCAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATATAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCt                               >  2:1136168/1‑150 (MQ=255)
                                                                                                                         agaTGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTa                              >  2:1765989/1‑150 (MQ=255)
                                                                                                                         agaTGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACGACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTa                              >  2:1176009/1‑150 (MQ=255)
                                                                                                                           attatCTTGCGTACTGTAAGAGCACTGATTTATCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACt                            <  2:1342346/147‑1 (MQ=255)
                                                                                                                           atgatCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACt                            <  2:5785/150‑1 (MQ=255)
                                                                                                                            tgatCTTGCGTACTGTCAGAGCACTGATTTGTCTGACAGTAGGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCATAGAAATTTATCAGTACTACAGTTCTACTg                           >  2:511431/1‑150 (MQ=255)
                                                                                                                             gatCTTGCGTACTGTCAGAGCACTGATTTATCTGACAGTAAGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCATATCTCATAGAGATAAGAGTAAACAAACCAAAAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACttg                          >  1:478252/1‑148 (MQ=255)
                                                                                                                              atCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGa                         >  2:1751129/1‑150 (MQ=255)
                                                                                                                              atCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGa                         <  2:1438182/150‑1 (MQ=255)
                                                                                                                                cTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGaaa                       >  2:1572770/1‑150 (MQ=255)
                                                                                                                                cTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGaaa                       >  1:840470/1‑150 (MQ=255)
                                                                                                                                  tGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAca                     <  1:649827/150‑1 (MQ=255)
                                                                                                                                     gTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACCTCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAACCCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattc                  >  1:99819/1‑150 (MQ=255)
                                                                                                                                     gTACTGTCAGAGAAAACAGTGCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattc                  >  1:1292832/1‑150 (MQ=255)
                                                                                                                                       aCTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattcat                <  2:131707/150‑1 (MQ=255)
                                                                                                                                        cTGTCGGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAACCCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattcatt               >  2:1870543/1‑150 (MQ=255)
                                                                                                                                            cAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTg           <  2:636636/150‑1 (MQ=255)
                                                                                                                                            aagagCACTGATTTCTCTGACAGTAAGTGCTGCTTGCAAAATCACTACATCAAATAAACAAGCCAAATCTCATAGAGATATGAGTAAACAAACCAACAGCAGCTAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTg           <  1:309841/149‑1 (MQ=255)
                                                                                                                                                  gcTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCCTTTTTGCCACCt     >  1:1519121/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTGCCACCTTTT  >  NC_000913/1516910‑1517208

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