Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,040 A→G 5.8% intergenic (‑427/+32) yhbX ← / ← leuU putative hydrolase, inner membrane/tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,322,0400AG 5.8% 371.2 / 9.1 120intergenic (‑427/+32)yhbX/leuUputative hydrolase, inner membrane/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (48/65);  new base G (4/3);  total (52/68)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 127 total aligned reads displayed.

AATAAAGTGCCAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGC  >  NC_000913/3321891‑3322189
                                                                                                                                                     |                                                                                                                                                     
aaTAAAGTGCCAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGGCCAGGACCAAAACGaaaaa                                                                                                                                                       >  2:919825/1‑150 (MQ=255)
 aTAAAGTGCCAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGaaaaaa                                                                                                                                                      >  1:1314567/1‑150 (MQ=255)
      gTGCCAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGCCCAGGACCAAAACGAAAAAAGACGc                                                                                                                                                 >  1:1427281/1‑150 (MQ=255)
      gTGCCAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGaaaaaa                                                                                                                                                      >  2:1236147/1‑145 (MQ=255)
      gTGCCAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGc                                                                                                                                                 >  1:1291272/1‑150 (MQ=255)
      gTGCCAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTACAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGAACAGGACCAAAACGAAAAAAGACga                                                                                                                                                 >  2:1406710/1‑149 (MQ=255)
       tGCCAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCt                                                                                                                                                <  1:902858/150‑1 (MQ=255)
        gCCAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCtt                                                                                                                                               <  2:1147653/150‑1 (MQ=255)
          cAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCtttt                                                                                                                                             >  1:144226/1‑150 (MQ=255)
           aaTTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTc                                                                                                                                            <  2:584242/150‑1 (MQ=255)
                 gTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtc                                                                                                                                      <  1:294366/150‑1 (MQ=255)
                  tCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTAAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtct                                                                                                                                     >  2:281020/1‑150 (MQ=255)
                     gTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCtt                                                                                                                                  <  2:1333063/150‑1 (MQ=255)
                        gAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAg                                                                                                                                                     >  1:220635/1‑128 (MQ=255)
                        gAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAg                                                                                                                                                     <  2:220635/128‑1 (MQ=255)
                         aaGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCt                                                                                                                              <  2:45089/150‑1 (MQ=255)
                             gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                          >  1:1064068/1‑150 (MQ=255)
                             gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                          >  1:419086/1‑150 (MQ=255)
                             gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                          >  1:798862/1‑150 (MQ=255)
                             gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                          >  1:1089186/1‑150 (MQ=255)
                             caaaaaTTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAGAGACGCTGAAAAGCGTCTTTTTTcgtttt                                                                                                                          <  1:135824/149‑6 (MQ=255)
                              aaaaaTTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAGAGACGCTGAAAAGCGTCTTTTTTcgttttg                                                                                                                         <  1:1117692/150‑7 (MQ=38)
                                   tttATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGt                                                                                                                    <  2:1499117/150‑1 (MQ=255)
                                   tttATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGt                                                                                                                    >  1:1343951/1‑150 (MQ=255)
                                     tataACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTAc                                                                                                                  >  2:1291878/1‑150 (MQ=255)
                                         aCTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAg                                                                                                              <  2:1372482/150‑1 (MQ=255)
                                         aCTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAg                                                                                                              <  1:83117/150‑1 (MQ=255)
                                          cTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAgg                                                                                                             >  1:710659/1‑150 (MQ=255)
                                          cTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAgg                                                                                                             <  2:1446369/150‑1 (MQ=255)
                                          cTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAgg                                                                                                             >  2:1010202/1‑150 (MQ=255)
                                            aaGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGAc                                                                                                           >  1:493854/1‑150 (MQ=255)
                                              ggCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACgg                                                                                                         <  2:556430/150‑1 (MQ=255)
                                               gCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACggg                                                                                                        <  2:798862/150‑1 (MQ=255)
                                                   aCCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACtt                                                                                                    >  1:93486/1‑150 (MQ=255)
                                                     cGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGa                                                                                                  <  2:224704/150‑1 (MQ=255)
                                                       gCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAAc                                                                                                >  2:933764/1‑150 (MQ=255)
                                                        cACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAAcc                                                                                               <  1:1459915/150‑1 (MQ=255)
                                                               gTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGc                                                                                        <  2:112713/150‑1 (MQ=255)
                                                                tttCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGcc                                                                                       <  1:1010202/150‑1 (MQ=255)
                                                                    aaGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTAt                                                                                   <  2:365342/150‑1 (MQ=255)
                                                                      ggTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTg                                                                                 >  1:586598/1‑150 (MQ=255)
                                                                      ggTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTg                                                                                 <  1:1143999/150‑1 (MQ=255)
                                                                        tACCAGCTACGAGTAAAGCAACGGGACGAGATCCCGATACCTGACACCCATTCCTCCGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGGCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTggg                                                                               <  2:1303569/150‑1 (MQ=255)
                                                                         aCCAGCTACGAGTAAAGCAACTTGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCAATTGGGc                                                                              <  2:1206977/150‑1 (MQ=255)
                                                                         aCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGc                                                                              <  2:767700/150‑1 (MQ=255)
                                                                         aCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGc                                                                              <  2:372105/150‑1 (MQ=255)
                                                                         aCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGc                                                                              <  2:1081311/150‑1 (MQ=255)
                                                                          ccAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGTGACTTGAACCCGTAAGCCCTATTGGGCa                                                                             <  1:62472/150‑1 (MQ=255)
                                                                          ccAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCa                                                                             >  2:718879/1‑150 (MQ=255)
                                                                           cAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTGTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTCAGCCCTATTGGGCAc                                                                            >  2:753339/1‑150 (MQ=255)
                                                                            aGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACt                                                                           <  2:419086/150‑1 (MQ=255)
                                                                               tACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTCTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTacc                                                                        <  2:1239053/150‑1 (MQ=255)
                                                                                aCGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTacca                                                                       >  1:1228521/1‑150 (MQ=255)
                                                                                 cGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTaccac                                                                      >  2:76901/1‑150 (MQ=255)
                                                                                  gAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGGACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTaccacc                                                                     >  1:1186505/1‑150 (MQ=255)
                                                                                   aGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCt                                                                    >  1:330399/1‑150 (MQ=255)
                                                                                     tAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCa                                                                  <  2:345151/150‑1 (MQ=255)
                                                                                      aaaGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCaa                                                                 >  2:1172102/1‑150 (MQ=255)
                                                                                        aGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAgg                                                               >  2:1024613/1‑150 (MQ=255)
                                                                                         gCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGt                                                              <  1:953285/150‑1 (MQ=255)
                                                                                             cTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCg                                                          <  1:1209407/150‑1 (MQ=255)
                                                                                               ggACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCgtg                                                        <  2:1308810/150‑1 (MQ=255)
                                                                                               ggACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCgtg                                                        <  2:1364743/150‑1 (MQ=255)
                                                                                               ggACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCgtg                                                        >  2:153535/1‑150 (MQ=255)
                                                                                                         cagatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAAtt                                              <  2:52983/150‑1 (MQ=255)
                                                                                                         cagatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAAtt                                              >  2:738327/1‑150 (MQ=255)
                                                                                                          agatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTGTTTTCTGGAATTTGGTACCGAGGACGGGAGTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTc                                             <  1:1406710/150‑1 (MQ=255)
                                                                                                          agatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTc                                             <  1:933764/150‑1 (MQ=255)
                                                                                                             tacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTccac                                          <  2:510268/150‑1 (MQ=255)
                                                                                                               ccTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTccacca                                        >  2:157925/1‑150 (MQ=255)
                                                                                                                 tGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCaccacc                                      <  2:722436/150‑1 (MQ=255)
                                                                                                                     aaCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCgg                                  <  1:718879/150‑1 (MQ=255)
                                                                                                                     aaCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCgg                                  <  1:1286237/150‑1 (MQ=255)
                                                                                                                      aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGc                                 >  1:671135/1‑150 (MQ=255)
                                                                                                                      aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGc                                 >  1:1288765/1‑150 (MQ=255)
                                                                                                                       ccATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCCCCTCAAGGTAGCGTGTCTACCAATTCCACCACCACGGCa                                >  1:1405568/1‑150 (MQ=255)
                                                                                                                        cATTCCTCAGACAAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCGCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCAc                               >  2:1111529/1‑150 (MQ=255)
                                                                                                                        aaTTCCTCAGACCAGTTCCAAAACTAAAAAAGACGGTTTTCAGAGTCTCTTTTCTGGAATTTGGTACCGAGGAGGGGAGTTGATCCCGTAAGGCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCAGCTCGGCag                               <  2:508389/149‑2 (MQ=255)
                                                                                                                          ttCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACgg                             <  2:713629/150‑1 (MQ=255)
                                                                                                                            ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                           <  2:12873/150‑1 (MQ=255)
                                                                                                                            ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                           >  1:348427/1‑150 (MQ=255)
                                                                                                                               cAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtact                        >  2:155102/1‑150 (MQ=255)
                                                                                                                                 gACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCCCCTCAAGGTAGCGTGTCTACCAAGTCCACCCCCTCGGCACGGAtactac                      >  1:154975/1‑150 (MQ=255)
                                                                                                                                 gACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtactac                      <  1:153535/150‑1 (MQ=255)
                                                                                                                                   ccAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACtt                    >  1:63533/1‑150 (MQ=255)
                                                                                                                                       gACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTa                <  2:493854/150‑1 (MQ=255)
                                                                                                                                         ccAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGt              >  1:1324361/1‑150 (MQ=255)
                                                                                                                                           ttccagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTc            >  1:1694/6‑150 (MQ=255)
                                                                                                                                           aaaaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTc            >  1:683382/1‑150 (MQ=255)
                                                                                                                                           aaaaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTc            >  1:64991/1‑150 (MQ=255)
                                                                                                                                           aaaaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAAACTACTTTTTAGTTc            >  1:600757/1‑150 (MQ=255)
                                                                                                                                             ccagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCgg          >  1:94015/4‑150 (MQ=255)
                                                                                                                                             aaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCgg          >  1:829397/1‑150 (MQ=255)
                                                                                                                                              cagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCggg         <  2:1694/148‑1 (MQ=255)
                                                                                                                                              aCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCggg         <  2:586598/150‑1 (MQ=255)
                                                                                                                                              aCGAAAAAAGACGCTTTGCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCggg         <  2:271436/150‑1 (MQ=255)
                                                                                                                                               cGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGa        <  1:157905/150‑1 (MQ=255)
                                                                                                                                                   aaaaGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATc    >  1:764478/1‑150 (MQ=255)
                                                                                                                                                   aaaaGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATc    >  2:92819/1‑150 (MQ=255)
                                                                                                                                                     aaGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGc  >  1:559725/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATAAAGTGCCAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGC  >  NC_000913/3321891‑3322189

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