Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,517,055 C→A 5.2% intergenic (+53/+44) 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,0550CA 5.2% 353.7 / 4.6 125intergenic (+53/+44)patD/yncLgamma‑aminobutyraldehyde dehydrogenase/stress‑induced small inner membrane enterobacterial protein
Reads supporting (aligned to +/- strand):  ref base C (62/54);  new base A (4/4);  total (66/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

GTGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTGC  >  NC_000913/1516911‑1517200
                                                                                                                                                |                                                                                                                                                 
gTGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGa                                                                                                                                              <  1:456054/150‑1 (MQ=255)
 tGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGAGCGTAAGAGTATGAGATGATATTGCGTACTGTCAGAGCACTGAt                                                                                                                                             >  2:378831/1‑150 (MQ=255)
    gCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACCCCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTtc                                                                                                                                          >  2:1609871/1‑150 (MQ=255)
    gCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTtc                                                                                                                                          >  2:110103/1‑150 (MQ=255)
      aGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTtctc                                                                                                                                        <  2:1143365/150‑1 (MQ=255)
       taaaCTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTctct                                                                                                                                       <  2:1775821/149‑1 (MQ=255)
       gAAACTTTCTGTTTACGGCATGTATATGTCACTTTATGGGCTGGAGGATTACACCTTCGTCCGCCACGTAATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTctct                                                                                                                                       <  2:1553208/150‑1 (MQ=255)
          aCTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCTCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGAc                                                                                                                                    <  1:986957/150‑1 (MQ=255)
          aCTTTCTGGTTACGGCAAGGAGATGTAACTTTATGGGCTGTAGGATTACACCGTCGTCCTCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGAc                                                                                                                                    <  2:1046249/150‑1 (MQ=255)
           cTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACa                                                                                                                                   >  1:461412/1‑150 (MQ=255)
           cTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACa                                                                                                                                   <  2:1013333/150‑1 (MQ=255)
            tttCTTGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAg                                                                                                                                  <  1:1334144/150‑1 (MQ=255)
            tttCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAg                                                                                                                                  <  2:916501/150‑1 (MQ=255)
              tCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTa                                                                                                                                <  1:1097514/150‑1 (MQ=255)
                tGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATATTTCGTACTGTCAGAGAACTGATTTCTCTTACAGTACg                                                                                                                              >  1:1194292/1‑150 (MQ=255)
                  gTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGtg                                                                                                                            >  1:1462075/1‑150 (MQ=255)
                   ttACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGtgc                                                                                                                           <  1:481980/150‑1 (MQ=255)
                    tACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGtgct                                                                                                                          >  2:1708633/1‑150 (MQ=255)
                     aCGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGtgctg                                                                                                                         >  2:148877/1‑150 (MQ=255)
                      cGGCAAGGATATGTCACATTATGGGCGGTAGGATAACACCGGCGTCCGCCACGTCATGGTTAAACATTAAGTATAATATTGCAGATCGTTAGAGTATGTGATGATCTTCCTTACTGGCAGAGAAATCATGGCTCTGACAGTACgctagat                                                                                                                        <  2:687355/150‑7 (MQ=255)
                           aGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCa                                                                                                                   <  1:586848/150‑1 (MQ=255)
                            ggATATGTCACTTTATGGTCTTGAGGATTACTACGTCGTCCGACACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGAGCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCaa                                                                                                                  <  2:1296435/150‑1 (MQ=255)
                                 tGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCa                                                                                                             >  1:2023802/1‑150 (MQ=255)
                                    cACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTa                                                                                                          >  2:1615456/1‑150 (MQ=255)
                                    cACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTa                                                                                                          >  1:97440/1‑150 (MQ=255)
                                    cACTTTATGGGCTGCAGGATTAGACCGTCGTCCGCCACGTCAGTTTTAAACATTAAGGATAATAATGCACATCGAAAGAGTAGGAGCGGATCTTGCGTACTGTCAGAGCACTGATTTCTCTCACAGTACGTGCTGCGTGCAAAATCACTa                                                                                                          <  2:318719/150‑1 (MQ=255)
                                      cTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGAATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACa                                                                                                        >  1:185606/1‑150 (MQ=255)
                                         tATTGGCTGGAGGATTACACCTGCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCa                                                                                                     <  1:1970102/150‑1 (MQ=255)
                                          aTGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCaa                                                                                                    >  2:1199354/1‑150 (MQ=255)
                                             ggCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTTAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAAtt                                                                                                 >  1:130801/1‑149 (MQ=255)
                                             ggCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCATATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATa                                                                                                 <  2:1612174/150‑1 (MQ=255)
                                               cTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGCTCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATaaa                                                                                               >  1:283677/1‑150 (MQ=255)
                                               cTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATaaa                                                                                               <  1:1591401/150‑1 (MQ=255)
                                                tGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAAc                                                                                              >  1:1097923/1‑150 (MQ=255)
                                                tGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGAGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAAc                                                                                              <  1:82501/150‑1 (MQ=255)
                                                   aggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCGGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAg                                                                                           >  1:613098/1‑150 (MQ=255)
                                                     gaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGcc                                                                                         <  1:1147944/150‑1 (MQ=255)
                                                     gaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGcc                                                                                         <  1:890461/150‑1 (MQ=255)
                                                       ttACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCaa                                                                                       >  1:277903/1‑150 (MQ=255)
                                                        tACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCaaa                                                                                      >  2:944644/1‑150 (MQ=255)
                                                        tACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCaaa                                                                                      >  2:2012491/1‑150 (MQ=255)
                                                           acCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAAtct                                                                                   >  1:1692123/1‑150 (MQ=255)
                                                            ccGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCATATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTATGACAGTACGTGCTGCTTGCAAAATCACTACATAAAATAAACAAGCCAAAtcta                                                                                  >  1:1308381/1‑149 (MQ=255)
                                                            ccGTCGTCCGCCACGTCATGGTTAAACATTAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCa                                                                                 >  1:1642956/1‑150 (MQ=255)
                                                             cgtcgtCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCa                                                                                 >  2:204257/1‑150 (MQ=255)
                                                                cgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATag                                                                              <  2:1583088/150‑1 (MQ=255)
                                                                actcCGCCAATTCATGGTTAACCATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCATGCGGACTGTAAGAGCAATGAGTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATag                                                                              <  1:1299201/148‑1 (MQ=255)
                                                                 gtcCGCCACGTCATGGTTAACCATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATaga                                                                             >  2:1345450/1‑150 (MQ=255)
                                                                 gtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATaga                                                                             >  1:1203095/1‑150 (MQ=255)
                                                                 gtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACAGATTTCTCTGACAGAACGTTCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATATAATaga                                                                             >  1:255165/1‑150 (MQ=255)
                                                                  tcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCGAATAAACCAGCCAAATCTCATagag                                                                            >  1:1478770/1‑150 (MQ=255)
                                                                   ccGCCACGTCATGGTTAAACATTAAGTATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGAGTACTGTCAGAGCACTGATTTCTCTGAAAGTACGTGCTGCGTGCAAAATAACTACATCAAATAAACCAGCCAAATCTCATAgaga                                                                           >  1:1672991/1‑150 (MQ=255)
                                                                    cGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGAt                                                                          <  2:1666277/150‑1 (MQ=255)
                                                                    cGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGAt                                                                          >  2:620489/1‑150 (MQ=255)
                                                                     gCCACGTCATGGGTAAACATTAAGGATAATATTGCCGATCGTAAGAGTAGGAGATGATCTTGCGTACTGTCAGAGCACTTAGTTCTCTGACAGTACGTGCTGCGTGCAAAATCACAACATCAAATCAACCAGCCAAATCTCATAGAGATa                                                                         <  2:1896951/150‑1 (MQ=255)
                                                                        aCGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAaga                                                                      <  2:700914/150‑1 (MQ=255)
                                                                          gTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGt                                                                    >  1:1986971/1‑150 (MQ=255)
                                                                          gTCAGGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGt                                                                    <  1:625027/150‑1 (MQ=255)
                                                                              tGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaac                                                                >  2:1668235/1‑150 (MQ=255)
                                                                              tGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaac                                                                <  1:669106/150‑1 (MQ=255)
                                                                                 ttAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaacaaa                                                             <  2:771395/150‑1 (MQ=255)
                                                                                 ttAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaacaaa                                                             >  1:1691356/1‑150 (MQ=255)
                                                                                   aaaCATTAAGTATAATATTGCAGATCGTAAGAGTATGAGATGTTCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTAATAGAGATAAGAGTATACAAAcc                                                           <  2:465132/150‑1 (MQ=255)
                                                                                    aaCATTAAGGATAATATTGCAGATCTTATGATTATGAGATGATATTTCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCCCTGCATCAAATAAACACGCCAAATCTCATAGAGAAAAGAGTAAACAAACCa                                                          >  2:1239870/1‑150 (MQ=255)
                                                                                     aCATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCaa                                                         >  1:694523/1‑150 (MQ=255)
                                                                                      cATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCCGCCAAATCTCATAGAGATAAGAGTAACCAAACCAAc                                                        >  2:1884788/1‑150 (MQ=255)
                                                                                      cATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAc                                                        >  1:836747/1‑150 (MQ=255)
                                                                                      cATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATAAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAc                                                        >  2:943470/1‑150 (MQ=255)
                                                                                      aaTTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAc                                                        <  2:917372/149‑1 (MQ=255)
                                                                                        ttAAGTATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGAACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACAGCAAATACACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAcag                                                      <  1:112062/150‑1 (MQ=255)
                                                                                             gATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAgcagca                                                 <  1:959906/150‑1 (MQ=255)
                                                                                             gATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAgcagca                                                 <  1:583621/150‑1 (MQ=255)
                                                                                                   atTGCAGATCGTAAGAGTAGGAGATGATCTTGCGTACTGTCAGAGCACTGATTGCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGaaa                                           <  2:1443863/150‑1 (MQ=255)
                                                                                                     tGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAAtt                                         <  1:1444374/150‑1 (MQ=255)
                                                                                                     tGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAAtt                                         >  2:1707183/1‑150 (MQ=255)
                                                                                                      gCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAAttt                                        >  1:1844410/1‑150 (MQ=255)
                                                                                                      gCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAAttt                                        >  2:2055244/1‑150 (MQ=255)
                                                                                                      gCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAAttt                                        >  2:269125/1‑150 (MQ=255)
                                                                                                       cAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTa                                       >  2:362306/1‑150 (MQ=255)
                                                                                                        aGATCGTAAGAGTATTAGATGATCTTGCGTACTGTCAGATCAATGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTAt                                      <  2:74797/150‑1 (MQ=255)
                                                                                                         gATCGTAAGAGTATGAGATGAGCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATc                                     <  2:340000/150‑1 (MQ=255)
                                                                                                          aTCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCa                                    <  1:110103/150‑1 (MQ=255)
                                                                                                          aTCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCa                                    >  1:587697/1‑150 (MQ=255)
                                                                                                            cGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGt                                  >  1:1951155/1‑150 (MQ=255)
                                                                                                             gTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGta                                 <  1:1616585/150‑1 (MQ=255)
                                                                                                               aaGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGtact                               <  1:717971/150‑1 (MQ=255)
                                                                                                                agagTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGtacta                              >  1:1204296/1‑150 (MQ=255)
                                                                                                                agagTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGtacta                              >  1:2043546/1‑150 (MQ=255)
                                                                                                                  agTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTCCGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATTAGAGTAAACAAACCGACAGCCGCTAAGAAATTTATCAGTACTACa                            >  1:70962/1‑150 (MQ=255)
                                                                                                                  agTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACa                            >  1:780391/1‑150 (MQ=255)
                                                                                                                   gTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAg                           >  1:1721564/1‑150 (MQ=255)
                                                                                                                        agaTGATCTTGCGTACTGTCAGAGCACTGATATCTCTGACAGTCCGTGCTGCGTGCAAAATCACTTCATCAAATAATCAAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACAACAGTTCTa                      >  1:289570/1‑150 (MQ=255)
                                                                                                                         tatgatCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTAc                     <  1:438763/149‑1 (MQ=255)
                                                                                                                           tgatCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCTAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTg                   >  2:763879/1‑150 (MQ=255)
                                                                                                                            gatCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTgg                  >  1:1430266/1‑150 (MQ=255)
                                                                                                                                       cTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATATATCAGTACTACAGTTCTACTGGAAAcattcatt       >  1:1619955/1‑150 (MQ=255)
                                                                                                                                         gTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCAtttt     <  2:2023802/150‑1 (MQ=255)
                                                                                                                                          tCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCAttttt    <  2:1053408/150‑1 (MQ=255)
                                                                                                                                           aattgCACTGAATTCTATGACAGTTCGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACATCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAATCATTCATTTTTg   <  1:1814165/146‑1 (MQ=255)
                                                                                                                                            agagCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCACCCAAATCTCATAGAGATAAGAGTAAAAAAACCAACACCAGCAAAGAAATTTATCAGAACTACAGTTCTAATGGAAACATTCATTATTGc  >  2:690177/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTGC  >  NC_000913/1516911‑1517200

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: