Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,797,391 A→C 6.2% L185R (CTG→CGG)  pheT ← phenylalanine‑‑tRNA ligase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,797,3910AC 6.2% 308.2 / 9.9 97L185R (CTG→CGG) pheTphenylalanine‑‑tRNA ligase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (46/45);  new base C (3/3);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTC  >  NC_000913/1797242‑1797539
                                                                                                                                                     |                                                                                                                                                    
tcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCa                                                                                                                                                      <  2:1075625/150‑1 (MQ=255)
tcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCa                                                                                                                                                      >  1:29464/1‑150 (MQ=255)
   tCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCt                                                                                                                                                   >  2:229579/1‑150 (MQ=255)
    ccACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTg                                                                                                                                                  >  2:805019/1‑150 (MQ=255)
    ccACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTg                                                                                                                                                  >  1:781914/1‑150 (MQ=255)
     cacaGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTgg                                                                                                                                                 <  1:1142263/150‑1 (MQ=255)
      acaGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGt                                                                                                                                                >  1:546919/1‑150 (MQ=255)
        aGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTc                                                                                                                                              <  2:1357896/150‑1 (MQ=255)
        aGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTc                                                                                                                                              <  1:1086513/150‑1 (MQ=255)
        aGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTc                                                                                                                                              <  1:662192/150‑1 (MQ=255)
          cGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAg                                                                                                                                            <  1:834106/150‑1 (MQ=255)
           ggAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaac                                                                                                                                           <  1:1264822/150‑3 (MQ=255)
            gAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                                                                          <  2:895040/150‑4 (MQ=255)
               ttGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACgg                                                                                                                                       <  1:1061754/150‑1 (MQ=255)
                tGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGc                                                                                                                                      >  1:1249087/1‑150 (MQ=255)
                          aCGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCAcgc                                                                                                                            >  1:1388571/1‑150 (MQ=255)
                            gTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                                                                                          >  1:1333979/1‑150 (MQ=255)
                                aTGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAacc                                                                                                                      >  2:1550651/1‑149 (MQ=255)
                                 tGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAacac                                                                                                                     <  1:319980/150‑1 (MQ=255)
                                  gCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACAcc                                                                                                                    >  2:585443/1‑150 (MQ=255)
                                     ttttACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                                                                                 >  2:145795/1‑150 (MQ=255)
                                      tttACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatg                                                                                                                <  2:597763/150‑1 (MQ=255)
                                         aCCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                                                                             <  1:1621951/150‑1 (MQ=255)
                                          ccACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGc                                                                                                            >  1:493035/1‑150 (MQ=255)
                                           cacaCGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGcc                                                                                                           >  2:187944/1‑150 (MQ=255)
                                           cacaCGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGcc                                                                                                           >  1:1142777/1‑150 (MQ=255)
                                             cacGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTa                                                                                                         >  1:752654/1‑150 (MQ=255)
                                               cGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAg                                                                                                       <  2:1094469/150‑1 (MQ=255)
                                                ggCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGc                                                                                                      <  2:1249087/150‑1 (MQ=255)
                                                ggCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGc                                                                                                      <  1:1610848/150‑1 (MQ=255)
                                                ggCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGc                                                                                                      <  1:1478423/150‑1 (MQ=255)
                                                     aGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTc                                                                                                  >  1:697206/1‑149 (MQ=255)
                                                     aGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTc                                                                                                  <  2:697206/149‑1 (MQ=255)
                                                          aCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACg                                                                                            >  2:174351/1‑150 (MQ=255)
                                                              gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGttt                                                                                        >  1:572920/1‑150 (MQ=255)
                                                                gCAGGCTTCCCGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTgg                                                                                      <  1:1130395/150‑1 (MQ=255)
                                                                  aGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCg                                                                                    <  1:1183387/150‑1 (MQ=255)
                                                                   ggCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCa                                                                                                                         <  2:354067/112‑1 (MQ=255)
                                                                   ggCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCa                                                                                                                         >  1:354067/1‑112 (MQ=255)
                                                                   ggCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGt                                                                                   >  1:1575360/1‑150 (MQ=255)
                                                                      ttCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCAc                                                                                <  1:1194683/150‑1 (MQ=255)
                                                                       tCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGc                                                                                                      >  2:385636/1‑127 (MQ=255)
                                                                       tCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGc                                                                                                      <  1:385636/127‑1 (MQ=255)
                                                                         cGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCa                                                                                                                   >  2:1232823/1‑112 (MQ=255)
                                                                         cGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCa                                                                                                                   <  1:1232823/112‑1 (MQ=255)
                                                                          ggCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTc                                                                                  >  1:1023958/1‑144 (MQ=255)
                                                                          ggCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTc                                                                                  <  2:1023958/144‑1 (MQ=255)
                                                                              cTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAttt                                                                        >  1:117432/1‑150 (MQ=255)
                                                                               ttCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTc                                                                       <  1:805019/150‑1 (MQ=255)
                                                                               ttCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTc                                                                       >  1:365062/1‑150 (MQ=255)
                                                                               ttCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTc                                                                       <  1:499072/150‑1 (MQ=255)
                                                                                tCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCg                                                                      <  1:321771/150‑1 (MQ=255)
                                                                                  gACTGTAATCGGGAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGAt                                                                    <  2:1260328/150‑1 (MQ=255)
                                                                                   aCTGTAATCGGCAGCGTGTCGTCGATGGTCACACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATg                                                                   <  1:445860/150‑1 (MQ=255)
                                                                                    cTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATgg                                                                  <  2:725006/150‑1 (MQ=255)
                                                                                      gTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCAc                                                                                             >  1:167964/1‑121 (MQ=255)
                                                                                      gTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCAc                                                                                             <  2:167964/121‑1 (MQ=255)
                                                                                      gTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg                                                                >  1:1119403/1‑150 (MQ=255)
                                                                                       tAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtgt                                                               >  1:1376034/1‑150 (MQ=255)
                                                                                         aTCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGttt                                                                                        <  1:844891/123‑1 (MQ=255)
                                                                                         aTCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGttt                                                                                        >  2:844891/1‑123 (MQ=255)
                                                                                         aTCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttg                                                             <  1:1111100/150‑1 (MQ=255)
                                                                                         aTCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttg                                                             >  1:514428/1‑150 (MQ=255)
                                                                                                   tgtCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGtt                                                   >  1:395308/1‑150 (MQ=255)
                                                                                                    gtcgtcGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGttt                                                  >  1:594315/1‑150 (MQ=255)
                                                                                                     tcgtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTc                                                 <  2:1142777/150‑1 (MQ=255)
                                                                                                      cgtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAtt                                                                         <  1:1170559/125‑1 (MQ=255)
                                                                                                      cgtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAtt                                                                         >  2:1170559/1‑125 (MQ=255)
                                                                                                        tcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAgg                                              <  2:1575360/150‑1 (MQ=255)
                                                                                                         cgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGt                                             <  1:229579/150‑1 (MQ=255)
                                                                                                         cgATGGTCGCACAAACCGGAACGATTTAAGGTTGAACCAGCTGCAGCTGGTTCAGCACGGCAACGTCACGCGAAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGt                                             >  2:1143036/1‑150 (MQ=255)
                                                                                                              gctggttcAACCGGAAATCGTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCa                                        >  2:185828/8‑150 (MQ=255)
                                                                                                              gTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCa                                        <  2:605114/150‑1 (MQ=255)
                                                                                                                 gCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACgg                                     >  2:644979/1‑150 (MQ=255)
                                                                                                                   aCCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGat                                   >  1:982115/1‑150 (MQ=255)
                                                                                                                   aCCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGat                                   <  2:374061/150‑1 (MQ=255)
                                                                                                                    tcAACCGGAAATCGTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGata                                  >  1:501588/2‑150 (MQ=255)
                                                                                                                    ccAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGata                                  <  1:965960/150‑1 (MQ=255)
                                                                                                                     cAACCGGAAATCGTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCCCGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGatat                                 >  2:1631598/1‑150 (MQ=255)
                                                                                                                     cAACCGGAAATCGTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGatat                                 >  2:218771/1‑150 (MQ=255)
                                                                                                                     cAACCGGAAATCGTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGatat                                 >  1:1611477/1‑150 (MQ=255)
                                                                                                                        ccGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCgg                              <  2:1376034/150‑1 (MQ=255)
                                                                                                                        ccGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCgg                              <  1:585443/150‑1 (MQ=255)
                                                                                                                          ggAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTg                            <  1:1246936/150‑1 (MQ=255)
                                                                                                                           gAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGc                           <  2:29464/150‑1 (MQ=255)
                                                                                                                                aTTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATc                      >  1:669486/1‑150 (MQ=255)
                                                                                                                                    ccGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGgcg                  <  2:554585/150‑1 (MQ=255)
                                                                                                                                     cGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGgcgc                 >  1:1603179/1‑150 (MQ=255)
                                                                                                                                       gctGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCAt               >  2:470498/3‑150 (MQ=255)
                                                                                                                                        ttGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATc              <  1:187944/150‑1 (MQ=255)
                                                                                                                                        ctGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATc              >  1:869069/2‑150 (MQ=255)
                                                                                                                                          gAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCg            >  1:836791/1‑150 (MQ=255)
                                                                                                                                           aaCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgc           <  2:365062/150‑1 (MQ=255)
                                                                                                                                           aaCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgc           >  1:947488/1‑150 (MQ=255)
                                                                                                                                            aCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCATTCGGCGCATCCgca          <  2:546919/150‑1 (MQ=255)
                                                                                                                                            aCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATcccag          >  2:517700/1‑147 (MQ=255)
                                                                                                                                            aCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgca          <  1:324150/150‑1 (MQ=255)
                                                                                                                                             ccAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgcag         >  1:1496464/1‑150 (MQ=255)
                                                                                                                                                    cAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTc  >  2:682469/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTC  >  NC_000913/1797242‑1797539

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