Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,517,061 T→G 5.1% intergenic (+59/+38) 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,0610TG 5.1% 368.5 / 8.3 119intergenic (+59/+38)patD/yncLgamma‑aminobutyraldehyde dehydrogenase/stress‑induced small inner membrane enterobacterial protein
Reads supporting (aligned to +/- strand):  ref base T (57/55);  new base G (2/4);  total (60/59)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

Only 100 of 123 total aligned reads displayed.

GGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTGCCACCTTTTAT  >  NC_000913/1516913‑1517210
                                                                                                                                                    |                                                                                                                                                     
gggCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGAtt                                                                                                                                                      <  2:430045/150‑1 (MQ=255)
 ggCAGAAACTTTCTGTTTACGGCAAGGATATGTCACTTTATGCGGTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTTCAGATCGTAAGAGTACGAGATGATCTTGCGTACTGTCAGAGCACTGAttt                                                                                                                                                     <  2:1770862/150‑1 (MQ=255)
    aGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTtctc                                                                                                                                                  >  1:664676/1‑150 (MQ=255)
         cTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACa                                                                                                                                             >  1:1159775/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)
           ttCTGGGTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAGt                                                                                                                                           <  1:207725/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: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)
                              atGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCAGGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGAACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATc                                                                                                                        <  1:294263/150‑1 (MQ=255)
                               tGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCa                                                                                                                       <  1:476783/150‑1 (MQ=255)
                                gTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCAc                                                                                                                      >  1:450052/1‑150 (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)
                                          gggCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAAt                                                                                                            <  1:1273768/150‑1 (MQ=255)
                                           ggCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATa                                                                                                           <  1:451258/150‑1 (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)
                                                          ccGTCGTCCGCCACGGCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAAtctc                                                                                            <  1:1370558/150‑1 (MQ=255)
                                                            gtcgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCAt                                                                                          <  1:697769/150‑1 (MQ=255)
                                                              cgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATag                                                                                        <  2:1486279/150‑1 (MQ=255)
                                                              cgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATag                                                                                        <  1:1380772/150‑1 (MQ=255)
                                                              cgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATag                                                                                        <  1:1799524/150‑1 (MQ=255)
                                                               gtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATaga                                                                                       <  2:619118/150‑1 (MQ=255)
                                                                tcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATagag                                                                                      <  2:506694/150‑1 (MQ=255)
                                                                tcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATagag                                                                                      <  2:1741449/150‑1 (MQ=255)
                                                                  cGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGAt                                                                                    <  1:970452/150‑1 (MQ=255)
                                                                   accACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATa                                                                                   >  1:135682/2‑150 (MQ=255)
                                                                     cACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAag                                                                                 <  2:669515/150‑1 (MQ=255)
                                                                       cGTCATGGGTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAagag                                                                               <  2:1663973/150‑1 (MQ=255)
                                                                         tCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTa                                                                             <  1:110758/150‑1 (MQ=255)
                                                                             ggTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAAACAGCCAAATCTCATAGAGATAAGAGTaaaca                                                                         >  1:1688209/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)
                                                                                     aTTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAca                                                                 <  2:1841595/150‑1 (MQ=255)
                                                                                     aTTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGAATTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAca                                                                 >  1:213859/1‑150 (MQ=255)
                                                                                         aGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAcagcag                                                             <  2:885127/150‑1 (MQ=255)
                                                                                          ggATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTTCAGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACagcagc                                                            >  2:249133/1‑150 (MQ=255)
                                                                                          ggATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACagcagc                                                            >  2:179152/1‑150 (MQ=255)
                                                                                           gATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTTCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTTAAAAAACAAACAgcagca                                                           >  2:319727/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)
                                                                                                     cAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTa                                                 <  2:1538488/150‑1 (MQ=255)
                                                                                                       gATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATc                                               <  2:1159775/150‑1 (MQ=255)
                                                                                                       gATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATc                                               >  1:1332172/1‑150 (MQ=255)
                                                                                                           gTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGta                                           >  2:381077/1‑150 (MQ=255)
                                                                                                            tAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGtac                                          >  1:1420938/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)
                                                                                                                           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)
                                                                                                                               tGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAca                       >  1:563245/1‑150 (MQ=255)
                                                                                                                                  gTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACCTCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAACCCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattc                    >  1:99819/1‑150 (MQ=255)
                                                                                                                                  gTACTGTCAGAGAAAACAGTGCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattc                    >  1:1292832/1‑150 (MQ=255)
                                                                                                                                    aCTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattcat                  <  2:131707/150‑1 (MQ=255)
                                                                                                                                      tGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAATATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCAttt                >  2:764968/1‑150 (MQ=255)
                                                                                                                                         cAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTg             <  2:636636/150‑1 (MQ=255)
                                                                                                                                               gcTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCCTTTTTGCCACCt       >  1:1519121/2‑150 (MQ=255)
                                                                                                                                                 ttaTTTATCTGATAGTAAGTTCTGATTGAAAAATCACTACATCAAATAATCCAGCTAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGGTCTACTGGAAACATTCATTTTTGCCACCttt     <  1:289842/148‑1 (MQ=255)
                                                                                                                                                  gttttCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTTAGCAAACCCACAGCAGCAAAGAAATTTATCAGTACTAGAGTTCTCCTGGAAACATACATTTTTGACACCattt    >  2:1229133/3‑150 (MQ=255)
                                                                                                                                                   aTTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTGCCACCTTTTa   <  2:970657/150‑1 (MQ=255)
                                                                                                                                                    tttCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACACTTCTACTGGAAACATTCATTTTTGCCACAtattat  >  2:600500/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTGCCACCTTTTAT  >  NC_000913/1516913‑1517210

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: