Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,307,294 A→T 9.0% L169Q (CTG→CAG)  alsK ← D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,307,2940AT 9.0% 261.4 / 12.9 99L169Q (CTG→CAG) alsKD‑allose kinase
Reads supporting (aligned to +/- strand):  ref base A (41/48);  new base T (5/4);  total (46/53)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-01

CGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGA  >  NC_000913/4307145‑4307441
                                                                                                                                                     |                                                                                                                                                   
cGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCa                                                                                                                                                     >  2:1306564/1‑150 (MQ=255)
  aaaGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgg                                                                                                                                                   >  2:1325727/1‑150 (MQ=255)
  aaaGGGGCCTGTTTCGCATGGACGAAAAGACCGCGCAAAGGGTAATTTACGGGCTGGTGTTCGTACCAGCGTCTTAGCGCCATTCCATAGCAATTGGTTTCGACGCACCCACGATCGCCACACCCGCAGTGTTTGGTCATATCTCCCAgg                                                                                                                                                   <  2:746295/150‑1 (MQ=255)
   aaGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggg                                                                                                                                                  <  2:1652615/150‑1 (MQ=255)
    aGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTGTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                 <  1:67329/150‑1 (MQ=255)
     ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                <  2:1741294/150‑1 (MQ=255)
     ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                <  2:1547929/150‑1 (MQ=255)
     ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                <  1:836148/150‑1 (MQ=255)
     ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTgggag                                                                                                                                                <  2:1068448/150‑3 (MQ=255)
     ggggCGGTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTTCCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                <  2:846948/150‑1 (MQ=255)
      gggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                               >  2:379164/1‑150 (MQ=255)
      gggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                               <  2:1712795/150‑1 (MQ=255)
      gggCGGTTTCCGCATGGACGACAAGCTCGCGCAATGGGGAATTTCGGGGCTGTTGTTCGTACCAGCGTCGTAGCGCCATTCCAGAGCAATTGGTTGCCCGGCACCCAGGATTGCCACACGCGCAGTGTGGGGTCATATCTCCCAGGGGGa                                                                                                                                               <  1:1766436/150‑1 (MQ=255)
       ggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGGTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGat                                                                                                                                              <  1:1133333/150‑1 (MQ=255)
        gCGTTTTCCGCATGGACGAAAAGATCTCTCAATGGGTACTTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGACATTCCAGAGCAATTGGTTTCCACGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGata                                                                                                                                             <  2:243066/150‑1 (MQ=255)
          gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATg                                                                                                                                           >  2:387899/1‑150 (MQ=255)
            tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                         <  1:390105/150‑1 (MQ=255)
             ttCCGCATGGACGACAAGATCGCGCAATGGGGAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAcc                                                                                                                                        <  2:1147160/150‑1 (MQ=255)
             ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTGCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAcc                                                                                                                                        <  2:1065200/150‑1 (MQ=255)
              tCCGCATGGACGAAAAGCTCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAccc                                                                                                                                       <  2:462974/150‑1 (MQ=255)
              tCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAccc                                                                                                                                       <  1:1789187/150‑1 (MQ=255)
              tCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAccc                                                                                                                                       <  1:1306564/150‑1 (MQ=255)
                 gCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                                                                                                    >  1:687470/1‑150 (MQ=255)
                  cATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGtt                                                                                                                                   >  1:1479766/1‑150 (MQ=255)
                     ggACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGc                                                                                                                                >  2:460343/1‑150 (MQ=255)
                       aCGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCt                                                                                                                              >  2:7266/1‑150 (MQ=255)
                        cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcccg                                                                                                                             <  2:1210573/150‑3 (MQ=255)
                        cGAAAAGATCGCGCAATGCGTAATATCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTGGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                             <  1:952455/150‑1 (MQ=255)
                        cGAAAAGATCGCGCAAAGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCTCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                             <  2:114124/150‑1 (MQ=255)
                          aacaGATCCCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGcc                                                                                                                           <  1:1140170/147‑1 (MQ=255)
                          aaaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGcc                                                                                                                           <  2:1791172/150‑1 (MQ=255)
                          aaaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACACAGTTCGCCTGcc                                                                                                                           >  2:1486489/1‑150 (MQ=255)
                          aaaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGCTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCCGGCACCCAGGAGTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGcc                                                                                                                           <  1:1374175/150‑1 (MQ=255)
                                   cgcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtgt                                                                                                                  <  1:1078893/150‑1 (MQ=255)
                                    gcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGtgtg                                                                                                                 <  2:947995/150‑1 (MQ=255)
                                    gcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGGCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGtgtg                                                                                                                 <  2:1783027/150‑1 (MQ=255)
                                     cactGGGGAATTTCGGGGCTGTTGTTCGTACCAGCGGCTTAGCGCCATTCCAGAGCAATTGGTTGCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGc                                                                                                                <  2:920762/147‑1 (MQ=255)
                                     cAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGc                                                                                                                <  1:1300393/150‑1 (MQ=255)
                                       aTGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAACCAGTTCGCCTGCCCCACCGTGTGCAc                                                                                                              >  1:832703/1‑150 (MQ=255)
                                            tAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTAGGGTCAAATCTCCCACGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTc                                                                                                         >  2:1391668/1‑150 (MQ=255)
                                                  cGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCAACACGCGCAGTGTTGTGTCATATCTCCCAGGGGGATATGACCAAGTTCGCCTGCCACACCGTGTTCACACTTCCACGgc                                                                                                   >  2:186430/1‑150 (MQ=255)
                                                      gCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCAc                                                                                               >  2:1677249/1‑150 (MQ=255)
                                                         gttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGt                                                                                            >  1:1552315/1‑150 (MQ=255)
                                                           tgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACCCGCGCAGTGTTGGGTCATATCTCCCAGCGGGATATGAACCAGTTCGCCTGCCACCCAGTCTGCCCCCGTCCACGGCGCACCGTTc                                                                                          >  2:703392/1‑150 (MQ=255)
                                                           tgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                          <  1:699899/150‑1 (MQ=255)
                                                           tgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                          >  1:1442904/1‑150 (MQ=255)
                                                            gttCGTACCAGCCTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                         >  2:499659/1‑150 (MQ=255)
                                                              gcGTCCCAGCGTCGTAGCGCCATTCCAGAGCACTTGGTTCCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGGCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGGTCATc                                                                                       <  2:302679/149‑1 (MQ=255)
                                                                gTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCca                                                                                     >  1:1219405/1‑150 (MQ=255)
                                                                   ccAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAAGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCccac                                                                                  >  2:159464/1‑150 (MQ=255)
                                                                   ccAGCGTATTAGCGCCATTCCAGAGCAATTGGATTCCAGGCACCCAAGATTGCCACACGCGCAGTCTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCTCATGCCACACCGTGTGCACCCGTACACGGCGCCCCGTTAATacacac                                                                                  >  2:159466/1‑150 (MQ=255)
                                                                   ccAGCGGCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                  <  1:224986/150‑1 (MQ=255)
                                                                    ccgCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACt                                                                                 <  1:415131/148‑1 (MQ=255)
                                                                    cAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTTCACCCGTACACGGCGCACCGTTCATCCACACt                                                                                 <  2:533518/150‑1 (MQ=255)
                                                                      gCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACAGTTCATCCACACTGc                                                                               >  1:201610/1‑150 (MQ=255)
                                                                      gCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGGGCACCCGTCCACGGCGCACCGTTCATCCACACTGc                                                                               <  1:1034507/150‑1 (MQ=255)
                                                                        gTCTTAGCGCCATTCCAGAGCACTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCACACCGTTCATCCACACTGCGa                                                                             >  1:285561/1‑150 (MQ=255)
                                                                        gTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGa                                                                             >  1:1171243/1‑150 (MQ=255)
                                                                          cTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCACGGGGATATGACCCAGTTCGCCTGCCACCCCGTGTGCACCCGTCCACGGCGCACCGATCATCCACACTGCGacc                                                                           >  1:371493/1‑148 (MQ=255)
                                                                             aGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAcccc                                                                        >  2:828386/1‑150 (MQ=255)
                                                                             aGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAcccc                                                                        >  2:474174/1‑150 (MQ=255)
                                                                                 ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATTACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATcc                                                                    >  2:1760392/1‑150 (MQ=255)
                                                                                  cATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATccc                                                                   >  1:1720659/1‑150 (MQ=255)
                                                                                       cAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGTCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAc                                                              >  2:1315594/1‑150 (MQ=255)
                                                                                       cAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCGCACCGTTCCTCCACACTTCGATCCCCCTCCCCGTAc                                                              >  2:1461877/1‑150 (MQ=255)
                                                                                        agagCAATTGGGTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACAGTTTTTGTTAATGATAAGGGGTCAACGTTGATCTACACTGCGAACCCCATCCCCGTAcc                                                             <  2:1603262/150‑1 (MQ=255)
                                                                                           gCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgag                                                          >  1:429693/1‑150 (MQ=255)
                                                                                           gCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgag                                                          >  1:505093/1‑150 (MQ=255)
                                                                                             aaTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGACCCCCCTCCCCGTACCGAGAt                                                        >  2:895901/1‑150 (MQ=255)
                                                                                              aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGACCCCCATCCCCGTACCGAGATa                                                       >  1:155197/1‑150 (MQ=255)
                                                                                              aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGAAATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATa                                                       >  1:149623/1‑150 (MQ=255)
                                                                                                    ttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccgc                                                 >  2:1472902/1‑150 (MQ=255)
                                                                                                     tCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccgcc                                                <  2:687470/150‑1 (MQ=255)
                                                                                                       cAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAg                                              >  1:1116491/1‑150 (MQ=255)
                                                                                                             cccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAg                                        >  1:1825703/1‑150 (MQ=255)
                                                                                                                aGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACAACGTGTGCACCCGTCCACGGCGCACCGTTCATCCACCCTGCGACCCCCAACCCCGTACCCAGAGAGGCCGCCAGACCCAgttg                                     >  2:819666/1‑150 (MQ=255)
                                                                                                                   aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCTCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGttgttg                                  <  2:285561/150‑1 (MQ=255)
                                                                                                                   aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGttgttg                                  <  2:1171243/150‑1 (MQ=255)
                                                                                                                   aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGttgttg                                  <  2:257638/150‑1 (MQ=255)
                                                                                                                     tGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCg                                <  1:1325727/150‑1 (MQ=255)
                                                                                                                             gcgcAGTGTTGGGTCATAGCTCCCCGGGGGATATGACCCAGTGCGCCTGCCACACCGGGTGCACCCGTCCACGGCGCACCGGTCCGCCCCACTGCGAACCCCATCCCCGTACCAAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCgg                        <  1:873976/150‑1 (MQ=255)
                                                                                                                                cAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGttt                     <  1:1088082/150‑1 (MQ=255)
                                                                                                                                  gtgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTc                   <  1:186430/150‑1 (MQ=255)
                                                                                                                                   tgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTct                  <  1:828386/150‑1 (MQ=255)
                                                                                                                                   tgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTct                  <  1:474174/150‑1 (MQ=255)
                                                                                                                                   aactGGGTCATATCCCCCTGGGAGATATTACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCAACGTTCATCCACACTGCGAACCCCTTCCCCGGACCGAGATAGGCCGCCAGACCCAGTTGTTGCGTAAGGCGGTTTTct                  >  1:1612460/4‑150 (MQ=255)
                                                                                                                                    gtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTcta                 <  1:606449/150‑1 (MQ=255)
                                                                                                                                    gtTGGGGCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGGCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTcta                 <  1:524595/150‑1 (MQ=255)
                                                                                                                                     ttGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctac                <  2:1825703/150‑1 (MQ=255)
                                                                                                                                      tGGGTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctact               >  2:403547/1‑150 (MQ=255)
                                                                                                                                      tGGGTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctact               >  2:454305/1‑150 (MQ=255)
                                                                                                                                       gggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGACCCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctacta              >  1:1737328/1‑150 (MQ=255)
                                                                                                                                       gggTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGTGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctacta              >  2:1478022/1‑150 (MQ=255)
                                                                                                                                         gTCATATCCCCCTGGGAGATATGACCCAGTTCGACTGCCCCCCCGTGTGCACCCGTCCTCGGCGCACCGGTAAGACACACTGCGAGCCCCATCCCCGTACCGAGATAGACCGCCAGACCCAGTTGATGCGTAAGGAGGGTGTGTAATACg            >  1:1372934/1‑150 (MQ=255)
                                                                                                                                             tatCTCCCAGGGGGATATGCCCCAGTTCGCCTACCACACCGTGTGCACCCGTCCACGGCGCGCCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGCACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc        >  2:966077/1‑150 (MQ=255)
                                                                                                                                               tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAg      <  1:612494/150‑1 (MQ=255)
                                                                                                                                               tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACAATGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAg      <  2:699020/150‑1 (MQ=255)
                                                                                                                                                  cccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgag   >  2:1028827/1‑150 (MQ=255)
                                                                                                                                                   ccAGGTGGATATGACCCAGTTCGCCTGCCACACAGTGTGCTACCGTCCACTGCGCACAGTTCATCCACACTGCGAACACCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgaga  <  2:188106/150‑1 (MQ=255)
                                                                                                                                                   ccAGGGGGATATGACGCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgaga  <  2:408241/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGA  >  NC_000913/4307145‑4307441

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: ‑