Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,075,232 G→A 5.2% A152T (GCA→ACA)  yeeR → CP4‑44 prophage; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,075,2320GA 5.2% 396.3 / 7.0 114A152T (GCA→ACA) yeeRCP4‑44 prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (43/65);  new base A (3/3);  total (46/68)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

CCGCTCTCGCCGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGACCTGTGCCATCAGCCAGACA  >  NC_000913/2075083‑2075381
                                                                                                                                                     |                                                                                                                                                     
 cGCTCTCGCCGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAGCGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATgcgc                                                                                                                                                      >  2:42751/1‑150 (MQ=255)
 cGCTCTCGCCGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATgcgc                                                                                                                                                      >  1:799221/1‑150 (MQ=255)
  gCTCTCGCCGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCa                                                                                                                                                     >  2:828288/1‑150 (MQ=255)
  gCTCTCGCCGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCa                                                                                                                                                     >  2:1668098/1‑150 (MQ=255)
      tcGCCGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCAt                                                                                                                                                 <  2:1013188/150‑1 (MQ=255)
      tcGCCGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCAt                                                                                                                                                 <  2:1013209/150‑1 (MQ=255)
         ccGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGa                                                                                                                                              >  2:708110/1‑150 (MQ=255)
          cGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGAt                                                                                                                                             <  1:1677650/150‑1 (MQ=255)
          cGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGAt                                                                                                                                             <  2:117920/150‑1 (MQ=255)
             tggtCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTGTCCTTATGGCGTGTTCTGCATGCGCACCATTGATcgc                                                                                                                                          <  1:1541195/150‑1 (MQ=255)
               gtCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGt                                                                                                                                        <  1:657444/150‑1 (MQ=255)
                tCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGtt                                                                                                                                       >  2:1340431/1‑150 (MQ=255)
                   gctggctgTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACa                                                                                                                                    <  1:343612/150‑1 (MQ=255)
                   gctggctgTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACa                                                                                                                                    <  1:343583/150‑1 (MQ=255)
                   gctggctgTTCACCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACa                                                                                                                                    <  1:1317824/150‑1 (MQ=255)
                       gctgTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTAc                                                                                                                                >  1:1377664/1‑150 (MQ=255)
                       gctgTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTAc                                                                                                                                <  2:1464435/150‑1 (MQ=255)
                              cccTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATAttt                                                                                                                         <  1:798956/150‑1 (MQ=255)
                              cccTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATAttt                                                                                                                         >  1:638612/1‑150 (MQ=255)
                              cccTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATAttt                                                                                                                         >  1:789325/1‑150 (MQ=255)
                              cccTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATAttt                                                                                                                         >  1:789346/1‑150 (MQ=255)
                                  gCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCAcc                                                                                                                     <  1:480683/150‑1 (MQ=255)
                                   cAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCg                                                                                                                    >  1:736314/1‑150 (MQ=255)
                                   cAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCg                                                                                                                    <  2:1630763/150‑1 (MQ=255)
                                   cAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCg                                                                                                                    <  2:1630764/150‑1 (MQ=255)
                                      aTCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCa                                                                                                                 >  1:859059/1‑150 (MQ=255)
                                      aTCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCa                                                                                                                 >  1:1297365/1‑150 (MQ=255)
                                      aTCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCa                                                                                                                 >  1:479654/1‑150 (MQ=255)
                                         aaaaaCATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTg                                                                                                              >  2:1147662/1‑150 (MQ=255)
                                            aaCATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGAtg                                                                                                           >  2:737799/1‑150 (MQ=255)
                                              cATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGAtgtg                                                                                                         <  2:1308263/150‑1 (MQ=255)
                                              cATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGAtgtg                                                                                                         <  2:1156289/150‑1 (MQ=255)
                                              cATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGAtgtg                                                                                                         <  2:342648/150‑1 (MQ=255)
                                              cATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTATTTATACTGCCTATTTCCTTATGGCGTGTTCTGCATGCGCACCCTTGATCGCGTTACAGTACGATATTTCACCGTCACTGAtgtg                                                                                                         <  2:342651/150‑1 (MQ=255)
                                                ttCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATgtggt                                                                                                       >  1:377625/1‑150 (MQ=255)
                                                ttCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATgtggt                                                                                                       >  2:1559138/1‑150 (MQ=255)
                                                ttCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATgtggt                                                                                                       >  2:1559133/1‑150 (MQ=255)
                                                 tCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTTACTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCTCCGTCACTGATgtggtg                                                                                                      <  1:1170800/150‑1 (MQ=255)
                                                    gcagAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTgggg                                                                                                   >  2:979720/1‑150 (MQ=255)
                                                     cggAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCACACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGc                                                                                                  >  2:87363/3‑150 (MQ=255)
                                                     cagAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTGTGCATGCACACCATTGATCGCGTTACAGTACGATATTTCACCGTCACGGAGGTGGTGGGGc                                                                                                  >  2:1339332/1‑150 (MQ=255)
                                                     cagAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCACACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGc                                                                                                  >  2:1339309/1‑150 (MQ=255)
                                                      agAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCg                                                                                                 <  1:538829/150‑1 (MQ=255)
                                                      agAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCg                                                                                                 >  1:926075/1‑150 (MQ=255)
                                                      agAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCg                                                                                                 <  1:1564985/150‑1 (MQ=255)
                                                       gAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCgg                                                                                                <  1:1546923/150‑1 (MQ=255)
                                                       gAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCgg                                                                                                <  2:1176565/150‑1 (MQ=255)
                                                       gAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCgg                                                                                                <  2:1604760/150‑1 (MQ=255)
                                                        aaCGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCggg                                                                                               <  1:166065/150‑1 (MQ=255)
                                                        aaCGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCggg                                                                                               >  1:1182608/1‑150 (MQ=255)
                                                             gATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGgttgtt                                                                                          <  2:1640794/150‑1 (MQ=255)
                                                             gATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGgttgtt                                                                                          <  2:473897/150‑1 (MQ=255)
                                                             gATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGgttgtt                                                                                          <  2:249627/150‑1 (MQ=255)
                                                                 acacCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTAc                                                                                      <  2:520544/150‑1 (MQ=255)
                                                                 acacCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTAc                                                                                      >  1:1578651/1‑150 (MQ=255)
                                                                     cAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTActggc                                                                                  >  2:983329/1‑150 (MQ=255)
                                                                     cAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTActggc                                                                                  <  1:1692390/150‑1 (MQ=255)
                                                                      agaagaTATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTActggct                                                                                 <  1:168777/150‑1 (MQ=255)
                                                                      agaagaTATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTActggct                                                                                 <  1:1041343/150‑1 (MQ=255)
                                                                          gaTATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTggctggct                                                                             <  2:352145/150‑1 (MQ=255)
                                                                          gaTATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTggctggct                                                                             <  2:352118/150‑1 (MQ=255)
                                                                          gaTATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTggctggct                                                                             <  2:788973/150‑1 (MQ=255)
                                                                              tCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCAt                                                                         >  1:1063814/1‑150 (MQ=255)
                                                                                aGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTa                                                                       >  1:779823/1‑150 (MQ=255)
                                                                                aGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTa                                                                       >  1:1115166/1‑150 (MQ=255)
                                                                                   aaGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTg                                                                    >  2:702776/1‑150 (MQ=255)
                                                                                           tAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACgctgct                                                            <  1:708110/150‑1 (MQ=255)
                                                                                           tAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACgctgct                                                            <  1:264820/150‑1 (MQ=255)
                                                                                                ggTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATggg                                                       <  2:211286/150‑1 (MQ=255)
                                                                                                ggTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATggg                                                       >  2:292961/1‑150 (MQ=255)
                                                                                                 gTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGggc                                                      >  2:85656/1‑150 (MQ=255)
                                                                                                   aTTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGggcgg                                                    <  1:617188/150‑1 (MQ=255)
                                                                                                      aCTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGgcca                                                 <  1:1115344/150‑1 (MQ=255)
                                                                                                      aCTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGgcca                                                 <  2:1377664/150‑1 (MQ=255)
                                                                                                      aCTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGgcca                                                 <  2:978291/150‑1 (MQ=255)
                                                                                                      aCTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGgcca                                                 <  2:104252/150‑1 (MQ=255)
                                                                                                        tgctgTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGgccagc                                               >  2:919968/1‑150 (MQ=255)
                                                                                                                      tGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCt                                 <  1:159319/150‑1 (MQ=255)
                                                                                                                      tGCCTTTGTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCt                                 <  1:1031787/150‑1 (MQ=255)
                                                                                                                       gCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCAGTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGGATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTg                                >  2:736739/1‑150 (MQ=255)
                                                                                                                       gCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTg                                >  2:666575/1‑150 (MQ=255)
                                                                                                                            tttCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGaaaaa                           >  1:326084/1‑150 (MQ=255)
                                                                                                                             ttCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGaaaaaa                          <  2:41847/150‑1 (MQ=255)
                                                                                                                             ttCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGaaaaaa                          <  1:404503/150‑1 (MQ=255)
                                                                                                                               ccTTATGGCGTGTTCTGCATGCACACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACt                        <  1:87363/150‑1 (MQ=255)
                                                                                                                               ccTTATGGCGTGTTCTGCATGCACACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACt                        <  1:1339309/150‑1 (MQ=255)
                                                                                                                               ccTTATGGCGTGTTCTGCATGCACACCATTGATCGCGGTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACt                        <  1:1339332/150‑1 (MQ=255)
                                                                                                                                 ttATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGa                      <  2:1195565/150‑1 (MQ=255)
                                                                                                                                 ttATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGa                      <  2:1555897/150‑1 (MQ=255)
                                                                                                                                 ttATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGa                      >  2:674418/1‑150 (MQ=255)
                                                                                                                                  tATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGAc                     <  2:638612/150‑1 (MQ=255)
                                                                                                                                      gCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGACCtgt                 <  2:120299/150‑1 (MQ=255)
                                                                                                                                       cGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCCGGCGCTGAAAAAACTGACCtgtg                >  1:263057/1‑150 (MQ=255)
                                                                                                                                       cGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGACCtgtg                <  2:918241/150‑1 (MQ=255)
                                                                                                                                            tCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGACCTGTGCCATc           >  1:612820/1‑150 (MQ=255)
                                                                                                                                               gCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGACCTGTGCCATCAGc        <  2:509870/150‑1 (MQ=255)
                                                                                                                                               gCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGACCTGTGCCATCAGc        <  2:509900/150‑1 (MQ=255)
                                                                                                                                                cATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGACCTGTGCCATCAGcc       >  1:667189/1‑150 (MQ=255)
                                                                                                                                                 aTGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGACCTGTGCCATCAGCCa      <  1:702776/150‑1 (MQ=255)
                                                                                                                                                     gcACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGACCTGTGCCATCAGCCAGaca  >  2:1354644/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCTCTCGCCGGTGGTCTGCTGGCTGTTCCCCTGCAGATCAAAAACATTCAGCAGAACGGGATCACACCAGAAGATATCAGCAAGGAAATTAACGGGTATTACTGCTGTTTTTATACTGCCTTTTTCCTTATGGCGTGTTCTGCATGCGCACCATTGATCGCGTTACAGTACGATATTTCACCGTCACTGATGTGGTGGGGCGGGTTGTTGTACTGGCTGGCTGCATTAGTGACGCTGCTATGGGCGGCCAGCCAGATCCAGGCGCTGAAAAAACTGACCTGTGCCATCAGCCAGACA  >  NC_000913/2075083‑2075381

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